Research on the Multilayer Free Damping Structure Design

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 dampingmaterials 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 throughABAQUS.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.


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 firstorder 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  1 , the change rate of the loss factors of damping materials  2 , and the change rate of the layer thickness ratio  3 and the mathematical model verification.Finally, with the given structural loss factor ( ≥ 2) and the layer number ( = 3, 4, 5, 6), 34 kinds of multilayer free damping structures are provided.

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 .
The normal tensile stress along the  direction is expressed as = 1, 2, . . .represents the basic layer and the multilayer free damping layers.Plugging formula (1) into formula (2), the following is gained: In pure bending, the force acted on the cross section of the beam along the  direction (see Figure 2) is equal to zero.
where  1 stands for the thickness of the basic layer and   ( = 2, . . ., ) 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 curvature (1/) ⋅ (/) is a constant.If the above formula is equal to zero, then Formula ( 5) can be expressed as By formula (6),  can be expressed as After  is determined, the bending moment  can be given as Plugging formula (9) into (8), formula (10) can be gained: Among them,   =   / 1 = (  / 1 )(1+  ) =   (1+  ) and ℎ  = (  / 1 ) ( = 2, . . ., ).Formula (7) can also be expressed as Plugging formula (11) into (10), the following can be given as Among them, where  is the composite bending stiffness of the beam of the multilayer free damping structure and  1 is the bending stiffness of the beam without damping layers.
Compared with the loss factor   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  = (1 + ),   =   (1 +   ),  = 2, . . ., , and  is the structural loss factor of the damping structure.So the following can be gained by plugging them into formula (12): Omitting     , the real component can be expressed as (15) Among them, Omitting     , the imaginary component can be expressed as Among them, Plugging formula (15) into ( 16) and omitting     , the formula for  can be expressed as Here ,  = 2, 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  1 , the elastic modulus ratio of the damping material of each layer   =   / 1 , the layer thickness ratio of each damping layer ℎ  =   / 1 , the loss factor of the damping material of each layer   .In this damping structure, from the 2nd to the th damping layer, the elastic moduli of the damping materials have the relation:  2 >  3 > ⋅ ⋅ ⋅ >   , the layer thicknesses have the relation:  2 >  3 > ⋅ ⋅ ⋅ >   , and the loss factors of the damping materials have the relation:  2 >  3 > ⋅ ⋅ ⋅ >   , so there are the relations: It is assumed that the elastic modulus ratio, the loss factor, and the layer thickness ratio of the damping layers have the relations: 2  2 , and ℎ  =  −2 3 ℎ 2 .Those are then plugged into formula (19), and the following formulas are gained: 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  1 under the condition with the layer thickness ratio ℎ 2 = 1, 3, 5, 10 and the layer number  = 3∼10.When ℎ 2 = 1,  is increasing with Shock and Vibration   1 , and the more the damping layers are, the more rapidly  increases.When ℎ 2 = 3,  is increasing with  1 , and the more the damping layers, the more slowly  increases.When ℎ 2 = 5, if  1 < 0.5,  is increasing with  1 ; if  1 > 0.5 and  ≤ 5,  is increasing with  1 , and the more the damping layers, the more slowly  increases;  > 5,  is decreasing with  1 , and the more the damping layers, the more rapidly  decreases.When ℎ 2 = 10,  is decreasing with  1 , and the more the damping layers, the more rapidly  decreases.
Figure 4 shows the change rule of the structural loss factor  along with the change rate of the loss factors of damping materials  2 under the same condition.When ℎ 2 = 1,  is increasing with  2 , and the more the damping layers, the more rapidly  increases.When ℎ 2 = 3, 5, 10,  is increasing with  2 , and, for those damping structures with  ≤ 5,  increases obviously faster.

Mathematical Model Verification of the Multilayer Free Damping Structure
In Section 3, the change rules of  along with  1 ,  2 , and  3 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/m 3 .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 ℎ 2 = 5 and the layer number  = 5 (as is shown in Figure 6).For the first damping structure,  2 and  3 are defined as fixed values and  1 is defined as variables.For the second damping structure,  1 and  3 are defined as fixed values and  2 is defined as variables.For the third damping structure,  1 and  2 are defined as fixed values and  3 is defined as variables.

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  2 and  3 as 1.1 and 0.9, respectively.If  1 = 0.9, its amplitude of vibration is reduced by 27% and its strain energy is reduced by 28.6%.If  1 = 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  1 = 0.7 is better than the one with  1 = 0.9.That means the structural loss factor  of the structure with  1 = 0.7 is greater than the one with  1 = 0.9.These results are consistent with the previous results obtained from MATLAB in Figure 3(c). 3 as 0.9.If  2 = 1.1, its amplitude of vibration is reduced by 25% and its strain energy is reduced by 20%.If  2 = 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  2 = 1.3 is better than the one with  2 = 1.1.

Finite Element Analysis for the
That means the structural loss factor  of the structure with  2 = 1.3 is greater than the one with  2 =1.1.These results are consistent with the previous results obtained from MATLAB in Figure 4(c).

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  1 and  2 as 0.9 and 1.1, respectively.If  3 = 0.9, its amplitude of vibration is reduced by 29.7% and its strain energy is reduced by 40%.
If  3 = 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  3 = 0.7 is better than the one with  3 = 0.9.That means the structural loss factor  of the structure with  3 = 0.7 is greater than the one with  3 = 0.9.These results are consistent with the previous results obtained from MATLAB in Figure 5(c).
With the three damping structures, the mathematical model is validated by the comparisons between their MAT-LAB results and ABAQUS results, which are found to have a good agreement.

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 ℎ 2 is 1, 3, 5, and 10, and the layer number  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).

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
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.

Figure 1 :
Figure1: 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;  is the transverse force;  is the mass per unit length of beam;  is the bending moment;  is the composite bending stiffness.

Figure 2 :
Figure 2: The sectional view of the multilayer free damping structure.

Figure 3 :
Figure 3: The change rule of  along with  1 .
rat e of ma ter ial los s fac tor La ye r qu an tit y ng e ra te of m at er ia l lo ss fa ct or L ay er q u an ti ty Structural rat e of ma ter ial los s fac tor L a y e r q u a n ti ty Structural rat e of m ate ria l los s fac tor L ay er q u an ti ty Structural damping loss factor(d) ℎ 2 = 10

Figure 4 :
Figure 4: The change rule of  along with  2 .
rat e of lay er thi ckn ess rat io La ye r qu an tit y rat e of lay er th ick ne ss rat io L ay er q u an ti ty Structural rat e of lay er thi ck ne ss rat io L ay er q u an ti ty Structural rat e of lay er thi ckn ess rat io La ye r qu an ti ty Structural damping loss factor(d) ℎ 2 = 10

Figure 5 :
Figure 5: The change rule of  along with  3 .

Figure 8 :
Figure 8: The amplitude and strain energy images of the second damping structure.

Figure 9 :
Figure 9: The amplitude and strain energy images of the third damping structure.

Table 1 :
Multilayer free damping structure parameters.