The aim of this paper is to put forward a design model for multilayer free damping structures. It sets up a mathematical model and deduces the formula for its structural loss factor η and analyzes the change rules of η along with the change rate of the elastic modulus ratio q1, the change rate of the loss factors of damping materials q2, and the change rate of the layer thickness ratio q3 under the condition with the layer thickness ratio h2=1,3,5,10 by software MATLAB. Based on three specific damping structures, the mathematical model is verified through ABAQUS. With the given structural loss factor (η≥2) and the layer number (n=3,4,5,6), 34 kinds of multilayer free damping structures are then presented. The study is meant to provide a more flexible and more diverse design solution for multilayer free damping structures.
1. Introduction
The damping structure is the structure with damping viscoelastic materials attached to mechanical parts to consume vibrational energy. When the machine vibrates, the damping layer bends and vibrates. Consequently, the tension stress and strain are generated within the damping material to dissipate the mechanical energy to achieve the vibration damping effect. Damping structures play a key role in vibration control, shock absorption, and noise reduction. They are commonly used in many domains such as aeronautics, aerospace, and mechanical and civil engineering. So many research works have been done on damping structures.
Moita et al. [1] developed a finite element model for vibration analysis of active-passive damped multilayer sandwich plates with a viscoelastic core sandwiched between elastic layers, including piezoelectric layers. Jin et al. [2] presented a unified formulation for vibration and damping analysis of a sandwich beam made up of laminated composite face sheets and a viscoelastic core with arbitrary lay-ups and general boundary conditions. Li and Narita [3] analyzed and did the optimal design for the damping loss factor of laminated plates under general edge conditions. Berthelot et al. [4] developed a synthesis of damping analysis of laminate materials and laminates with interleaved viscoelastic layers and sandwich materials. Araújo et al. [5] formulated a finite element model by using a mixed layer-wise approach for anisotropic laminated plates with a viscoelastic core. They considered two different deformation theories, and the optimal design and parameter estimation were also discussed. Akoussan et al. [6] proposed a high order continuous sensitivity analysis of the damping properties of viscoelastic composite plates according to their layer thicknesses. Yang et al. [7] presented an accurate solution approach based on the first-order shear deformation theory for the free vibration and damping analysis of thick sandwich cylindrical shells with a viscoelastic core under arbitrary boundary conditions. Chen et al. [8] presented an order-reduction-iteration approach for vibration analysis of viscoelastically damped sandwiches.
In light of earlier studies, it appears that great attention has been paid to investigations pertaining to the sandwich damping structure with viscoelastic materials, namely, the damping structure with the constraint layer. For that structure in vibration, shear deformation occurs between the constraint layer and the damping layer, and the vibration and deformation of the damping layer are limited and affected by the constraint layer, so the mechanical energy of vibration cannot be fully converted into the material internal energy, which influences the vibration damping effect, whereas the multilayer free damping structure in this paper is the structure composed of the base layer and the added multilayers of viscoelastic damping materials with different properties. During vibration, there is no shear deformation among those layers, and the vibration and deformation of each damping layer are not limited or affected by the upper layer, and then more mechanical energy of vibration can be converted into the material internal energy, which contributes to a better vibration damping effect. Moreover, for the multilayer free damping structure, we can choose different damping layer materials according to the characteristics of the base layer and adjust the thickness and the number of damping layers to achieve the desired vibration damping effect.
The following sections will introduce the mathematical model for the damping structure, the derivation of the structural loss factor η, the change rules of η along with the change rate of the elastic modulus ratio q1, the change rate of the loss factors of damping materials q2, and the change rate of the layer thickness ratio q3 and the mathematical model verification. Finally, with the given structural loss factor (η≥2) and the layer number (n=3,4,5,6), 34 kinds of multilayer free damping structures are provided.
2. Multilayer Free Damping Structure Model
When the multilayer free damping beam is bending (as shown in Figure 1), the elongation ε of both the basic layer and the multilayer free damping layers along the longitudinal direction varies linearly with y.(1)ε=ydθdx=yiω∂ω∂x.
The vibration parameter description of the bending beam of the multilayer free damping structure. Note. ω is the angular velocity; v is the velocity of transverse vibration; θ is the displacement of the bending angle; P is the transverse force; m is the mass per unit length of beam; M is the bending moment; B- is the composite bending stiffness.
The normal tensile stress along the x direction is expressed as(2)σn=Enεn.n=1,2,… represents the basic layer and the multilayer free damping layers.
Plugging formula (1) into formula (2), the following is gained: (3)σn=Eniωy∂ω∂x.
In pure bending, the force acted on the cross section of the beam along the x direction (see Figure 2) is equal to zero. (4)∫-H1-ζ∑a=2nHa+ζσdy=0,where H1 stands for the thickness of the basic layer and Ha (a=2,…,n) is the thickness of the damping layers. ζ is the distance from the center line of the barycenter of the damping beam to the basic layer surface.
The sectional view of the multilayer free damping structure.
The curvature (1/iω)·(∂ω/∂x) is a constant. If the above formula is equal to zero, then(5)∫-H1-ζ∑a=2nHa+ζEydy=0.
Formula (5) can be expressed as(6)∑k=2n-1∫∑b=2kHb+ζ∑a=2k+1Ha+ζE-k+1ydy+∫ζH2+ζE-2ydy+∫-H1-ζζE1ydy=0.
By formula (6), ζ can be expressed as(7)ζ=12E1H12-E-2H22-∑k=2n-1E-k+1Hk+12∑a=2kHa+Hk+1E1H1+∑a=2nE-aHa.
After ζ is determined, the bending moment M can be given as(8)M=∫-H1-ζ∑a=2nHa+ζyσdy.
Plugging formula (9) into (8), formula (10) can be gained:(9)M=B-iω∂ω∂x(10)B-=∫-H1-ζ∑a=2nHa+ζEy2dy=E1H133∑k=2n-1e-k+1hk+13∑a=2khaha+hk+1+6∑a=2,b=2,a≠bkhahb+hk+12+3ζhk+1+2∑a=2kha+3ζ2+e-2h23+3h22ζH1+3h2ζ2H12+1-3ζH1+3ζ2H12.
Among them, e-a=E-a/E1=(Ea/E1)1+iβa=ea1+iβa and ha=(Ha/H1)(a=2,…,n).
Formula (7) can also be expressed as (11)ζ=H121-e-2h22-∑k=2n-1e-k+1hk+12∑a=2kha+hk+11+∑a=2ne-aha.
Plugging formula (11) into (10), the following can be given as(12)B-=∫-H1-ζ∑a=2nHa+ζEy2dy=B1∑k=2n-1e-k+1hk+112∑a=2khaha+hk+1+24∑a=2,b=2,a≠bkhahb+4hk+12+6H1chk+1+2∑a=2kha+3H12c2+1+2e-22h2+3h22+2h23+e-22h241+e-2h2.
Among them, (13)c=1-e-2h22-∑k=2n-1e-k+1hk+12∑a=2kha+hk+11+∑a=2ne-aha,where B- is the composite bending stiffness of the beam of the multilayer free damping structure and B1 is the bending stiffness of the beam without damping layers.
Compared with the loss factor βa of the viscoelastic damping materials of the damping layers, the loss factor β1 of the basic layer is so small that it can be ignored. And also B-=B1+iη, e-a=ea1+iβa, a=2,…,n, and η is the structural loss factor of the damping structure. So the following can be gained by plugging them into formula (12):(14)BB11+iη=∑k=2n-1ek+11+iβk+1hk+112∑a=2khaha+hk+1+24∑a=2,b=2,a≠bkhahb+4hk+12+6H1hk+1+2∑a=2kha1-e21+iβ2h22-∑k=2n-1ek+11+iβk+1hk+12∑a=2kha+hk+11+∑a=2nea1+iβaha+3H121-e21+iβ2h22-∑k=2n-1ek+11+iβk+1hk+12∑a=2kha+hk+11+∑a=2nea1+iβaha2+1+2e21+iβ22h2+3h22+2h23+e221+iβ22h241+e21+iβ2h2.
Omitting βaβb, the real component can be expressed as(15)BB1=∑k=2n-1ek+1hk+112∑a=2khaha+hk+1+24∑a=2,b=2,a≠bkhahb+4hk+12+3p1H1p222hk+1+4∑a=2kha+H1p1p22+1+2e22h2+3h22+2h23+e22h241+e2h2.
Among them, p1=p21-e2h22-∑k=2n-1ek+1hk+12∑a=2kha+hk+1 and p2=1+∑a=2neaha.
Omitting βaβb, the imaginary component can be expressed as(16)iBB1η=i∑k=2n-1ek+1hk+1βk+112∑a=2khaha+hk+1+24∑a=2,b=2,a≠bkhahb+4hk+12+3p1H1p222hk+1+4∑a=2kha+H1p1p22-6H1p3p22∑k=2n-1ek+1hk+1hk+1+2∑a=2kha+H1p1p22+iβ2e2h23+6h2+4h22+2e2h23+e22h241+e2h22.
Among them,(17)p3=p2e2β2h22+∑k=2n-1ek+1hk+1βk+12∑a=2kha+hk+1+∑a=2neaβa1-e2h22-∑k=2n-1ek+1hk+12∑a=2kha+hk+1.
Plugging formula (15) into (16) and omitting eaeb, the formula for η can be expressed as(18)η=6H1p6+H1p4+e2h2β23+6h2+4h226H1p7+2+2e2h22+3h2+2h22.
Here a,b=2,3,…,n.(19)p4=e2β2h22+∑a=2neaβa+p6.p5=hk+1+2∑a=2khap6=∑k=2n-1ek+1hk+1βk+1p5p7=∑k=2n-1ek+1hk+1p5
3. Analysis of Influencing Factors of the Structural Loss Factor η
From formula (18), it is known that there are several factors to influence the structural loss factor η. They are the thickness of the basic layer H1, the elastic modulus ratio of the damping material of each layer ea=Ea/E1, the layer thickness ratio of each damping layer ha=Ha/H1, the loss factor of the damping material of each layer βa. In this damping structure, from the 2nd to the nth damping layer, the elastic moduli of the damping materials have the relation: E2>E3>⋯>En, the layer thicknesses have the relation: H2>H3>⋯>Hn, and the loss factors of the damping materials have the relation: β2>β3>⋯>βn, so there are the relations: e2>e3>⋯>en and h2>h3>⋯>hn. It is assumed that the elastic modulus ratio, the loss factor, and the layer thickness ratio of the damping layers have the relations: en=q1n-2e2, βn=q2n-2β2, and hn=q3n-2h2. Those are then plugged into formula (19), and the following formulas are gained:(20)p4=e2β2h22+q1q221-q1q2n-21-q1q2+p6p5=h2q3k+2q321-q3k-21-q3p6=e2β2h2q1q2q3221-q1q2q32n-21-q1q2q32+21-q31-q1q2q3n-21-q1q2q3-q31-q1n-2q2n-2q32n-51-q1q2q32p7=e2h2q1q3221-q1q32n-21-q1q32+21-q31-q1q3n-21-q1q3-q31-q1n-2q32n-51-q1q32.
Plugging formulas (20) into (18), the formula for the structural loss factor η can be worked out. Then it can be analyzed and calculated by software MATLAB.
Figure 3 shows the change rule of the structural loss factor η along with the change rate of the elastic modulus q1 under the condition with the layer thickness ratio h2=1,3,5,10 and the layer number n=3~10. When h2=1, η is increasing with q1, and the more the damping layers are, the more rapidly η increases. When h2=3, η is increasing with q1, and the more the damping layers, the more slowly η increases. When h2=5, if q1<0.5, η is increasing with q1; if q1>0.5 and n≤5, η is increasing with q1, and the more the damping layers, the more slowly η increases; n>5, η is decreasing with q1, and the more the damping layers, the more rapidly η decreases. When h2=10, η is decreasing with q1, and the more the damping layers, the more rapidly η decreases.
The change rule of η along with q1.
h2=1
h2=3
h2=5
h2=10
Figure 4 shows the change rule of the structural loss factor η along with the change rate of the loss factors of damping materials q2 under the same condition. When h2=1, η is increasing with q2, and the more the damping layers, the more rapidly η increases. When h2=3,5,10, η is increasing with q2, and, for those damping structures with n≤5, η increases obviously faster.
The change rule of η along with q2.
h2=1
h2=3
h2=5
h2=10
Figure 5 shows the change rule of the structural loss factor η along with the change rate of the layer thickness ratio q3 under the same condition. When h2=1, if 0.1<q3<0.7, η is decreasing slowly with q3; if 0.7<q3<0.9, η is increasing rapidly with q3. When h2=3,5,10, if q3<0.5, η is decreasing slowly with q3; if q3>0.5, η is decreasing rapidly with q3.
The change rule of η along with q3.
h2=1
h2=3
h2=5
h2=10
4. Mathematical Model Verification of the Multilayer Free Damping Structure
In Section 3, the change rules of η along with q1, q2, and q3 are analyzed based on the mathematical model which has to be verified. Obviously, since there are a large number of design solutions, it is impossible to evaluate all of them. Therefore, this section will only evaluate three specific damping structures by ABAQUS to verify the model. For the three structures, the length of the basic layer is given as 50 cm, the width as 10 cm, and the thickness as 1 cm. Its material density is given as 7800 kg/m3. The material parameter uses isotropic material with a Young’s modulus of 210 GPa and a Poisson’s ratio of 0.3. And they are all with the layer thickness ratio h2=5 and the layer number n=5 (as is shown in Figure 6). For the first damping structure, q2 and q3 are defined as fixed values and q1 is defined as variables. For the second damping structure, q1 and q3 are defined as fixed values and q2 is defined as variables. For the third damping structure, q1 and q2 are defined as fixed values and q3 is defined as variables.
The damping structure.
4.1. Finite Element Analysis for the First Damping Structure
The results of the finite element analysis of the first damping structure are shown in Figure 7. We set q2 and q3 as 1.1 and 0.9, respectively. If q1=0.9, its amplitude of vibration is reduced by 27% and its strain energy is reduced by 28.6%. If q1=0.7, its amplitude of vibration is reduced by 31.3% and its strain energy is reduced by 83.3%. The damping effect of the structure with q1=0.7 is better than the one with q1=0.9. That means the structural loss factor η of the structure with q1=0.7 is greater than the one with q1=0.9. These results are consistent with the previous results obtained from MATLAB in Figure 3(c).
The amplitude and strain energy images of the first damping structure.
q1=0.9 amplitude of vibration
q1=0.7 amplitude of vibration
q1=0.9 strain energy
q1=0.7 strain energy
4.2. Finite Element Analysis for the Second Damping Structure
The results of the finite element analysis of the second damping structure are shown in Figure 8. We set both q1 and q3 as 0.9. If q2=1.1, its amplitude of vibration is reduced by 25% and its strain energy is reduced by 20%. If q2=1.3, its amplitude of vibration is reduced by 50% and its strain energy is also reduced by 50%. The damping effect of the structure with q2=1.3 is better than the one with q2=1.1. That means the structural loss factor η of the structure with q2=1.3 is greater than the one with q2=1.1. These results are consistent with the previous results obtained from MATLAB in Figure 4(c).
The amplitude and strain energy images of the second damping structure.
q2=1.1 amplitude of vibration
q2=1.3 amplitude of vibration
q2=1.1 strain energy
q2=1.3 strain energy
4.3. Finite Element Analysis for the Third Damping Structure
The results of the finite element analysis of the third damping structure are shown in Figure 9. We set q1 and q2 as 0.9 and 1.1, respectively. If q3=0.9, its amplitude of vibration is reduced by 29.7% and its strain energy is reduced by 40%. If q3=0.7, its amplitude of vibration is reduced by 35.5% and its strain energy is reduced by 50%. The damping effect of the structure with q3=0.7 is better than the one with q3=0.9. That means the structural loss factor η of the structure with q3=0.7 is greater than the one with q3=0.9. These results are consistent with the previous results obtained from MATLAB in Figure 5(c).
The amplitude and strain energy images of the third damping structure.
q3=0.9 amplitude of vibration
q3=0.7 amplitude of vibration
q3=0.9 strain energy
q3=0.7 strain energy
With the three damping structures, the mathematical model is validated by the comparisons between their MATLAB results and ABAQUS results, which are found to have a good agreement.
5. Application Selection of the Multilayer Free Damping Structure
Let us suppose that the structural loss factor η of the damping structure is equal to or greater than 2, the layer thickness ratio h2 is 1, 3, 5, and 10, and the layer number n is 3, 4, 5, and 6. 34 kinds of design solutions are obtained by combining different damping materials of different thicknesses (as is shown in Table 1).
Multilayer free damping structure parameters.
Layers n
Layer thickness ratio h2
Data range
3
1
q1=0.9, q2=1.1, 0.1<q3<0.9
3
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.215<q2<1.3
5
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.19<q2<1.3
10
q1=0.9, q3=0.9, 1.175<q2<1.3
q2=1.1, q3=0.9, 0.1<q1<0.9
q1=0.9, q2=1.1, 0.1<q3<0.9
4
1
q1=0.9, q2=1.1, 0.1<q3<0.9
3
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.205<q2<1.3
5
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.185<q2<1.3
10
q1=0.9, q3=0.9, 1.17<q2<1.3
q2=1.1, q3=0.9, 0.1<q1<0.9
q1=0.9, q2=1.1, 0.1<q3<0.9
5
1
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.29<q2<1.3
3
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.19<q2<1.3
5
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.17<q2<1.3
10
q1=0.9, q3=0.9, 1.16<q2<1.3
q2=1.1, q3=0.9, 0.1<q1<0.9
q1=0.9, q2=1.1, 0.1<q3<0.9
6
1
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.28<q2<1.3
3
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.17<q2<1.3
5
q1=0.9, q2=1.1, 0.1<q3<0.9
q1=0.9, q3=0.9, 1.15<q2<1.3
10
q1=0.9, q3=0.9, 1.155<q2<1.3
q2=1.1, q3=0.9, 0.1<q1<0.9
q1=0.9, q2=1.1, 0.1<q3<0.9
In Table 1, for example, with the layer n=3 and the layer thickness ratio h2=10, three sets of data ranges of q1, q2, and q3 are gained. In Figure 3(d), when q2=1.1 and q3=0.9, q1 should meet the condition 0.1<q1<0.9. In Figure 4(d), when q1=0.9 and q3=0.9, q2 should meet the condition 1.175<q2<1.3. In Figure 5(d), when q1=0.9 and q2=1.1, q3 should meet the condition 0.1<q3<0.9. In line with the practical request, the relative parameters, such as the elastic modulus, the loss factors of damping materials, and the layer thickness, are further determined to achieve the desired vibration damping effect.
6. Conclusion
In this paper, a mathematical model for the multilayer free damping structure is presented and verified. And the relations of its structural loss factor η with the change rate of the elastic modulus ratio q1, the change rate of the loss factors of damping materials q2, and the change rate of the layer thickness ratio q3 are studied by MATLAB. Then 34 design solutions of the damping structures are listed out for choice.
In practical application, only if some parameters (such as the structural loss factor, the base layer material property) are known can some other parameters (such as the damping material, the layer thickness, and the layer number) be set and then be compared. Finally, an optimal combination is selected to achieve the requested vibration damping effect.
As future work, it is necessary to improve the theory by taking into account the temperature and frequency which are the most important environmental factors affecting the dynamic properties of damping materials. Moreover, for optimization, an intensive study is needed into the relationships between the total layer thickness, the layer number, and the structural loss factor to obtain a more careful design solution.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors wish to express their sincere thanks to Liu Shizhong, Yan Bijuan, and Song Yong of Damping Vibration Attenuation Laboratory of Taiyuan University of Science and Technology for their general guidance and technical advice.
MoitaJ. S.AraújoA. L.MartinsP.Mota SoaresC. M.Mota SoaresC. A.A finite element model for the analysis of viscoelastic sandwich structures20118921-22187418812-s2.0-8005292586110.1016/j.compstruc.2011.05.008JinG.YangC.LiuZ.Vibration and damping analysis of sandwich viscoelastic-core beam using Reddy's higher-order theory20161403904092-s2.0-8496046021310.1016/j.compstruct.2016.01.017LiJ.NaritaY.Analysis and optimal design for the damping property of laminated viscoelastic plates under general edge conditions20134519729802-s2.0-8486950474910.1016/j.compositesb.2012.09.014BerthelotJ.-M.AssararM.SefraniY.MahiA. E.Damping analysis of composite materials and structures20088531892042-s2.0-4434909308610.1016/j.compstruct.2007.10.024AraújoA. L.MartinsP.SoaresC. M. M.SoaresC. A. M.HerskovitsJ.Damping optimisation of hybrid active-passive sandwich composite structures2012466974AkoussanK.BoudaoudH.El-MostafaD.CarreraE.Vibration modeling of multilayer composite structures with viscoelastic layers2015221-21361492-s2.0-8490766329910.1080/15376494.2014.907951YangC.JinG.LiuZ.WangX.MiaoX.Vibration and damping analysis of thick sandwich cylindrical shells with a viscoelastic core under arbitrary boundary conditions2015921621772-s2.0-8492051737410.1016/j.ijmecsci.2014.12.003ChenX.ChenH. L.HuX. L.Damping predication of sandwich structures by order-reduction-iteration approach199922258038122-s2.0-001245052710.1006/jsvi.1998.2131