Research Article  Open Access
Analytical Modeling of HardCoating Cantilever Composite Plate considering the Material Nonlinearity of Hard Coating
Abstract
Due to the material nonlinearity of hard coating, the coated structure produces the nonlinear dynamical behaviors of variable stiffness and damping, which make the modeling of hardcoating composite structure become a challenging task. In this study, the polynomial was adopted to characterize this material nonlinearity and an analytical modeling method was developed for the hardcoating composite plate. Firstly, to relate the hardcoating material parameters obtained by test and the analytical model, the expression of equivalent strain of composite plate was derived. Then, the analytical model of hardcoating composite plate was created by energy method considering the material nonlinearity of hard coating. Next, using the NewtonRaphson method to solve the vibration response and resonant frequencies of composite plate and a specific calculation procedure was also proposed. Finally, a cantilever plate coated with MgO + Al_{2}O_{3} hard coating was chosen as study case; the vibration response and resonant frequencies of composite plate were calculated using the proposed method. The calculation results were compared with the experiment and general linear calculation, and the correctness of the created model was verified. The study shows the proposed method can still maintain an acceptable precision when the material nonlinearity of hard coating is stronger.
1. Introduction
In recent years, vibration control using hard coating which is prepared by the metal, ceramic, or their mixture is receiving a reasonable amount of attention. Compared with traditional surface treatments, such as viscoelastic dampers, the biggest technical advantage is that the hard coating can maintain the consistent damping characteristics in the high temperature or high corrosion environment. Thus, this vibration control technology is expected to apply in the power equipment, such as the turbine engine. From the characterizing test studies [1, 2] about the coated composite structure, some interesting phenomenon has been found; for example, the resonant frequency decreased with the increase of excitation amplitude and the force response of frequency were not symmetric about the resonant frequency. But for the uncoated structure, the abovementioned phenomenon did not happen, which proves that the hard coating makes the composite structure have the nonlinear vibration phenomenon. Further studies show that this nonlinearity comes from the strain dependent characteristic of coating material, that is, the storage modulus (or Young’s modulus) and loss modulus of hard coating change with the strain response amplitude of structure, and belongs to the material nonlinear behavior.
To effectively implement the hardcoating damping technology, the reliable analysis model, which can reproduce the nonlinear dynamic phenomenon observed in the experiment, need be created. However, to the authors’ opinions, there is much bigger gap away from this goal. As far as the current research status is concerned, the modeling about hardcoating composite structure is mainly based on the micro material science. For example, Tassini et al. [3] created a phenomenological model to characterize the nonlinear mesoscopic elastic behaviors of hardcoating material. Torvik [4] developed a slip damping model to describe the idealized microstructure characteristic of hard coatings. AlRub and Palazotto [5] proposed a representative volume elementbased micromechanical model to observe the nonlinear damping in plasma sprayed hard ceramic coatings. All the above studies show that the internal friction in the coating creates the energy dissipation of composite structure, being very important to understand the nonlinear dynamics mechanism or vibration reduction mechanism of hard coating. However, for the hardcoating damping design, it is not realistic that entirely depending on these micro material models and developing the dynamics model based on the macro vibration science are necessary. Furthermore, the challenge of developing this model comes from how to consider the material nonlinearity of hard coating.
To introduce the material nonlinearity of hard coating into the analysis model, the researchers have proposed several methods. Considering the hardcoating dissipate energy by the internal friction of coating particles, then, some researchers attempt to use the Iwan model (shown in Figure 1) to simulate the strain dependent characteristic of hard coating.
In this model, the storage and loss modulus of hard coating are expressed as spring stiffness, , and friction forces, , and both of them are the functions of strain amplitude. is the externally applied force, is the mass of the system, and is the displacement of the mass. Some typical examples show as follows. Green and Patsias [6] applied the Iwan model to describe the response of a coated beam and analyzed the dissipation mechanism of hard coating. Reed [7] did a similar research and used the singlefreedom nonlinear adjusted Iwan model to analyze the free response and forced response of hardcoating beam and verified the rationality of relative damping identification method. To effectively use the Iwan model, some key parameters, such as spring stiffness, friction forces, mass, and applied force, should be determined firstly. But it is not easy and the above studies did not show yet how these parameters were selected. Another method is using the stored energy per unit volume, , and dissipated energy per unit volume, , to describe the material nonlinearity of hard coating, which is given bywhere and are the storage modulus and the loss modulus of hard coating, respectively, and is the equivalent strain. The original intention of this expression is to satisfy the need of deducing the identification formula of the material parameters with strain dependent characteristic. For example, Torvik and Hansel [8, 9] utilized this expression and obtained the identification formula of material properties of a ceramic coating. The last characterizing idea is applying the polynomial to describe the material parameters of hard coating and this characterizing method is also adopted in this study, which will be introduced in the following section. The benefits of using a polynomial to characterize the material parameters of hard coating are mainly reflected in two aspects. On the one hand, the polynomial can more accurately describe the material parameters of hard coating obtained by experiment. For example, all the discrete point values of hardcoating material parameters gotten by Patsias et al. [10], Reed et al. [11], and Tassini et al. [12] can be described by the polynomial. Furthermore, it is found that high order polynomial should be chosen to fit these discrete point values. On the other hand, the material parameters expressed by polynomial can be easily introduced into the analytical model during the macro vibration analysis. Recently, using the finite element method (FEM) to solve the dynamics of hardcoating composite structure introducing the polynomial expression of hardcoating material parameters has been able to achieve. For example, Filippi and Torvik [13] adopted the linear polynomial to simulate the hardcoating parameters with strain dependent characteristic and analyzed the vibration of compressor blade coated hard coating. However, to the authors’ knowledge, the analytical modeling of hardcoating composite structure with the polynomial expression of hardcoating material parameters cannot be found.
On the basis of considering the material nonlinearity of hard coating, the purpose of this study was to develop the analytical modeling method for the hardcoating cantilever thin plate structure. This paper is organized as follows. In Section 2, to relate the material parameters with strain dependent characteristic of hard coating obtained by test and the analytical modeling of hardcoating composite plate, the expression of equivalent strain of composite plate was derived. In Section 3, the analytical model of hardcoating composite plate was created by energy method considering the material nonlinearities of hard coating. In Section 4, NewtonRaphson method was used to solve the vibration response and resonant frequencies of composite plate and a specific calculation procedure was proposed. In Section 5, a cantilever plate coated with MgO + Al_{2}O_{3} hard coating was chosen as study case, and the vibration response and resonant frequencies of composite plate were calculated using the proposed analytical model and calculation procedures. The calculation results were compared with the experiment and linear calculation, and the correctness of the created model was verified. Finally, some important conclusions about this study were listed in Section 6.
2. Derivation of Equivalent Strain for the HardCoating Thin Plate Structure
As above mentioned, during the study of identification of material parameters for the hard coating, the relative curve between the strain amplitude and the storage modulus or loss modulus is determined. Furthermore, these relative curves can be fitted by polynomial, and then the material parameters with strain dependent characteristic are described effectively. Both the finite element modeling and analytical modeling are needed to introduce these mechanical parameters with polynomial expression correctly. For the beam structure, only a unidirectional strain can characterize the dynamics of beam, and corresponding to this unidirectional strain, the material parameters expressed by polynomial can be introduced easily. However, for the plate structure in this study, even the classical plate theory is adopted, and at least three kinds of strain still need to be adopted to describe the dynamics of plate. Then, how to relate the material parameters obtained by test and the analytical model? In this section, the equivalent strain of hardcoating composite plate was chosen to relate the measured material parameters expressed by polynomial and the analytical model. In the following, a brief introduction about the polynomial expression of hardcoating material was given firstly, and then the equivalent strain of hardcoating thin plate structure was derived.
The elastic modulus of hard coating is expressed as complex modulus and the following formula can be used to describe strain dependent characteristic:where refers to complex value.
When the complex modulus of hard coating is expressed by polynomial, the storage and loss modulus in (2) change aswhere and are the storage and loss modulus without considering the strain dependent characteristic of hard coating, and are the specific coefficients of strain dependent storage modulus and loss modulus. Therefore, integrating (2)(3), the complex modulus of hard coating can be also expressed aswhere
In (4), the material parameters of hard coating can be divided two parts: one is , and the other is . For , it does not include the strain dependent characteristic, so it can be named as linear part of hardcoating materials. On the contrary, can be called nonlinear part of hardcoating materials.
The schematic diagram of hardcoating composite plate is shown in Figure 2, and the thin plate is in cantilever status and is excited by base. The reason of choosing this boundary condition is that the base excitation was adopted to test the vibration characteristics in many studies [1, 2]. For the most coating structures, of course including the hardcoating composite plate structure of this study, all of them belong to thin shell structure; that is, the thickness is far less than the length and width or other geometry dimensions. So, herein the hardcoating thin plate can be analyzed according to the classic thin plate theory. Only the material nonlinearity of hard coating should be introduced into the analytical model. Since the classic thin plate theory is suitable to deal with this problem, then the movement of any point in the composite plate can be described by the transverse displacement (or deflection) of neutral surface.
(a) Hardcoating composite plate
(b) Cross section of composite plate
The coordinate plane is set to locate at the neutral surface of composite plate; the length and width of plate are and , respectively. The deflection in the neutral surface of composite plate is . The expression of base excitation is , where is the level of base excitation and is the angular frequency of excitation. Figure 2(b) shows a small part of the plate section. Here, and represent the thickness of the hard coating and the metal substrate, respectively, and is the distance between the interface of coatingsubstrate and the neutral surface.
According to the classic thin plate theory, the strain of any point in the composite plate can be expressed aswhere , are the normal strain of and direction, respectively, is the shear strain in the plane, and is the distance between the point of the cross section and neutral surface.
Furthermore, according to the physical equations, the stress of any point in the composite plate can be given:where , are the normal stress of and direction, respectively, is the shear strain in the plane, , are equivalent Young’s modulus and Poisson’s ratio, respectively (it should be noted that Poisson’s ratio in this work is thought as independent with the strain amplitude), and is the equivalent shear modulus and the relation with equivalent Young’s modulus and Poisson’s ratio can be expressed asHere, the equivalent strain of hardcoating composite plate was determined according to the principle of equal strain energy density; the expression iswhere is the equivalent stress and satisfies
Equations (6), (7a), (7b), and (7c), and (10) can be substituted into (9) and then yieldThus, the equation of solution of equivalent strain can be obtained and shown as
It should be noted that because of introducing the complex response in the subsequent analytical modeling of composite plate, the deflection of composite plate is also a complex. For simplification of expression, only material parameters were marked as complex in this study. The equivalent strain amplitude should be a real, so the solution formula is changed aswhere the operator refers to solving the modulus of complex values.
The calculation in (13) is very complex, which will reduce the computational efficiency obviously for the subsequent iterations, so a relatively simple expression of equivalent strain amplitude is further given:The final expression is corresponding to the maximum value of equivalent strain amplitude in (13) and practice has shown that this simplification does not introduce big error but improves the computational efficiency greatly.
3. Analytical Modeling of HardCoating Composite Plate
In the following, the energy method was adopted to create the analytical model of the hardcoating composite plate (shown in Figure 2) under the base excitation. Usually, the loss factor of hard coating is higher than the common metals, but it is much less than the viscoelastic materials, so the damping of hard coating belongs to medium or low levels. Then, the damping contribution of both hard coating and metal substrate was considered in this study. Here, the complex modulus model was used to describe both hard coating and metal substrate. The complex modulus expression of hard coating can be seen from (2) to (5), and the modulus of metal substrate is expressed aswhere , are the storage modulus and the loss modulus of metal substrate, and both of them were set as constant in this study.
The distance between the interface of coatingsubstrate and the neutral surface can be determined as [14]
Due to undertaking the base excitation, the displacement of any point in the composite plate can be expressed asIt can be noted from (17) that if the deflection of any point is known, then the displacement can be determined.
The energy in the hardcoating thin plate system should include strain energy, kinetic energy, and damping dissipation energy, and these energies will be solved below. According to the classic thin plate theory, the strain energy of hardcoating plate system can be given aswhere refers to the volume of composite plate and and are named as the linear and nonlinear strain energy, respectively, which are corresponding to the linear and nonlinear part of hardcoating material parameters.
Introducing the constitutive model of substrate and linear part of hard coating, respectively, and integrating along axis, the solution formula of linear strain energy can be yielded: where refers to the area of composite plate, and the variables, such as , and , can be expressed aswhere , are Poisson’s ratios of metal substrate and hard coating. Similarly, the strain dependent characteristics of these Poisson ratios are not considered in the analytical model. , are shear modulus of metal substrate and hard coating corresponding to the linear part of hardcoating material.
Furthermore, the nonlinear part of hardcoating material is introduced into the solution of nonlinear strain energy of hardcoating plate. To correspond to the polynomial expression of hardcoating material, the equivalent strain is adopted to calculate the nonlinear strain energy, and the solution formula can be described aswhere is the volume of hard coating and the expressions of other parameters are shown as follows:where and refer to the function of extracting real and imaginary part. It is worth noting that in (21) the maximum order of the polynomial expression of hardcoating materials is 3. The practice has shown that for most conditions the 3order polynomial can describe the material parameters with strain dependent characteristic accurately. In addition, it can be known from (18) to (22a), (22b), (22c), (22d), and (22e) that the material damping of both hard coating and metal substrate has included the above strain energy expression.
The kinemics energy of hardcoating composite plate can be described aswhere and are density of substrate and hardcoating material, respectively. Furthermore, one can define , which presents the mass of the unit area.
For the hardcoating composite plate structure, besides material damping of hard coating and substrate, the boundary conditions damping and fluid damping in air also should be contained and this damping can be equivalent to viscous damping, named as the remaining equivalent viscous damping in this work. The remaining equivalent viscous damping can also dissipate the energy of composite system and the expression is shown as follows:where is the coefficient of remaining equivalent viscous damping of composite plate.
Because only steadystate response of composite plate was considered in this study, Galerkin discrete method was adopted to reduce the order of the above nonlinear problem. In fact, in [13], similar method is also adopted to deal with the nonlinearity of hardcoating composite structure. It is assumed that the deflection expression of composite plate is where is the order modal shape, which can be assumed consistent with the modal shape of uncoated plate, is the order modal participation factor, and is the number of considered orders.
Substituting (25) into (19), the expression of linear strain energy can be changed aswhere the meanings of operators of , , and are , , and , respectively.
Moreover, to obtain the newer expression of the nonlinear strain energy, (25) is substituted into (21). It can be known from (22a), (22b), (22c), (22d), and (22e) that only (22d) and (22e) need be transformed. Because is real, the differential about it should be also real, so the imaginary part only exists in the , and then the expression of and can be given as
Similarly, substituting (25) into (23) and (24), the newer expressions of kinemics energy and dissipation energy by the remaining equivalent viscous damping can be listed as follows:
Arriving here, the expressions, such as strain energy , kinemics energy , and dissipation energy by the remaining equivalent viscous damping, are obtained by Galerkin discretion for the hardcoating plate structure. Substitute the abovementioned expression into the following Lagrange equation:and the final movement equation of hardcoating composite plate will be gotten.
Next, the derivatives of various energies to time , and , are solved. It should be noted that the equivalent strain is independent of time and the variable in (26) does not include the time item yet, so they can be thought as constants during the differential calculation. Furthermore, if only the fundamental frequency vibration is considered, it can be assumed that , which is introduced in the differential calculation. Then, by rearranging, the motion equation of the following style can be yielded:
It can be noted that using this simple expression of solution the involved nonlinear problem has been transformed into solving the response vector independent of time. The other item in the above equation, , is a nonlinear stiffness matrix considering the material nonlinearity of hard coating, refers to the material damping matrix, refers to the remaining equivalent viscous damping matrix, is the mass matrix, and is the excitation vector.
The matrixes , , , and are all order square matrix and the vectors and are order vector. Additionally, can be considered as one item and expressed as complex matrix , and it is also an order square matrix.
The element of the row and column of complex matrix can be shown as
The element of the row and column of mass matrix is
The element of the row and column of remaining equivalent viscous damping matrix is
The damping coefficient can be simplified; for the order, the damping coefficient can be expressed bywhere is the order natural frequency and is the order modal damping ratio which can be obtained by test.
For convenience, if ignoring the nondiagonal elements, the equivalent viscous damping matrix can be further simplified. Then the matrix became an order diagonal matrix and the element of the row and column isHere, it was assumed that the boundary damping and air damping do not change before and after coating; then the modal damping ratio can be introduced according to the damping of uncoated plate; similarly, refers to the natural frequency of uncoated plate.
The element of the excitation vector can be described as
In summary, the dynamics model of composite plate considering the strain material nonlinearity of hard coating was created. In this model, both the material damping and the remaining equivalent viscous damping are considered. Based on this model, the resonant frequencies of composite plate with the material nonlinearity and the nonlinear vibration response can be solved.
4. Solution of Vibration Response and Resonant Frequency of Composite Plate
4.1. Solution of Vibration Response
Due to the material nonlinearity of hard coating, (30) became a nonlinear algebraic equation. Here, the NewtonRaphson method was adopted to solve the vibration response of hardcoating composite plate. Equation (30) can be transformed asBecause the residual value vector contains the and items, also, each element of vector can be expressed as , so the Jacobi matrix of vector is needed and easily separated into the real part and imaginary part, which is shown aswhere is an order square matrix, and the element of row and column isHere,
Similarly, the element of row and column of matrix isHere,
Extract the real and imaginary part of the vectors and , respectively, and form two new vectors, which are shown as follows:
Then, the iterative formulas of NewtonRaphson method for (37) can be gotten:where is the iteration number.
Using the iterative calculation of (44a) and (44b), the solution with certain precision can be obtained and the precision of solution can be controlled by setting the maximum norm value of vector . When attaining termination conditions, the vector should be transformed into the complex solution , and then the deflection and strain of plate can be solved by (25) and (6). Finally, the vibration response of composite plate under base excitation is obtained by (17). The whole solution procedure is shown in Figure 3 and some steps are explained in the same figure.
For step 1, it is necessary to confirm the iterative initial value , which can be obtained according to the following equation:where the stiffness matrix does not contain the strain dependent characteristic of hardcoating material; that is, the linear part of mechanical parameter of hard coating is introduced in the model according to . The element of the row and column of stiffness matrix is
For step 8, the convergence condition is defined as norm 2 of vector less than the set value, that is, .
In the last step, the output values such as the deflection, strain, and displacement response of hardcoating composite plate are obtained by calculating the modules of the relative complex solution.
4.2. Solution of Resonant Frequency
When calculating resonant frequency of hardcoating composite plate, the resonant response of composite should be obtained in advance, so the calculation method of vibration response was described firstly at above section.
Neglecting the damping and excitation force items in (30), the characteristic equation considering the material nonlinearity of hard coating can be expressed as Here, the iterative method is also adopted to solve the resonant frequency of composite plate. The calculating procedure is shown in Figure 4.
For step 1, to obtain the initial value of natural frequency, the material nonlinearity of hard coating is also not considered; the solution expression is
For step 2, corresponding to the initial value , the calculation equation of resonant response is
Similar to the response calculation, in step 6, the convergence condition is defined as norm 2 of the difference between two adjacent frequencies; the expression is .
5. Study Case and Discuss
In the following, the titanium plate coated with MgO + Al_{2}O_{3} hard coating on one side was considered as a study case. The developed method was demonstrated to solve the dynamics characteristics considering the material nonlinearity of hard coating.
5.1. Characterizing the Material Parameters of MgO + Al_{2}O_{3} Hard Coating by Polynomial
The material parameters of MgO + Al_{2}O_{3} hard coating are not constants and have strain dependent characteristic. In this study, the identification results of hardcoating material in [11] were introduced in the developed model. The storage modulus and loss modulus of MgO + Al_{2}O_{3} are shown in Figure 5.
(a) Storage modulus
(b) Loss modulus
It can be seen that these measured values cannot be fitted by linear function and have to be fitted by high order polynomial. Here, the cubic polynomial was chosen and the fitting curve is also shown in Figure 5.
The final expressions of material parameters of MgO + Al_{2}O_{3} hard coating are
5.2. Vibration Experiment of HardCoating Composite Plate
The Ti6Al4V cantilever plate was chosen as study object and one side of plate was coated with MgO + Al_{2}O_{3} hard coating. The method of preparing hard coating was Air Plasma Spraying (APS) and the proportion of coating materials and depositing technology were similar to the ones mentioned in [1, 11]. The geometry parameters of substrate and hard coating are listed in Table 1 and the densities of substrate and hard coating are 4420 kg/m^{3} and 2565 kg/m^{3}, respectively. The material parameters of hard coating have the strain dependent characteristic, which are considered according to (50).

The experimental system is shown in Figure 6. In the system, the electromagnetic vibration shaker is base excitation equipment and the Doppler laser vibrometer is vibration picking equipment. The vibration picking point locates at the lower part of composite plate, and the coordinate values are mm, mm, according to the coordinates shown in Figure 6. Table 2 lists the instruments used in this test in detail. It should be noted that an acceleration sensor was placed in the fixture and this sensor was used to monitor the base excitation amplitude of fixture. Another acceleration sensor was used to control the amplitude of vibration shaker. It has been found that the measurement values are different between the two accelerators. In this study, the measurement value in the fixture was chosen as the base excitation amplitude. To ensure having consistent boundary conditions for the uncoated and coated plate, the torque wrench was chosen and the tightening torque of bolt is 34 N·m in the fixture.

Herein, the 6order resonant frequency and resonant response calculation of hardcoating composite plate were chosen to display the feasibility of the proposed method. The reason of choosing the 6order vibration characteristic was that the effect of boundary condition on the dynamics characteristics of higher order of composite plate is small. According to the needs of analytical calculation and verifying the analytical calculation results, the first 8order modal damping ratios of uncoated plate, the natural frequencies of composite plate, and the 6order resonant frequency and resonant response under different excitation levels (0.1 g, 0.3 g, 0.4 g, 0.5 g, 0.6 g, and 0.7 g) are tested. Among them, the modal damping ratios of uncoated plate were obtained by halfpower bandwidth method from the sweeping response of composite plate under 0.1 g. The natural frequencies of composite plate were gotten from the 3D waterfall figure of sweeping response under 0.1 g. In this study, 0.1 g excitation level was thought to be very low, which can not lead to the appearance of material nonlinearity of hard coating, so the natural frequencies obtained by experiment were thought as linear natural frequencies. The 6order resonant frequency and resonant response were also acquired by the 3D waterfall figure of sweeping response and the sweeping direction was from high frequency to low frequency. All the tested results are listed in Table 3 to Table 5.

It can be seen from Table 5 that the 6order resonant frequency of composite plate decreases with the increase of excitation levels. The sweeping response in the considered frequency range is shown in Figure 7. It should be explained that the labels of excitation levels in Figure 7 were the setting amplitude of vibration table, which are used to control the vibration shaker. Also, it can be seen from Table 5 that the excitation level in the fixture was less than the excitation level in the table. Maybe, the reason is that the energy is dissipated when transferring from table to fixture.
5.3. Analytical Calculation
The strain dependent characteristic of hard coating is considered for the proposed method, so the solving process using the proposed method is named as nonlinear calculation. In contrast, if not considering the strain dependent characteristic, the solving process is named as linear calculation. To compare with the proposed method, the linear calculation was also done for the hardcoating composite plate. For the linear calculation, the material parameters of hard coating were taken as follows: storage modulus 50.439 GPa and loss modulus 0.918 GPa, which are corresponding to the 0 strains in Figure 5. The picking vibration point was consistent with experiment. Using these parameters, the resonant frequencies and the 6order resonant response under different excitation levels were solved. For the nonlinear calculation, the convergence conditions of both vibration response and resonant frequency calculation were defined as . The relative results are listed or shown in Tables 4 and 6 and Figure 8, and so forth. It can be noted from Table 4 that the maximum difference of the first 8order natural frequencies between the experiment and linear calculation is less than 1%, and the discussion about the other results can be seen in the following section.



(a) Linear calculation
(b) Nonlinear calculation
5.4. Results Comparison and Difference Analysis
In this section, the linear and nonlinear calculation results for the 6order resonant frequency and resonant response were compared with the experimental values. The comparison of resonant response is listed in Table 7 and the comparison of resonant frequency is listed in Table 8.


From Table 7, a general rule can be still drawn; that is, the difference between linear calculation of resonant response and experiment values becomes bigger with the increase of excitation levels, although there are exception conditions when the excitation levels are 0.3 g and 0.6 g. For example, when the excitation level is less than 0.4 g, the maximum difference is 7.49%, and when the excitation level is bigger than this value, the difference increases obviously and the maximum value has achieved 19.59% during this calculation. Compared with the linear calculation, the nonlinear calculation proposed in this study can always maintain an acceptable precision and the maximum difference with experiment value is less than 7.5%. For the resonant frequency calculation listed in Table 8, similar conclusion can be drawn; that is, when the excitation level is bigger, the nonlinear calculation should be adopted, although for this case the precisions of both the linear and nonlinear calculation are acceptable. By the above comparison, the rationality of the calculation method considering the material nonlinearity of hard coating was verified.
6. Conclusions
On the basis of considering the material nonlinearity of hard coating, the dynamics modeling and analysis for a hardcoating composite plate were studied, and some important conclusions are listed as follows.(1)Using the high order polynomial to characterize the material parameters of hard coating is an effective method. On the one hand, the high order polynomial can describe the change of mechanical parameter with the strain amplitude; on the other hand, the mechanical parameters expressed by polynomial can be easily introduced into the analytical model.(2)The created analytical model and solving method in this study can correctly solve the dynamics characteristics. The practice for the 6order calculation of hardcoating composite plate shows that the maximum difference of the resonant response calculation is less than 7.5% and the maximum difference of the resonant frequency calculation is less than 0.93%.(3)The resonant response values between linear and nonlinear calculation are significantly different and the difference between linear calculation and experiment values becomes bigger with the increase of excitation levels. So, when the excitation level is bigger, the material nonlinearity of hard coating must be considered; using the nonlinear calculation is needed to acquire the dynamics characteristics of hardcoating composite structure.(4)By comparison of the resonant responses obtained by linear and nonlinear calculation for the hardcoating composite plate, it can be found that the resonant response obtained by nonlinear calculation is less than the result of linear calculation, which maybe means that the material nonlinearity of hard coating is good to the vibration reduction effect.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This project was supported by National Natural Science Foundation of China (Grant no. 51375079) and the Fundamental Research Funds for the Central Universities of China (Grant no. N140301001). In addition, the first author (Wei Sun) would like to thank the University of Akron for providing a Visiting Professor position to carry out this research work.
References
 C. Blackwell, A. Palazotto, T. J. George, and C. J. Cross, “The evaluation of the damping characteristics of a hard coating on titanium,” Shock and Vibration, vol. 14, no. 1, pp. 37–51, 2007. View at: Publisher Site  Google Scholar
 F. Ivancic and A. Palazotto, “Experimental considerations for determining the damping coefficients of hard coatings,” Journal of Aerospace Engineering, vol. 18, no. 1, pp. 8–17, 2005. View at: Publisher Site  Google Scholar
 N. Tassini, S. Patsias, and K. Lambrinou, “Ceramic coatings: a phenomenological modeling for damping behavior related to microstructural features,” Materials Science and Engineering A, vol. 442, no. 12, pp. 509–513, 2006. View at: Publisher Site  Google Scholar
 P. J. Torvik, “A slip damping model for plasma sprayed ceramics,” Journal of Applied Mechanics, vol. 76, no. 6, Article ID 061018, 2009. View at: Publisher Site  Google Scholar
 R. K. A. AlRub and A. N. Palazotto, “Micromechanical theoretical and computational modeling of energy dissipation due to nonlinear vibration of hard ceramic coatings with microstructural recursive faults,” International Journal of Solids and Structures, vol. 47, no. 16, pp. 2131–2142, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Green and S. Patsias, “A preliminary approach for the modeling of a hard damping coating using friction elements,” in Proceedings of the 7th National Turbine Engine High Cycle Fatigue Conference, pp. 7c1–7c9, Palm Beach, Fla, USA, May 2002. View at: Google Scholar
 S. A. Reed, Development of experimental, analytical, and numerical approximations appropriate for nonlinear damping coatings [Ph.D. thesis], Air Force Institute of Technology, Greene County, Ohio, USA, 2007.
 P. J. Torvik and J. Hansel, “Mechanical properties of a ceramic coating with vem infiltration,” Journal of Engineering Materials and Technology, vol. 131, no. 3, Article ID 031003, 2009. View at: Publisher Site  Google Scholar
 P. J. Torvik, “Determination of mechanical properties of nonlinear coatings from measurements with coated beams,” International Journal of Solids and Structures, vol. 46, no. 5, pp. 1066–1077, 2009. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Patsias, C. Saxton, and M. Shipton, “Hard damping coatings: an experimental procedure for extraction of damping characteristics and modulus of elasticity,” Materials Science and Engineering A, vol. 370, no. 12, pp. 412–416, 2004. View at: Publisher Site  Google Scholar
 S. A. Reed, A. N. Palazotto, and W. P. Baker, “An experimental technique for the evaluation of strain dependent material properties of hard coatings,” Shock and Vibration, vol. 15, no. 6, pp. 697–712, 2008. View at: Publisher Site  Google Scholar
 N. Tassini, K. Lambrinou, I. Mircea, M. Bartsch, S. Patsias, and O. Van der Biest, “Study of the amplitudedependent mechanical behaviour of yttriastabilised zirconia thermal barrier coatings,” Journal of the European Ceramic Society, vol. 27, no. 23, pp. 1487–1491, 2007. View at: Publisher Site  Google Scholar
 S. Filippi and P. J. Torvik, “A methodology for predicting the response of blades with nonlinear coatings,” Journal of Engineering for Gas Turbines and Power, vol. 133, no. 4, Article ID 042503, 2011. View at: Publisher Site  Google Scholar
 G. Parthasarathy, C. V. R. Reddy, and N. Ganesan, “Partial coverage of rectangular plates by unconstrained layer damping treatments,” Journal of Sound and Vibration, vol. 102, no. 2, pp. 203–216, 1985. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2015 Wei Sun 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.