Research Article  Open Access
Nonlinear Free Vibration for Viscoelastic Moderately Thick Laminated Composite Plates with Damage Evolution
Abstract
The nonlinear free vibration for viscoelastic crossply moderately thick laminated composite plates under considering transverse shear deformation and damage effect is investigated. Based on the TimoshenkoMindlin theory, strainequivalence hypothesis, and Boltzmann superposition principle, the nonlinear free vibration governing equations for viscoelastic moderately thick laminated plates with damage are established and solved by the Galerkin method, Simpson integration, NewtonCotes, Newmark, and iterative methods. In the numerical results, the effects of transverse shear, material viscoelasticity, spanthickness ratio, aspect ratio, and damage effect on the nonlinear free vibrating frequency of the viscoelastic crossply moderately thick laminated plates are discussed.
1. Introduction
The structure will present the resonance phenomenon when the external excitation frequency is near to a certain natural frequency of the structure during the service life of the structures. Structure destructions caused by the resonance are prevalent in the practical engineering. The damage will emerge in the viscoelastic composite structures during the process of vibration and lead to the change of the dynamic behavior. When the damage develops, the structure will probably enter into the resonant state. As soon as the resonance appears, the stress values in the structure will increase, which will cause the development of the damage accelerate. Therefore, it is a very important research field to investigate the nonlinear dynamic behavior of viscoelastic laminated plates with damage effect.
Extensive studies have been made in dynamics of viscoelastic homogeneous structures. On the basis of the linear theory and the concept of the Lyapunov exponents, Aboudi and Cederbaum [1] investigated the dynamic stability of viscoelastic rectangular plates. Librescu and Chandiramani [2] analyzed the dynamic stability of transversely isotropic viscoelastic plates. Sun and Zhang [3] investigated the chaotic behaviors of viscoelastic rectangular plates subjected to an inplane periodic load and pointed out that the stability of the structure could be increased by adjusting the material parameters. Chen et al. [4] analyzed the steadystate response of the parametrically excited axially moving string constituted by the Boltzmann superposition principle. T. W. Kim and J. H. Kim [5] applied finite element analysis and the method of multiple scales to investigate the nonlinear vibrating frequency of viscoelastic laminated plates. Yu and Huang [6] presented a mathematical model for the vibration of a threelayered sandwich circular plate with viscoelastic core and discussed the effect of viscoelasticity on the frequency and amplitude. Relatively, few works have been devoted to study the effects of local damage and defects on the static and dynamic behavior of plates. Prabhakara and Datta [7, 8] analyzed the effect of the structural flaw on the natural frequency and buckling load of elastic plates subjected to a uniform inplane load. Laura and Gutierrez [9] presented the linear fundamental frequency of transverse vibration for a damaged circular annular plate. However, the materials of the composite laminated plates have the property of viscoelasticity with the apparent creep phenomenon and relaxation characteristic, so it is very necessary to examine the influences of the damage effect on the nonlinear dynamics of viscoelastic laminated plates. Sheng and Cheng [10] used the history curve, phase trajectory diagram, Poincare map, bifurcation figure, and power spectrum to analyze the nonlinear dynamical properties of viscoelastic thick plate with damage. Fu et al. [11, 12] studied the nonlinear dynamic response of viscoelastic composite plate with transverse matrix cracks based on Schapery’s 3D constitutive relationship. To author’s work, Zheng and Fu [13] have studied the effect of local damage on the bifurcation and chaos of viscoelastic isotropic plates, and the nonlinear dynamic properties of viscoelastic isotropic plates and laminated plates with considering damage evolution [14–16].
In the present study, the nonlinear free vibration equations of the viscoelastic crossply moderately thick laminated composite plates with damage effect are established by applying TimoshenkoMindlin theory, strain equivalence hypothesis, and Boltzmann superposition principle. By employing the standard linear solid model to express the viscoelastic material properties, Kachanov’s approach to describe the damage evolution, and using the Galerkin method, Simpson integration, NewtonCotes, Newmark method, and iterative procedure, the solutions of the problem are obtained. Numerical results are presented for different parameters.
2. Basic Equations
Consider a viscoelastic crossply rectangular plate having length in the direction, width in the direction, and thickness in the direction. The middle plane of the undeformed plate contains the axes and the origin of the coordinate system is taken at the upper left corner of the plate. Based on TimoshenkoMindlin kinematic hypotheses taking into account the transverse normal deformation, the displacement components , and that include the effect of transverse shear deformation may be described by the following expressions [17]: where is the time, , and are the values of , and at the middle surface, and and are rotation angles of the normal to the middle surface in the and planes, respectively. The nonlinear straindisplacement relationship can be written as where a comma denotes partial differentiation with respect to the corresponding coordinates and where
By applying the loading equivalent principle and assuming that the internal forces acting on any damaged section are the same as the ones before damage, the relationship between the effective stresses and the Cauchy stresses is given as [18] where the anisotropic damage variables and are similarly defined as in [15, 16].
The above relation is expressed as follows:
Employing the strain energy equivalence principle [19] and Boltzmann superposition principle [20], the stressstrain constitutive equations of the coupled viscoelastic/damage crossply laminated plates for the th layer can be obtained in the following form: in which is the timedependent relaxation function and is the initial Young's modulus of materials for the layer. For orthotropic viscoelastic materials and considering that Poisson ratio contains constant, is the strain transformation relation for the layer, having where and . is the angle measured from the fibre direction to coordinate axis for the th layer; here or
For simplification, (2.6) is rewritten as where the symbol () is the Stieltjes convolution operation symbol, which is defined as The nonzero elements in and are given as
As the classical plate theory, the stress resultants and couples are defined by where, , and are the membrane stress resultants per unit length, and are the transverse shear stress resultants per unit length, , and are the bending and twisting moments per unit length, and is the shear correction factor taking . Substituting (2.8) into (2.11), we can obtain where and the elements in the coupled damaged stiffness tensors , , , and are determined as
Neglecting the effects of inplane inertia, rotary inertia, and coupled normalrotary inertia, the nonlinear equilibrium equations for moderately thick laminated plates are [21] where is the mass of unit volume. By substituting (2.12) into (2.15), and introducing the following dimensionless parameters: then, the dimensionless equilibrium equations of crossply laminated plates with the coupled effects of viscoelasticity and damage can be expressed as
Suppose that all of the boundary conditions of the plate are simply supported. In such case, the dimensionless boundary conditions can be written as
3. Solution Methodology
A solution for (2.17) in conjunction with the boundary condition (2.18) is sought in the following separable form:
Substituting (3.1) into the governing equations (2.17) and making use of the oneterm approximation of the Galerkin method, we can transform the nonlinear integralpartial differential equations into the nonlinear integralordinary differential equations in terms of , and . The domain is divided by square mesh into divisions and the time is equally divided into small time segments . The Simpson integral formula is used to compute the integrations with respect to the spatial coordinates and the NewtonCotes trapezoidal rule is used to compute the integrations with respect to time. Moreover, using the Newmark scheme, the acceleration item and velocity item can be expressed as follows: where denotes the value of at the time . Through these treatings, (2.17) can be transformed into the nonlinear algebraic equations only with time. And the whole equations are iterated to seek solutions. At each step of the iteration, the nonlinear items are linearized. For example, at the step , the nonlinear items may be transformed to where is the average value of those obtained in the preceding two iterations. For the initial step of the iteration, it can be determined by using the quadratic extrapolation, that is, And for the different step of the iteration, the coefficients , and can be expressed as follow: For every time step, the iteration lasts until the difference of the present value and the former is smaller than ; then continue the calculation of the next step.
4. Damage Evolution Equation
In the present research, the following damage evolution equation is employed [19, 22]: Here and are the material constants. is the stress threshold value in the direction at which the damage begins to grow, and
Assume that there is no damage at initial time; thus, the damage values of all points are and . When stopping iterative step , we can obtain the five function values , and . Before continuing the next iterative step, it must determine whether the damage values of an arbitrary point of the plate develop. Therefore, firstly, the displacements of an arbitrary point of the plate can be obtained by using (3.1). Calculate the strains of an arbitrary point of the plate by using (2.2) and (2.3). Then, using (2.5), (2.8) and (4.2), and can be obtained. Finally, use (4.1) to determine whether the damage grows. Suppose that is the damage value of an arbitrary point for the th layer in the direction at the time . If , then the damage value at the time is and if , then
5. Numerical Results
To ensure the accuracy and effectiveness of the present method, the fundamental frequencies of a threelayer symmetric crossply laminated elastic plate without considering damage effect are solved firstly and the materials properties [23, 24] are given as , , , and . Define , in which and are the dimensionless and dimensional linear frequencies of laminated elastic plates without damage, respectively. The fundamental frequencies are calculated and compared in Table 1 with those of [23, 24]. Table 1 shows that the present result approximately agrees with the result in [24] by the firstorder sheardeformation theory (FOT). And for the moderately thick laminated plates, the FOT solution has small difference compared with the secondorder sheardeformation theory (SOT) solution and thirdorder sheardeformation theory (TOT) solution, but large difference is compared with the classical plate theory (CPT). It is noted that the effect of the transverse shear deformation cannot be neglected for the moderately thick laminated plates.
Consider the amplitudes and the frequencies varying with time and the viscoelasticity and damage developing with the increase of time, so we must elect the later vibrating period to be analyzed. Hence, the relationship of the average frequency and the maximum amplitude in the 20th period are presented. Also, assume that every layer has the same material and the identical thickness. The parameters used in numerical calculation are taken as
Define where and are the dimensionless and dimensional nonlinear free vibrating frequency of viscoelastic plates with damage, respectively.
Table 2 shows the effect of transverse shear deformation on the nonlinear free vibrating frequency of the viscoelastic crossply laminated square plate with considering damage effect or without considering damage effect, in which is the dimensionless maximum vibration amplitude of the center point of the plate (similarly below), is the tracing constant which represents the influence of transverse shear when and the effect is neglected when , as well as is the tracing constant which represents the influence of damage when and the effect is neglected when . From Table 2 it may be observed that the thicker the thickness of the laminated plate, the more significant the effect of transverse shear on the nonlinear free vibrating frequency of the laminated plate. And in all cases transverse shear results in a decrease of the vibrating frequency. Therefore, the effect of transverse shear deformation is considered in the following calculated examples.

Figure 1 shows the effect of spanthickness ratio on the nonlinear free vibrating amplitudefrequency response curves of the viscoelastic crossply laminated square plate. In Figure 1, the solid line denotes the undamaged condition, while the dashed line denotes the condition with damage evolution. It can be observed that with the decrease of the spanthickness ratio, namely, the increase of thickness of the plate, the nonlinear free vibrating frequency of the plate becomes bigger under a given vibration amplitude. Meanwhile, the frequency increases with increase of the vibration amplitude under a fixed spanthickness ratio. It can be also seen that the two curves with/without damage agree very well where the vibration amplitude is small. But with the increase of the vibration amplitude, the damage begins to appear and the accumulation velocity of the damage increases, which results in the stiffness of the plate becoming smaller. Accordingly, the nonlinear vibration frequency considering damage reduces remarkably than the result neglecting damage. Moreover, the less the spanthickness ratio of the plate is, the more the damage becomes, then the more the frequency descends.
Figure 2 displays the effect of aspect ratio on the nonlinear free vibrating amplitudefrequency response curves of the viscoelastic crossply laminated plate. With the increase of the aspect ratio, that is, the width decreases under the condition that the length remains constant, the nonlinear free vibrating frequency of the plate becomes bigger. Similarly, under the given aspect ratio , the increase of the amplitude will expedite the accumulation velocity of the damage which will cause the more reduction in the stiffness of the plate. Consequently, the reduction of frequency caused by the damage will be more remarkable. The more the aspect ratio is, the more the frequency reduces.
The effect of material viscoelasticity parameter on the nonlinear free vibrating frequency of the viscoelastic crossply laminated square plate is listed in Table 3. The nonlinear free vibrating frequency of the laminated plate decreases with the higher value of the material viscoelasticity parameter. Due to the fact that viscoelastic material possesses dissipative nature and it acts as damping in the dynamic problems, it can improve the stability of the structure. Therefore, at the same condition, the higher viscoelasticity parameter can suppress the emergence of damage and cause the difference of the vibrating frequency between the damaged plate and undamaged plate decrease.

6. Conclusions
The nonlinear free vibration for viscoelastic crossply laminated composite plates under considering transverse shear deformation and damage effect has been investigated. The effects of transverse shear, material viscoelasticity, spanthickness ratio, aspect ratio, and damage effect on the nonlinear vibration of laminated plates have been discussed. In summary, the transverse shear effect on the nonlinear vibration of the viscoelastic laminated plates with damage is significant, especially, for the laminates with large thickness. This effect decreases the nonlinear frequency but does not change the general behavior of the nonlinear vibration in all cases. With the decrease of spanthickness ratio or increase of aspect ratio, the vibrating frequency increases, as well as the difference of vibration frequency between considering damage and neglecting damage becomes bigger. In addition, with the increase of the vibration amplitude, the damage begins to appear, the accumulation velocity of the damage increases, and the reduction of the vibrating frequency becomes more significant. The larger the material viscoelasticity, the smaller the reduction of the vibrating frequency.
Acknowledgments
The project was supported by the Young Technological Talents’ Innovation Foundation of Fujian Province (2006F3077) and Scientific and Technical Development Foundation of Fuzhou University (2007XQ20).
References
 J. Aboudi and G. Cederbaum, “Dynamic stability analysis of viscoelastic plates by Lyapunov exponents,” Journal of Sound and Vibration, vol. 139, no. 3, pp. 459–467, 1990. View at: Google Scholar
 L. Librescu and N. K. Chandiramani, “Dynamic stability of transversely isotropic viscoelastic plates,” Journal of Sound and Vibration, vol. 130, no. 3, pp. 467–486, 1989. View at: Google Scholar
 Y. X. Sun and S. Y. Zhang, “Chaotic dynamic analysis of viscoelastic plates,” International Journal of Mechanical Sciences, vol. 43, no. 5, pp. 1195–1208, 2001. View at: Publisher Site  Google Scholar
 L. Q. Chen, J. W. Zu, and J. Wu, “Steadystate response of the parametrically excited axially moving string constituted by the Boltzmann superposition principle,” Acta Mechanica, vol. 162, no. 1–4, pp. 143–155, 2003. View at: Publisher Site  Google Scholar
 T. W. Kim and J. H. Kim, “Nonlinear vibration of viscoelastic laminated composite plates,” International Journal of Solids and Structures, vol. 39, no. 10, pp. 2857–2870, 2002. View at: Publisher Site  Google Scholar
 S. C. Yu and S. C. Huang, “Vibration of a threelayered viscoelastic sandwich circular plate,” International Journal of Mechanical Sciences, vol. 43, no. 10, pp. 2215–2236, 2001. View at: Publisher Site  Google Scholar
 D. L. Prabhakara and P. K. Datta, “Vibration and static stability characteristics of rectangular plates with a localized flaw,” Computers and Structures, vol. 49, no. 5, pp. 825–836, 1993. View at: Google Scholar
 D. L. Prabhakara and P. K. Datta, “Parametric instability characteristics of rectangular plates with localized damage subjected to inplane periodic load,” Structural Engineering Review, vol. 1, pp. 71–79, 1993. View at: Google Scholar
 P. A. A. Laura and R. H. Gutierrez, “Transverse vibration of a damaged circular annular plate with a free edge,” Journal of Sound and Vibration, vol. 241, no. 4, pp. 720–724, 2001. View at: Publisher Site  Google Scholar
 D. F. Sheng and C. J. Cheng, “Dynamical behaviors of nonlinear viscoelastic thick plates with damage,” International Journal of Solids and Structures, vol. 41, no. 26, pp. 7287–7308, 2004. View at: Publisher Site  Google Scholar
 Y. Fu, P. Li, and Y. Zheng, “Analysis of nonlinear dynamic response for viscoelastic composite plate with transverse matrix cracks,” Acta Mechanica Solida Sinica, vol. 17, no. 3, pp. 230–238, 2004. View at: Publisher Site  Google Scholar
 Y. Fu and Y. H. Lu, “Analysis of the nonlinear dynamic response of viscoelastic symmetric crossply laminated plates with transverse matrix crack,” Composite Structures, vol. 72, no. 4, pp. 469–476, 2006. View at: Publisher Site  Google Scholar
 Y. F. Zheng and Y. Fu, “Effect of damage on bifurcation and chaos of viscoelastic plates,” International Journal of Nonlinear Sciences and Numerical Simulation, vol. 6, no. 1, pp. 87–92, 2005. View at: Google Scholar
 Y. F. Zheng and Y. Fu, “Effect of damage on nonlinear dynamic properties of viscoelastic rectangular plates,” Applied Mathematics and Mechanics, vol. 26, no. 3, pp. 319–326, 2005. View at: Google Scholar
 Y. F. Zheng and Y. Fu, “Nonlinear dynamic analysis of viscoelastic/damage behaviour for symmetric crossply laminated plates,” Advances in Vibration Engineering, vol. 3, pp. 185–197, 2004. View at: Google Scholar
 Y. F. Zheng and Y. Fu, “Analysis of nonlinear vibration for symmetric angleply laminated viscoelastic plates with damage,” Acta Mechanica Sinica, vol. 21, no. 5, pp. 459–466, 2005. View at: Publisher Site  Google Scholar
 R. D. Mindlin, “Influences of rotatory inertia and shear inflexural motion of isotropic, elastic plates,” Journal of Applied Mechanics, vol. 18, pp. 1031–1036, 1951. View at: Google Scholar
 C. L. Chow and J. Wang, “An anisotropic theory of elasticity for continuum damage mechanics,” International Journal of Fracture, vol. 33, no. 1, pp. 3–16, 1987. View at: Publisher Site  Google Scholar
 L. M. Kachanov, Introduction to Continuum Damage Mechanics, Martinus Nijhoff, Dordrecht, The Netherlands, 1986.
 A. C. Pipkin, Lectures on Viscoelasticity Theory, Springer, New York, NY, USA, 1986.
 Y. Nath and K. K. Shukla, “Nonlinear transient analysis of moderately thick laminated composite plates,” Journal of Sound and Vibration, vol. 247, no. 3, pp. 509–526, 2001. View at: Publisher Site  Google Scholar
 W. H. Zhang, Numerical analysis of continuum damage mechanics, Ph.D. thesis, University of New South Wales, Sydney, Australia, 1992.
 A. A. Khdeir and J. N. Reddy, “Free vibrations of laminated composite plates using secondorder shear deformation theory,” Computers and Structures, vol. 71, no. 6, pp. 617–626, 1999. View at: Publisher Site  Google Scholar
 L. Librescu, A. A. Khdeir, and D. Frederick, “A shear deformable theory of laminated composite shallow shelltype panels and their response analysis I: free vibration and buckling,” Acta Mechanica, vol. 76, no. 12, pp. 1–33, 1989. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2010 Y. F. Zheng and L. Q. Deng. 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.