Research Article  Open Access
Hui Li, Yi Niu, Chao Mu, Bangchun Wen, "Identification of Loss Factor of FiberReinforced Composite Based on Complex Modulus Method", Shock and Vibration, vol. 2017, Article ID 6395739, 13 pages, 2017. https://doi.org/10.1155/2017/6395739
Identification of Loss Factor of FiberReinforced Composite Based on Complex Modulus Method
Abstract
The identification of the loss factor of fiberreinforced composite based on complex modulus method is presented. Firstly, the damping model of fiberreinforced composite plate is established, and the relation between each loss factor and modal damping ratio is deduced based on the complex modulus method. Then, the least square relative error function is formed by using the modal damping ratio obtained in the experimental test, and the appropriate stepsize is selected in the range of 0~10% to calculate the loss factor. Next, the identification procedure of loss factor of such composite material is summarized, and the corresponding identification procedure is realized based on selfdesigned MATLAB program. Finally, TC300 carbon/epoxy composite plate is taken as an example to carry out a case study, and its loss factors along the longitudinal, transverse, and shear direction are identified by the complex modulus method. By comparing the measured damping results obtained in this paper and the calculated damping results based on the AdamsBacon model with the same loss factor, it is found that the corresponding maximum deviation between them is less than 15%, so the correctness of such identification method has been verified indirectly, which can be used to identify loss factor of fiberreinforced composite with high precision and efficiency.
1. Introduction
Fiberreinforced composite has excellent mechanical properties, good thermal stability, and capability for weight reduction, which is widely used in the area of aeronautics, astronautics, automotive, naval vessel, and weapon industry [1, 2]. Currently, there are a large number of such composite thin walled structures, such as solar panels, aircraft engine fan blades, and large wind turbine blades. Loss factor is a key parameter for such type of composite, which has significant influence on vibration control, noise suppression, and fatigue life of composite structures, and usually it can be divided into longitudinal, transverse, and shear direction due to its anisotropy property. But, because of the complexity of the microstructure and different dynamic loads on the composite structures [3, 4], for instance, the loss factor is not only closely related to fiber orientation, fiber/matrix interface conditions, material preparation process, environmental temperature, and and other factors, but also closely related to the external excitation frequency and amplitude [5, 6]. Therefore, it is of great importance to study the damping property of fiberreinforced composite material and its structure, especially the corresponding damping analysis methods and identification techniques.
At present, great progress has been made in the identification of the loss factor of fiberreinforced composites. Adams and Bacon [7] firstly established the macroscopic model of fiberreinforced composite plate, which was called AdamsBacon damping model. They considered that the energy consumed by each thin layer in the composite plate was the sum of its longitudinal tensile, transverse tensile, and shear stresses. Therefore, the damping capacity of the thin layer can be defined as the ratio of the dissipated energy to the stored strain energy. Rao and He [8] presented closedform solutions for resonance frequencies and modal loss factors of laminated composite beams with multiple viscoelastic layers by the energy formulation and Ritz method. Rikards et al. [9] established the damping model of the composite plate by introducing the complex elastic moduli, and the frequency dependence of the loss factor in the longitudinal, transverse, and shear directions of composite plates was studied, but they only obtained bending loss factor at different frequencies and also did not give the detailed identification procedure of loss factors , , and . De Visscher et al. [10] proposed a mixed numericalexperimental method for the identification of the material damping properties of fiberreinforced polymer composites, and the complex elastic moduli were introduced, so that the modal strain energy can be expressed as the sum of four partial energies which were longitudinal tensile, transverse tensile, shear stresses, and thickness direction. Dalenbring [11] proposed a threedimensional material damping estimation methodology for planar isotropic material symmetry by using a constitutive viscoelastic vibration model, and an experimental study of the loss factors of the fiberreinforced composite was carried out. Melo and Radford [12] established an analytical model based on only five independent dynamic stiffness parameters and three independent damping loss factors, and, by using this model, a method was developed based on energy equations. Kostopoulos and Korontzis [13] introduced a mixed analyticalexperimental approach for the measurement of the loss factors with frequency dependence of the single lamina, and the effects of experimental errors on the evaluation of the viscoelastic properties of the single lamina were discussed. Chandra et al. [14] measured all the six loss factors for glass fiberreinforced composite (three in normal and three in shear) by the free decay method in the vacuum chamber, and different types of specimen such as beam, tubular, and cuboidal in shape were tested under different loading conditions. Berthelot and Sefrani [15] presented an extensive analysis of the damping of unidirectional fiber composite as function of frequency and fiber orientation, and the longitudinal loss factor can be identified for 0° fiber orientation while the transverse loss factor can be identified for 90° fiber orientation. Berthelot [16] also presented an analysis of the damping of unidirectional fiber composites and different laminates, and damping modeling of unidirectional composites and laminates was developed using the Ritz method for describing the flexural vibrations of beams or plates. Kumar et al. [17] proposed threephase mathematical model for the evaluation of damping incorporating the effect of fiber packing, and the longitudinal loss factor, transverse loss factor, transverse shear loss factor, and longitudinal shear loss factor were predicted under different interphase volume fractions. Shokrieh and Najafi [18] measured the specific damping capacity and loss factors of laminated polymer matrix composite under the frequency of 0.2 Hz–4 Hz, and its viscoelastic behavior was characterized by using a modified classical lamination theory. Matter et al. [19] presented a mixed numericalexperimental identification method for estimating the loss factor of composite plates and shells, and the natural frequencies, modal damping factors, and mode shapes of the specimen were measured with an optimized contactfree experimental setup, which used a loudspeaker for exciting the structure and a scanning laser interferometer for measuring the dynamic response. Devalve and Pitchumani [20] investigated the damping effects of carbon nanotubes embedded in the matrix of fiberreinforced composite material, and dynamic mechanical analysis (DMA) system was employed in the dualcantilever mode of the composite beam to measure its material loss factor.
Although many researchers have studied the identification of the loss factor of fiberreinforced composite plates, most of them usually regard the elasticity modulus as a real number, and few literatures have studied the damping characteristics of composite structures with complex modulus method, which regard the elasticity modulus as an imaginary number. At present, the complex modulus theory has been widely used to describe the damping properties of viscoelastic material and structure [21, 22], and it has a series of advantages, such as simple theoretical principle, high calculation accuracy and low cost of time expenditure. Therefore, in this research, such theory method is extended to identify the loss factor of fiberreinforced composites, which is a new try to find out the identification solution of damping property of fiberreinforced composite material.
In this research, the identification of loss factor of fiberreinforced composite based on complex modulus method is presented. Firstly, in Section 2, the damping model of fiberreinforced composite plate is established as an example, and the total strain energy and dissipation energy expressions are deduced, so that the relation between each loss factor and the modal damping ratio can be determined based on the complex modulus method. Besides, the least square relative error function is formed by using the modal damping ratio obtained by experimental test, and the appropriate stepsize is selected in the range of 0~10% to calculate the loss factor. Then, in Section 3, the identification procedure of loss factor of such composite laminates is summarized, and the corresponding identification procedure is realized based on selfdesigned MATLAB program. Finally, in Section 4, by taking the TC300 carbon/epoxy composite plate as the research object, the damping parameters are measured in experiment, and the loss factor is also identified by the complex modulus method. And, by comparing the measured damping results obtained in this paper and the calculated damping results based on the AdamsBacon model with the same loss factor, the correctness of such identification method has been verified indirectly.
2. Theoretical Principle of Identifying Loss Factor of FiberReinforced Composite Based on Complex Modulus Method
In this section, complex modulus method is proposed to identify loss factor of fiberreinforced composites, which is based on the established damping model of composite thin plate, and its theoretical principle is explained in the following.
2.1. Damping Model of Composite Thin Plate
The fiberreinforced composite thin plate, which is made of fiber and matrix material with orthotropic layers, is introduced in Figure 1. Firstly, the coordinate system is set up at the middle surface, in which the length, width, and thickness of the composite plate can be expressed as , , and , while the fiber direction within a layer is defined as from the axis of coordinate system . Besides, each layer of the composite plate is located at and along the axis with equal thickness of . In this theoretical model, “1” represents the direction parallel to the fiber, while “2” is the direction perpendicular to the fiber, and “3” is direction perpendicular to the 12 surface. Assuming and as the complex elastic moduli of the layer parallel and perpendicular to fiber, respectively, as the complex shear modulus in the 12 surface, and as Poisson’s ratios which are induced by the stress in “1” and “2” direction, and as the density of composite plate, the complex elastic modulus in each direction can be respectively expressed aswhere, , , and are the loss factors of fiber longitudinal, transverse, and shear direction.
The concerned fiberreinforced composite plate is symmetrical between the middle surfaces, so its inner and outer displacements are decoupled. According to the classical laminate theory, its displacement field can be expressed aswhere “” is complex, , , and are the displacements of any point of composite plates, and , , and are the displacement in the midplane. Besides, is the frequency and is the time.
In the fiber coordinate system, the constitutive equation of the th layer can be written as where is complex stiffness matrix coefficients of the composite plates with and , whose expression is as follows:
The coefficient of complex stiffness includes the real and the imaginary part, therefore, it can be expressed as follows:where and are the real part and imaginary part of the complex stiffness matrix coefficient.
Based on the assumed displacement field of the classical laminate theory, the normal strain and shear strains and of composite plate are equal to zero; that is, . Considering the relationship between strain and displacement, the strain of any point of composite plates can be obtained as
The strain energy stored in the laminate structure of the th layer can be expressed aswhere, , , and are strain energy of the layer fiber of reinforced plate in the longitudinal, transverse, and thickness direction of fiber, which can be expressed as
Substitute (3) into relation (8) which yieldswhere
The strains , , and are related to the strains , , and expressed in the geometric directions according to the strain transformation relation:where and .
Substitution of (11) into relation (10) results in
The total inplane strain energies stored in the layers can thus be obtained by summation of strain energies of all the layers:
The relationship between the total strain energy and strain energy values of fiberreinforced composite plate , , , and can be written as
The total dissipated energy of fiberreinforced composite plates is expressed aswhere the expression is as follows:
According to the definition of the specific damping capacity of fiberreinforced composite plate, it can be expressed as
Generally, the relationship between the specific damping capacity and structural damping ratio can be expressed as
By substituting (18) into (17), the relationship between the strain energy and the damping ratio of the fiberreinforced composite plate can be expressed as
In this way, the relationship between the material loss factors , , and and the damping ratio is established, based on (19).
2.2. Theoretical Calculation of Modal Damping Ratio of Composite Thin Plate
We assume the midplane displacement of composite plate is , and it can be expressed aswhere is the excitation frequency and is modal shape function which can be defined as where are the eigenvectors which need to be solved and and are the orthogonal polynomials.
The orthogonal polynomials can be obtained by implementing orthogonalization operation on polynomial function which should satisfy the boundary condition of composite plate, and these polynomials have the following expressions:where and are the coefficient functions and their expressions can be written aswhere is the weighting function and generally . Moreover, and are polynomial functions which can satisfy different boundary conditions, such as the clamped, simply supported, and free boundary, and they can be expressed as
Considering that the studied composite plate is under free boundary condition, the four parameters are given by , , , and . The damping can be neglected, considering that it has little influence on natural frequency; therefore, only the real part is calculated in the stressstrain constitutive equation (3). The kinetic energy of bending vibration of thin plates can be expressed as
Substitute (19) into (15), then the expression of kinetic energy expressed by undetermined parameters can be obtained. Set the harmonic components in kinetic energy to be 1, and the expression of the maximum kinetic energy can be derived as
Similarly, the strain energy in (14) is also expressed by the undetermined parameters . Set the harmonic component of the strain energy to be 1, and the expression of the maximum strain energy can be obtained. Next, the Lagrangian energy function is expressed as
In the same way to solve the inherent characteristics of fiberreinforced composite thin plate, all pending parameters that making the Lagrangian energy function with respect to have the minimum value can be solved; then the following expression is obtained:
Consequently, the result obtained leads to equations in frequency domain which can be written as the matrix form:where and are the stiffness matrix and mass matrix, respectively, and is the generalized displacement vector.
The elements in and are coefficients of all pending parameters. In order to make (32) have nonzero solution or nontrivial solution, the determinant of the coefficient matrix should be equal to zero:
In this way, the natural frequency of each mode can be obtained. In (21), the larger the values of , , the more accurate the calculated results of natural frequency, and can usually satisfy sufficient accuracy. Meanwhile, substitution of the eigenvectors which correspond to each natural frequency into (21) can lead to the modal shapes. Moreover, substitution of the modal shape function into (20) can lead to the midplane displacement of the composite plate.
Substitution of (6) into (13) leads to the relationship expression between midplane displacement and strain energy of the composite plate.
Substitution of (6) into (13) leads to the relationship expression between midplane displacement and dissipated energy of the composite plate.
The strain energy and dissipation energy of the composite plate are obtained by substituting midplane displacement which is related to each natural frequency into (31) and (32). Then, substitute the above results into (19), and we can determine the functional relation between the material loss factors , , and and the modal damping ratio.
2.3. Identification of Loss Factor of Composite Thin Plate Based on Least Square Method
Because the dissipated energy of the composite plate is a function of loss factors , , and , according to the relationship between damping ratio, strain energy, and dissipation energy in (19), it can be found that the modal damping ratio of composite plates is also the function of loss factors. The least square method is used to solve the loss factor because its algorithm is simple and highly efficient. Firstly, the damping ratio of composite plate is obtained by (19) and the experimental test within a certain frequency band, and the relative error function is established aswhere is the number of modes included in the frequency band and is the th difference for damping ratio between theoretical calculation and experimental results, while is the th damping ratio obtained by experimental test.
The value of loss factors of fiberreinforced composite plate is generally less than , and therefore, the material loss factor can be given as
Select the appropriate stepsize (e.g., ) to construct the material loss factor vectors , , and , which can be expressed as where , and .
Then, by iteratively solving the material loss factor in a permutation and combination manner, we can identify the desired loss factor along the longitudinal, transverse, and shear direction of fiberreinforced composite when the relative error function takes the minimum value.
3. Identification Procedure
In Section 2, the identification principle of loss factor of fiberreinforced composite thin plate by using complex modulus method is gradually clarified. In this section, the identification procedure of loss factor of such composite laminates is summarized, and the corresponding identification procedure is realized based on selfdesigned MATLAB program, which can be divided into the following steps.
(1) Input the Geometric Parameters and Material Parameters of the Composite Plate. Firstly, the geometric parameters such as length, width, thickness, and total number of layers of fiberreinforced composite plate are required. Then, the complex elastic modulus, Poisson’s ratio, and density of the composite thin plate are individually input to prepare for the subsequent calculation of strain energy and dissipated energy.
(2) Calculate the Mass Matrix and Complex Stiffness Matrix of Composite Plate. The mass matrix and complex stiffness matrix of composite thin plates are calculated, with considering the influence of fiber orientation based on classical laminated plate theory. The value of the complex stiffness matrix of the th layer of the fiberreinforced composite plate can be obtained by (4) and substitution of it into (3) leads to the constitutive equation of the composite plate.
(3) Solve the Natural Frequencies and Midplane Displacement of Composite Thin Plate by Using the Ritz Method. The natural frequencies can be obtained by solving the characteristic equation of free vibration of fiberreinforced composite plates. The stiffness matrix and mass matrix gained from step (2) can lead to the characteristic equation (29). In the case of nonzero solution of (29), the determinant of the coefficient matrix should satisfy (30). Then the natural frequencies can be obtained by solving (30). The mode shapes can be obtained by substituting the eigenvectors which correspond to each natural frequency into (21). Moreover, substitute the mode shape function into (20), and then the midplane displacement of composite thin plate can be obtained.
(4) Obtain Total Strain Energy, Dissipation Energy, and Function of Damping Ratio and Loss Factor. Substitute the midplane displacement which corresponds to each natural frequency into (31) and (32), and then the strain energy values , , , and and dissipation energy values , , , and can be derived. Next, inputting the above results into (19), the functional relation between the damping ratio of various modes and the loss factors , , and can be obtained.
(5) Identify Loss Factor of Composite Thin Plate by Least Square Method. In view of the fact that the value of the loss factor is less than 10%, the appropriate stepsize (e.g., is selected to construct the material loss factor vectors and identify the loss factor of composite thin plate by least square method based on the damping ratio obtained in a certain frequency band. The value of loss factor is regarded as the best predictor when the least squares error function takes the minimum value, and identification process of the loss factor of fiberreinforced composite thin plate is given in Figure 2.
4. A Study Case
In this section, a TC300 carbon/epoxy composite plate is taken as a research object, as shown in Figure 3. Its damping parameters are measured under free boundary condition in experiment, and the loss factor of such material is also identified by the complex modulus method. The studied TC300 carbon/epoxy composite plate is symmetrically laid, which has total 21 layers with laminate configuration of [(0/90)_{5}/0/(90/0)_{5}], and each layer has the same thickness and fiber volume fraction with longitudinal elastic modulus of GPa, transverse elastic modulus of GPa, shear modulus of GPa, Poisson’s ratio of 0.33, and density of 1780 kg/m^{3}. Besides, the length, width, and thickness of composite plate are 200 × 130 × 2.36 mm.
4.1. Damping Test System and Test Results
The damping test system of composite thin plate is set up, and its connection schematic can be seen in Figure 3. The instruments used in the test are as follows: (I) BK 4517001 lightweight accelerometer; (II) LMS SCADAS 16channel mobile frontend; (III) Dell M7700 notebook computer and LMS Test.Lab 10B software; (IV) PCB 086C01Hammer. After the repeated testing and comparison, the accelerometer is glued firmly at reference point P by super glue 502, and the point P is 40 mm from the left edge and 30 mm from the downside edge of the composite plate. Then, the following setups and parameters are chosen: (I) sampling frequency of 3200 Hz; (II) frequency resolution of 0.125 Hz; (III) forceexponential window for excitation signal and exponential window for response signals. Consequently, we can employ pulse excitation technique to carry out modal test, and, by analyzing the response signal in frequency domain with LMS software, the natural frequency and the damping ratio of each mode can be obtained with the halfpower bandwidth method. The actual photograph of damping test of fiberreinforced composite thin plate can be seen in Figure 4, and Figure 5 also lists the measured frequency response functions at three different measuring points. In order to improve test accuracy of frequency and damping results, the final results are obtained by averaging the results at these measuring points, which are listed in Table 1.

4.2. Identification Process and Results
By using the identification process in Section 2, the natural frequency of composite thin plate can be calculated according to the elastic modulus and other material parameters given by the manufacturer. Table 2 gives the natural frequencies obtained by such method and also gives the corresponding errors compared with the experimental results. Next, the functional relation between the damping ratio and the loss factor is derived based on (19), and loss factor vector is constructed by selecting the appropriate stepsize in the case that loss factor value is less than 10%. Finally, the loss factor of composite thin plate can be identified with least square method. When the least square error function (such as ) takes the minimum, the corresponding loss factor can thus be regarded as the desired result.

Because the loss factor of composite thin plate is closely related to the stepsize, we set the six different stepsizes, namely, 0.001, 0.0005, 0.0002, 0.0001, 0.00005, and 0.00001, to calculate the concerned loss factor, and Table 3 gives the identified loss factor along the longitudinal, transverse, and shear direction of fiberreinforced composite obtained by different stepsizes as well as their time expenditure. It can be found that when the value of stepsize is less than 0.0001, the loss factor results are basically the same with each other, which can meet the requirement of calculation accuracy. However, the time costs of stepsize of 0.00005 and 0.00001 are more than 30 times compared with the corresponding result of stepsize of 0.0001, so we finally choose 0.0001 as the best stepsize to identify each loss factor.

Besides, in order to investigate the sensitivity of the identified loss factors with respect to the calculated maximum error of natural frequencies, we deliberately adjust the thickness, elastic modulus, Poisson’s ratio and so forth, and the identified loss factors with respect to the different maximum error of natural frequencies are listed in Table 4, from which we can see that, with the increase of the calculated maximum error of the first three natural frequencies, the deviations between the identified loss factors would show an increasing tendency. So, in order to ensure that the identified loss factors along the longitudinal, transverse, and shear directions can have better accuracy, we need to improve the calculated accuracy of natural frequencies of composite plate as much as possible, and, in this research, the calculated maximum error of natural frequency based on the established damping model in Section 2 is less than 8%, so it is appropriate to employ this theoretical model to identify the concerned loss factors.

4.3. Verification and Discussion
In order to verify the correctness of the proposed identification method, we also establish the AdamsBacon damping model based on the theoretical principle described in literature [23] to calculate the modal damping ratio of such composite plate, which can be realized by introducing the loss factor results identified in this paper into the algorithm in literature [23]. And in Table 4 the calculated damping results in literature [23] are listed and, for better comparison of damping results, their deviations between the tested damping ratios in this paper and calculated damping ratios in literature [23] are also given in Table 5.

It can be seen from Table 5 that the modal damping ratio results calculated by the AdamsBacon damping model can have a good agreement with the corresponding experimental results, and the maximum deviation of modal damping ratios is less than 15%. Therefore, the correctness and effectiveness of the complex modulus method have been verified indirectly, and the complex modulus method can be applied to the analysis and prediction of damping property of fiberreinforced composites. However, it is still necessary to find out the reasons for the above deviations, which probably result from both theoretical modeling and experimental testing aspects, for example, in the theoretical modeling process, which may contain the following calculation errors: (I) without considering the damping resulting from interfacial friction between layers and the damping resulting from interface defect; (II) without considering the effects of residual stress; (III) without considering the influence of the dispersion of composite material parameters; (IV) without considering the influence of the accumulation of the truncation error and rounding error on the damping results in the calculation process. Meanwhile, in the experimental measurement, the following factors may also lead to some errors: (I) without considering the influence of boundary conditions in the measurement; (II) without considering the influence of additional mass of the lightweight accelerometer; (III) without considering the influence of air damping. Besides, it should be noted that because the premise of identification of loss factor of fiberreinforced composite by complex modulus method is based on the composite plate model, we need to cut a plate test specimen from the fiber composite material, and if the preparation process parameters of such composite material is changed, including fiber angle, layup form, and other variation in configuration, then we need to identify the new loss factors once again, which would somewhat restrict the application of such method in the engineering field. However, the geometrical size of fiberreinforced composite plate indeed will not affect the identification of loss factor results.
5. Conclusion
This research has investigated the identification technique of the loss factor of fiberreinforced composite, and the following conclusions are obtained.
(1) With the comparison between the measured damping results obtained in this paper and the calculated damping results based on the AdamsBacon model with the same loss factor, the corresponding maximum deviation between them is less than 15%, so the correctness of complex modulus method has been verified indirectly, which can be used to identify the loss factor of fiberreinforced composite with high precision and efficiency.
(2) Complex modulus method mainly contains the following key steps: (I) input the geometric parameters and material parameters; (II) calculate the mass matrix and complex stiffness matrix; (III) solve the natural frequencies and midplane displacement; (IV) obtain total strain energy, dissipation energy, and function of damping ratio and loss factor; (V) identify loss factor by least square method.
(3) In order to ensure that the identified loss factors along the longitudinal, transverse, and shear directions can have better accuracy, the calculated accuracy of natural frequencies of composite plate need to be improved as much as possible. Besides, after comparing with the identified loss factor of fiberreinforced composite obtained by different stepsizes as well as their time expenditure, it has been found out that the stepsize of 0.0001 is the best stepsize.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This study was supported by the National Natural Science Foundation of China Grant no. 51505070, the Fundamental Research Funds for the Central Universities of China Grant no. N150304011 and no. N160313002, and Key Laboratory of Vibration and Control of AeroPropulsion System Ministry of Education, Northeastern University, Grant no. VCAME201603.
References
 R. Chandra, S. P. Singh, and K. Gupta, “Damping studies in fiberreinforced composites—a review,” Composite Structures, vol. 46, no. 1, pp. 41–51, 1999. View at: Publisher Site  Google Scholar
 A. Treviso, B. Van Genechten, D. Mundo, and M. Tournour, “Damping in composite materials: Properties and models,” Composites Part B: Engineering, vol. 78, pp. 144–152, 2015. View at: Publisher Site  Google Scholar
 D. Montalvão, R. A. L. D. Cláudio, A. M. R. Ribeiro, and J. DuarteSilva, “Experimental measurement of the complex Young's modulus on a CFRP laminate considering the constant hysteretic damping model,” Composite Structures, vol. 97, pp. 91–98, 2013. View at: Publisher Site  Google Scholar
 L. UlkeWinter, M. Klaerner, and L. Kroll, “Determining the damping behavior of fiber reinforced composites: a new approach to find mathematical relationships in data sets,” Composite Structures, vol. 100, pp. 34–39, 2013. View at: Publisher Site  Google Scholar
 L. Dunxiang, N. Ronggen, and R. D. Adams, “The finite element technique for predicting the natural frequencies, mode shapes and damping values of filamentary composite plates,” Applied Mathematics and MechanicsEnglish Edition, vol. 7, no. 2, pp. 197–213, 1986. View at: Publisher Site  Google Scholar
 J.M. Berthelot, M. Assarar, Y. Sefrani, and A. E. Mahi, “Damping analysis of composite materials and structures,” Composite Structures, vol. 85, no. 3, pp. 189–204, 2008. View at: Publisher Site  Google Scholar
 R. D. Adams and D. G. C. Bacon, “Effect of Fibre Orientation and Laminate Geometry on the Dynamic Properties of CFRP,” Journal of Composite Materials, vol. 7, no. 4, pp. 402–428, 1973. View at: Publisher Site  Google Scholar
 M. D. Rao and S. He, “Dynamic analysis and design of laminated composite beams with multiple damping layers,” AIAA Journal, vol. 31, no. 4, pp. 736–745, 1993. View at: Publisher Site  Google Scholar
 R. Rikards, A. Chate, A. K. Bledzki, and V. Kushnevsky, “Numerical modelling of damping properties of laminated composites,” Mechanics of Composite Materials, vol. 30, no. 3, pp. 256–266, 1994. View at: Publisher Site  Google Scholar
 J. De Visscher, H. Sol, W. P. De Wilde, and J. Vantomme, “Identification of the Damping Properties of Orthotropic Composite Materials Using a Mixed Numerical Experimental Method,” Applied Composite Materials, vol. 4, no. 1, pp. 13–33, 1997. View at: Publisher Site  Google Scholar
 M. Dalenbring, “Experimental material damping estimation for planar isotropic laminate structures,” International Journal of Solids and Structures, vol. 39, no. 19, pp. 5053–5079, 2002. View at: Publisher Site  Google Scholar
 J. D. D. Melo and D. W. Radford, “Viscoelastic characterization of transversely isotropic composite laminae,” Journal of Composite Materials, vol. 37, no. 2, pp. 129–146, 2003. View at: Publisher Site  Google Scholar
 V. Kostopoulos and D. T. Korontzis, “A new method for the determination of viscoelastic properties of composite laminates: A mixed analyticalexperimental approach,” Composites Science and Technology, vol. 63, no. 10, pp. 1441–1452, 2003. View at: Publisher Site  Google Scholar
 R. Chandra, S. P. Singh, and K. Gupta, “Experimental evaluation of damping of fiberreinforced composites,” Journal of Composites, Technology and Research, vol. 25, no. 2, pp. 96–107, 2003. View at: Google Scholar
 J.M. Berthelot and Y. Sefrani, “Damping analysis of unidirectional glass and Kevlar fibre composites,” Composites Science and Technology, vol. 64, no. 9, pp. 1261–1278, 2004. View at: Publisher Site  Google Scholar
 J.M. Berthelot, “Damping analysis of laminated beams and plates using the Ritz method,” Composite Structures, vol. 74, no. 2, pp. 186–201, 2006. View at: Publisher Site  Google Scholar
 P. Kumar, R. Chandra, and S. P. Singh, “Interphase effect on damping in fiberreinforced composites,” ICCES, vol. 4, pp. 67–72, 2007. View at: Google Scholar
 M. M. Shokrieh and A. Najafi, “Damping characterization and viscoelastic behavior of laminated polymer matrix composites using a modified classical lamination theory,” Experimental Mechanics, vol. 47, no. 6, pp. 831–839, 2007. View at: Publisher Site  Google Scholar
 M. Matter, T. Gmür, J. Cugnoni, and A. Schorderet, “Numericalexperimental identification of the elastic and damping properties in composite plates,” Composite Structures, vol. 90, no. 2, pp. 180–187, 2009. View at: Publisher Site  Google Scholar
 C. Devalve and R. Pitchumani, “Experimental investigation of the damping enhancement in fiberreinforced composites with carbon nanotubes,” Carbon, vol. 63, pp. 71–83, 2013. View at: Publisher Site  Google Scholar
 R. F. Gibson and R. Plunkett, “Dynamic Mechanical Behavior of FiberReinforced Composites: Measurement and Analysis,” Journal of Composite Materials, vol. 10, no. 4, pp. 325–341, 1976. View at: Publisher Site  Google Scholar
 R. F. Gibson and R. Plunkett, “A forcedvibration technique for measurement of material damping  Internal damping of beam specimens in resonant flexural vibration is found from measurements of input acceleration and resulting specimenbending strain,” Experimental Mechanics, vol. 17, no. 8, pp. 297–302, 1977. View at: Publisher Site  Google Scholar
 R. D. Adams and M. R. Maheri, “Damping in advanced polymermatrix composites,” Journal of Alloys and Compounds, vol. 355, no. 12, pp. 126–130, 2003. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2017 Hui Li et al. 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.