International Journal of Aerospace Engineering
Volume 2015 (2015), Article ID 246352, 14 pages
http://dx.doi.org/10.1155/2015/246352

Research Article

Nonlinear Vibrations of FGM Cylindrical Panel with Simply Supported Edges in Air Flow

Y. X. Hao,¹ W. Zhang,² S. B. Li,³ and J. H. Zhang¹

¹College of Mechanical Engineering, Beijing Information Science and Technology University, Beijing 100192, China
²College of Mechanical Engineering, Beijing University of Technology, Beijing 100124, China
³College of Science, Civil Aviation University of China, Tianjin 300300, China

Received 6 June 2014; Accepted 9 December 2014

Academic Editor: Mahmut Reyhanoglu

Copyright © 2015 Y. X. Hao 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.

Abstract.

Chaotic and periodic motions of an FGM cylindrical panel in hypersonic flow are investigated. The cylindrical panel is also subjected to in-plane external loads and a linear temperature variation in the thickness direction. The temperature dependent material properties of panel which are assumed to be changed through the thickness direction only can be determined by a simple power distribution in terms of the volume fractions. With Hamilton’s principle for an elastic body, a nonlinear dynamical model based on Reddy’s first-order shear deformation shell theory and von Karman type geometric nonlinear relationship is derived in the form of partial equations. A third-order piston theory is adopted to evaluate the hypersonic aerodynamic load. Here, Galerkin’s method is employed to discretize this continuous nonlinear dynamic system to ordinary differential governing equations involving two degrees of freedom. The chaotic and periodic response are studied by the direct numerical simulation method for influences of different Mach number and the value of in-plane load. The bifurcations, Poincare section, waveform, and phase plots are presented.

1. Introduction

With the continuous variation of the material properties along the thickness, functionally graded materials (FGM) can be used in high temperature gradient environments especially when they are made of metal and ceramic. The metal can keep a certain extent of toughness and ceramics have superior heat resistant ability. So they usually act as thermal protection structures in spacecraft and other structural components in high temperature environments [1, 2]. It is well known that, due to the combined load of airflow and heating, the flexible panels might exhibit large aerothermal deflections [3]. In an extensive search of panel flutter literature, a number of investigations were dedicated to FGM plates with supersonic or hypersonic flow regimes. Considered a curved skin panel with geometrical imperfection, Abbas et al. [4] gave its flutter in unsteady flow by numerical simulation and Galerkin method. With the help of structural nonlinear and the third-order piston theory, the governing equations were derived. The effects of the system parameters on the flutter were discussed in detail. Using the linear approach, the flutter of rectangular flat plates in supersonic flow and thermal environment was studied by Prakash and Ganapathi [5]. The plate was subjected to the two-dimensional aerodynamic force. They showed that under real flight conditions heating caused by aerodynamics is enormous. Sohn and Kim [6] took a static and dynamic stability study on the panel under aerodynamic force as well as thermal loads. Ibrahim et al. [7, 8] investigated the thermal buckling and nonlinear flutter of thin FGM panels under the action of aerothermoelasticity by the finite element method. Then they presented the results for different factors. Accounting for both the geometric and aerodynamic nonlinearities, Prakash et al. [9] studied the nonlinear flutter of FGM plates under high supersonic airflow in frequency domain and time domain, respectively. The influence of various parameters including the geometrics and physics on the flutter of FGM plates was discussed. Hosseini and Fazelzadeh [10] used the numerical and analytical methodologies to study postcritical and vibration behaviors of the FGM panels in a supersonic air flow. Ibrahim et al. [8] analyzed nonlinear flutter and thermal buckling of an FGM panel under the combined effect of elevated temperature conditions and aerodynamic loading. Lee and Kim [11] dealt with the aerothermopostbuckling behaviors of the FGM panel in supersonic air flow. Limit-cycle vibration was found in this study by Newmark time integration method and Guyan reduction technique. Mey et al. [12] found that, in many studies which are about the flutter characteristics of FGM plates, the linear or nonlinear theory in conjunction with the first-order piston theory used to approximate the aerodynamic pressure was considered. The first-order piston theory is valid for sufficiently high supersonic Mach numbers on surfaces with small geometric characteristics.

Since the first report of flutter instability for circular cylindrical shells, the studies of the aeroelastic stability of cylindrical shells in axial flow received extensive attention [13]. Numerous studies on the cylindrical shells focused on their flutter. Marzocca et al. [14] reviewed the corresponding researches that analyse the dynamic behavior of curved and flat panels exposed to supersonic flow fields. Hosseini et al. [15] analyzed the nonlinear response of FGM curved panels in high temperature supersonic air flows. The effects of curved panel height-rise and volume fraction index on the nonlinear dynamical behavior of the panel which is subjected to aerothermoelastic loads are investigated.

Librescu et al. [16] presented a theoretical investigation of the flutter and postflutter of the long thin-walled circular cylindrical panels in supersonic/hypersonic flow field. For the problems of the flutter boundaries about the simply supported functionally graded truncated conical shell subjected to supersonic air flow, Mahmoudkhani et al. [17] made aerothermoelastic analysis to predict the flutter boundaries. The flutter boundaries were obtained for the FGM conical shells with different semivertex cone angles, different temperature distributions, and different volume fraction indices.

When it comes to FGM cylindrical panel identification, damage detection, and the control of the dynamics, it is necessary to investigate their complex nonlinear flutter in hypersonic air flow in great detail [18]. However, to the best of the authors’ knowledge, works on dynamic instability of FGM cylindrical panel subjected to supersonic/hypersonic flow, including the effects of thermal load and in-plane loads, appear to be scarce in the open literature. Many interesting researches used the linear shell theory. There are a few literatures on the nonlinear dynamic behavior of FGM cylindrical panel taking into account the aerodynamic nonlinearities.

In the present research, the bifurcations and chaotic dynamics of the hypersonic FGM cylindrical panel subjected to thermal and mechanical loads are investigated by applying geometrical nonlinear and the third-order piston theory. Materials properties of the constituents are graded in the thickness direction according to a power law distribution. Only transverse nonlinear oscillations of the FGM cylindrical panel are considered; the equations of motion can be reduced into a two-degree-of-freedom nonlinear system. By the numerical method, the nonlinear dynamical equations are analyzed to find the nonlinear responses of the system.

2. Theoretical Formulation

2.1. Model of the FGM Cylindrical Panel

Consider a simply supported hypersonic FGM circular cylindrical panel of a length , thickness , midsurface radius , and angular width . This panel is also subjected to the in-plane harmonic excitations. Cartesian coordinate is adopted to describe the deformations of the FGM circular cylindrical panel. The coordinates and are in longitudinal and tangential directions, respectively. The coordinate is taken to be positive outward radially, as shown in Figure 1. The displacements of an arbitrary point are denoted by , and . Assume that , and are the displacements of a point in the middle plane in the axial, circumferential, and radial directions, respectively. The rotations of the transverse normal to the midplane about and axes are assumed to be and , respectively. The in-plane excitation of the FGM panel distributed uniformly along the direction at the end of panels and is of the form , where is the static in-plane preload and is part of time dependent that has the form . Here, is the frequency of excitation along axial direction.

Figure 1: The model of an FGM cylindrical panel with simply supported edges and the coordinate system.

Assume that the panel material is made of a composite of the ceramics and metals. The material properties of constituents of the panel including density , elastic moduli , and thermal coefficient of expansion are temperature dependent and can be expressed as [19]

where , and are temperature dependent coefficients and is the environment temperature.

The effective material properties vary continuously in the direction by following a simple power law in terms of the volume fractions and can be expressed as

where and indicate, respectively, the properties of the ceramic and metal and and are their volume fractions and have the following relationship:

the metal volume fraction can be written as

where is the volume fraction exponent which can depict the material variation profile through the thickness. According to Zhang et al. [20], linear temperature change is considered which has the form

where and indicate the temperature of the top and bottom surfaces of the panel, respectively. The FGM circular cylindrical panel is subjected to a uniform temperature variation , where is the reference temperature.

2.2. Aerodynamics Loading

The linear piston theory is valid for Mach numbers changing from to 5 and for higher Mach numbers the nonlinear piston theory must be used in Mei et al. [21]. To study the nonlinear oscillation of functionally graded material cylindrical panel under a hypersonic air flow, one should use third-order piston theory. Piston theory aerodynamics, which is used in problems of oscillating airfoils advanced by Lighthill [22] and later used as an aeroelastic tool by Ashley and Zartarian [23], is a popular modeling technique for supersonic and hypersonic aeroelastic analyses.

For a cylindrical panel, which is exposed to an external hypersonic flow field parallel to the centerline of the panel on the surface, the aerodynamic pressure using the third-order piston theory is expressed by Amabili [13] and Cao and Zhao [24] as

where is the adiabatic exponent and and are the Mach number and time, respectively. The free-stream static pressure is given as

where , and are the free-stream air density, velocity, and the free-stream speed of sound, respectively.

2.3. Geometry and Constitutive Relations

The displacement components for the FGM cylindrical panel based on Reddy’s first-order shear deformation theory [25] can be represented as

It is assumed that the transverse normal stress is negligible and normals are not vertical to the midplane after deformation. Substituting displacement components into Von Karman nonlinear strains-displacement relations, the strains in terms of middle-surface displacements are given as

The constitutive relations of the panel in which the thermal effects due to temperature difference are considered can be written as

where are the elastic constants which can be expressed as

The thermal expansion coefficient can be given as

2.4. Equations of Motion

By using Hamilton’s principle, the motion equations in terms of midplane displacements are obtained as follows [20]:

where is the damping coefficient, all kinds of the stiffness elements , and of the FGM cylindrical panel are denoted by

and various mass inertia terms in (12) can be defined as

Thermal force resultants due to temperature rise are functions of the temperature and coefficient of thermal expansion equation. They can be calculated by

Here, the shear correction factor that is introduced by Reddy [25] and Kadoli and Ganesan [26] is equal to .

For simply supported hypersonic FGM circular cylindrical panel with rectangular base, the first two mode shapes that satisfy the boundary conditions are assumed to be

where and are the time dependence amplitudes of the first two modes. The constant can be represented by .

Compared to the transverse inertia term, the influences of the in-plane and rotary inertia terms on the vibration of the panel are small and can be neglected; see [27, 28]. Following the work given in [20], the displacement components of can be expressed in terms of . Applying the Galerkin method on (12), one can obtain a coupled set of nonlinear ordinary differential equations in time that can take the form as follows:

where and . The coefficients and are the results of thermal stress resultants and the static components of in-plane preloads. The variables and are related to the aerodynamic pressure; they can be expressed as

All the coefficients in (18a) and (18b) can be affected by geometric and physical parameters and temperature field of the panel. And the coefficients in (19a) and (19b) are dependent on geometric and free-stream air density and free-stream speed of sound. They are too long to be listed out in the paper for abbreviation. This is a nonlinear dynamic system which includes cubic and quadratic terms. To study the nonlinear aeroelastic behavior of the FGM cylindrical panel in hypersonic flow, the obtained nonlinear ordinary differential equations can be solved by the fourth-order Runge-Kutta method.

3. Numerical Results and Discussions

In order to validate the numerical results presented in this study, a comparison is shown in Figure 2 for a simply supported intact aluminum-zirconia FGM square plate with length and thickness . A suddenly uniform transverse load of intensity of is applied on it. The material properties vary in the thickness by following a simple power law and the power law exponent is taken as . It means that the plate is zirconia-rich on the top surface and aluminum-rich at the bottom surface. In Figure 2, temporal evolution curve of center deflection of the square plate is given by Reddy’s finite element results (see [29]) and the results of this study. The dimensionless time and center deflection are defined as and , respectively. According to this figure, close agreements between present methodology and Reddy’s finite element results can be observed. Although the plate is subjected to the aerodynamic force in this study and not uniform transverse load, as external pressures, both of them play the same role on the plate.

Figure 2: Comparison of nondimensional center deflection of the square plate with simply supported edges under a suddenly applied uniform load by Reddy’s FEM results [29] and present (—: present; ■: Reddy’s).

In this study, Al₂O₃ and Ti-6Al-4V are chosen as the two constituent materials of the FGM cylindrical panel. The properties for these two constituent materials can be found in Shen [30]. For simplicity, it is assumed that Poisson’s ratio is a constant as the value of and the reference temperature remained a constant as . The simply supported FGM cylindrical panel has the following geometries: length , radius of curvature , thickness , and angular width . In all of the numerical analyses, the temperature difference through the thickness of the panel was assumed to have attained a steady state . In all of the calculations, the initial conditions are fixed at the value of , , and . In addition, the flow field characteristics are as follows: the free-stream air density , the free-stream speed of sound , and the adiabatic exponent ; see Librescu et al. [31]. In given conditions, the influences of the in-plane excitation and Mach numbers on the nonlinear dynamic responses of the FGM cylindrical panel in hypersonic flow are investigated. In order to perform analysis of the chaotic and periodic responses for the FGM cylindrical panel, the bifurcation diagrams are depicted. The dimensionless transverse amplitudes are defined as here. To study convenience, the overbar is dropped in the following research.

Firstly, the Mach number and frequency of the in-plane excitation are taken to be 5.0 and , respectively. And the volume fraction index is as . Figure 3 plots the bifurcation diagrams of transverse amplitude of the first two modes by changing the in-plane loads with other parameters fixed. In Figures 3(a) and 3(b), the bifurcation diagrams of the () versus and are given. Chaotic and periodic motions have been detected, as indicated in Figure 3, which shows the very complex and rich nonlinear dynamics.

(a)

(b)

(c)

(d)

Figure 3: (a) and (b) depict the bifurcation diagrams of versus the first two modes and which are given when the volume fraction index is taken to be 0.5. The Mach number is taken to be 5.0. The magnitude and frequency of the in-plane excitation are taken to be and 161 Hz, respectively. (c) and (d) are partially enlarged details for (a) and (b), respectively. The value of horizontal coordinate is defined as .

As the amplitude of excitation increases within region of , periodic motion occurs for the FGM cylindrical panel. Furthermore, a region of multiperiodic motion exists when the amplitude of excitation varied from to . The system appears to undergo complicated bifurcations where the system is in the instability region. Within the apparently chaotic region where the amplitude of excitation changes from to , there are some regions where the response is periodic. Figures 3(c) and 3(d) are partially enlarged details for Figures 3(a) and 3(b), respectively.

To better understand nonlinear dynamical behaviors, phase portraits, time responses, the Poincare section, and three-dimensional phase portrait are illustrated. Figure 4 illustrates nonlinear dynamic response for the FGM cylindrical panel when the forcing excitation is . The two-dimensional phase plane portraits are given in Figures 4(a) and 4(c) by (). Figures 4(b) and 4(d) denote the time responses on the planes (), respectively. Here subscript “” indicates the first two modes and takes the value . Figures 4(e) and 4(f) represent the three-dimensional phase portrait in space and the Poincare section on plane (), respectively (see [20, 32, 33]). It is known that if the Poincare section displays a finite number of points, the motion can be periodic. Since there is single point in Poincare sections given in Figure 4(f), it can be concluded that there exists periodic oscillations motion for the FGM cylindrical panel in this case. Figure 5 suggests that the multiple-periodic motion of the FGM cylindrical panel occurs when the in-plane excitationchanges to . Figure 6 illustrates the existence of the chaotic motion for this system when the forcing excitation is .

(a)

(b)

(c)

(d)

(e)

(f)

Figure 4: The periodic response of the FGM cylindrical panel occurs when the magnitude of the in-plane excitation is taken to be . The volume fraction index is taken to be 0.5. The Mach number is taken to be 5.0. (a) The phase portrait on plane (); (b) the time response on the planes (); (c) the phase portrait on plane (); (d) the responses on the planes (); (e) three-dimensional phase portrait in space (); (f) the Poincare section on plane ().

(a)

(b)

(c)

(d)

(e)

(f)

Figure 5: The multiple-periodic motion of the FGM cylindrical panel occurs when the magnitude of the in-plane excitation is taken to be . The volume fraction index is taken to be 0.5. The Mach number is taken to be 5.0.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 6: The chaotic motion of the FGM cylindrical panel occurs when the magnitude of the in-plane excitation is taken to be . The volume fraction index is taken to be 0.5. The Mach number is taken to be 5.0.

Figure 7, where in-plane excitation is also taken as the bifurcation parameter, shows the effect of in-plane excitation on the dynamics of the FGM cylindrical panel when the Mach is 10. Compared to the bifurcation diagrams given in Figure 3 whose Mach is 5, we can see that the appearance of the bifurcation point is almost in the same locations. In the instability region, the chaotic and multiperiodic motion may be detected. Figure 8 indicates that the chaotic motion of the FGM cylindrical panel occurs.

(a)

(b)

Figure 7: (a) and (b) depict the bifurcation diagrams of versus and which are given when the Mach number is taken to be 10.0. The value of horizontal coordinate has the same definition as that in Figure 3.

(a)

(b)

(c)

(d)

(e)

(f)

Figure 8: The chaotic motion of the FGM cylindrical panel occurs when the magnitude of the in-plane excitation is taken to be . The volume fraction index is taken to be 0.5. The Mach number is taken to be 10.0.

4. Conclusions

The nonlinear dynamics of an FGM cylindrical panel under a hypersonic flow are presented. The material properties are graded continuously throughout the thickness of the panel according to the power law function and are temperature dependent. A third-order piston theory is applied for the hypersonic aerodynamic load. The Von Karman larger deflection theory in conjunction with energy approach is used to obtain the equations of motion. The bifurcation diagrams, phase portraits, time responses, and the Poincare section are employed to understand the periodic and chaotic motions of the cylindrical panel. It is obtained that, with the change of in-plane load parameter, the different nonlinear dynamics of the FGM panel occur. Different nonlinear dynamic behaviors alternate from stable periodic motion to instable chaotic motions in hypersonic flow. In addition, from the phase portraits and time responses, it can be seen that when the chaotic motion occurs, their phase portraits resemble each other in appearance for different Mach.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors acknowledge the financial support of the National Natural Science Foundation of China through Grants nos. 11272063, 11102226, and 11472298, the Science Foundation of Beijing Municipal Education Commission through Grant no. cit&tcd201304112, and Foundation of Tianjin City through Grant no. 13JCQNJC04400.

References

E. Efraim and M. Eisenberger, “Exact vibration analysis of variable thickness thick annular isotropic and FGM plates,” Journal of Sound and Vibration, vol. 299, no. 4-5, pp. 720–738, 2007.

R. A. Arciniega and J. N. Reddy, “Large deformation analysis of functionally graded shells,” International Journal of Solids and Structures, vol. 44, no. 6, pp. 2036–2052, 2007.

X. Guo and C. Mei, “Application of aeroelastic modes on nonlinear supersonic panel flutter at elevated temperatures,” Computers and Structures, vol. 84, no. 24-25, pp. 1619–1628, 2006.

L. K. Abbas, X. Rui, P. Marzocca, M. Abdalla, and R. de Breuker, “A parametric study on supersonic/hypersonic flutter behavior of aero-thermo-elastic geometrically imperfect curved skin panel,” Acta Mechanica, vol. 222, no. 1-2, pp. 41–57, 2011.

T. Prakash and M. Ganapathi, “Supersonic flutter characteristics of functionally graded flat panels including thermal effects,” Composite Structures, vol. 72, no. 1, pp. 10–18, 2006.

K.-J. Sohn and J.-H. Kim, “Structural stability of functionally graded panels subjected to aero-thermal loads,” Composite Structures, vol. 82, no. 3, pp. 317–325, 2008.

H. H. Ibrahim, M. Tawfik, and M. Al-Ajmi, “Non-linear panel flutter for temperature-dependent functionally graded material panels,” Computational Mechanics, vol. 41, no. 2, pp. 325–334, 2008.

H. H. Ibrahim, M. Tawfik, and M. Al-Ajmi, “Thermal buckling and nonlinear flutter behavior of functionally graded material panels,” Journal of Aircraft, vol. 44, no. 5, pp. 1610–1618, 2007.

T. Prakash, M. K. Singha, and M. Ganapathi, “A finite element study on the large amplitude flexural vibration characteristics of FGM plates under aerodynamic load,” International Journal of Non-Linear Mechanics, vol. 47, no. 5, pp. 439–447, 2012.

M. Hosseini and S. A. Fazelzadeh, “Aerothermoelastic post-critical and vibration analysis of temperature-dependent functionally graded panels,” Journal of Thermal Stresses, vol. 33, no. 12, pp. 1188–1212, 2010.

S. L. Lee and J. H. Kim, “Thermal post-buckling and limit-cycle oscillation of functionally graded panel with structural damping in supersonic airflow,” Composite Structures, vol. 91, no. 2, pp. 205–211, 2009.

C. Mey, K. Abdel-Motagaly, and R. Chen, “Review of nonlinear panel flutter at supersonic and hypersonic speeds,” Applied Mechanics Reviews, vol. 52, no. 10, pp. 321–332, 1999.

M. Amabili, Nonlinear Vibrations and Stability of Shells and Plates, Cambridge University Press, New York, NY, USA, 2008.

P. Marzocca, S. A. Fazelzadeh, and M. Hosseini, “A review of nonlinear aero-thermo-elasticity of functionally graded panels,” Journal of Thermal Stresses, vol. 34, no. 5-6, pp. 536–568, 2011.

M. Hosseini, S. A. Fazelzadeh, and P. Marzocca, “Chaotic and bifurcation dynamic behavior of functionally graded curved panels under aero-thermal loads,” International Journal of Bifurcation and Chaos, vol. 21, no. 3, pp. 931–954, 2011.

L. Librescu, P. Marzocca, and W. A. Silva, “Supersonic/hypersonic flutter and postflutter of geometrically imperfect circular cylindrical panels,” Journal of Spacecraft and Rockets, vol. 39, no. 5, pp. 802–812, 2002.

S. Mahmoudkhani, H. Haddadpour, and H. M. Navazi, “Supersonic flutter prediction of functionally graded conical shells,” Composite Structures, vol. 92, no. 2, pp. 377–386, 2010.

P. P. Friedmann, J. J. McNamara, B. J. Thuruthimattam, and I. Nydick, “Aeroelastic analysis of hypersonic vehicles,” Journal of Fluids and Structures, vol. 19, no. 5, pp. 681–712, 2004.

Y. S. Touloukian, Thermophysical Properties of High Temperature Solid Materials, Mac Millian, New York, NY, USA, 1967.

W. Zhang, Y. X. Hao, and J. Yang, “Nonlinear dynamics of FGM circular cylindrical shell with clamped-clamped edges,” Composite Structures, vol. 94, no. 3, pp. 1075–1086, 2012.

C. Mei, K. Abdel-Moltagaly, and R. Chen, “Review of nonlinear panel flutter at supersonic and hypersonic speeds,” Applied Mechanics Reviews, vol. 52, no. 10, pp. 321–332, 1999.

M. J. Lighthill, “Oscillating airfoils at high Mach number,” Journal of the Aeronautical Sciences, vol. 20, no. 6, pp. 402–406, 1953.

H. Ashley and G. Zartarian, “Piston theory—a new aerodynamic tool for the aeroelastician,” vol. 23, pp. 1109–1118, 1956.

D. Cao and N. Zhao, “Active control of supersonic/hypersonic aeroelastic flutter for a two-dimensional airfoil with flap,” Science China Technological Sciences, vol. 54, no. 8, pp. 1943–1953, 2011.

J. N. Reddy, Mechanics of Laminated Composite Plates and Shells: Theory and Analysis, CRC Press, 2004.

R. Kadoli and N. Ganesan, “Buckling and free vibration analysis of functionally graded cylindrical shells subjected to a temperature-specified boundary condition,” Journal of Sound and Vibration, vol. 289, no. 3, pp. 450–480, 2006.

A. Bhimaraddi, “Large amplitude vibrations of imperfect antisymmetric angle-ply laminated plates,” Journal of Sound and Vibration, vol. 162, no. 3, pp. 457–470, 1993.

A. Nosier and J. N. Reddy, “A study of non-linear dynamic equations of higher-order shear deformation plate theories,” International Journal of Non-Linear Mechanics, vol. 26, no. 2, pp. 233–249, 1991.

J. N. Reddy, “Analysis of functionally graded plates,” International Journal for Numerical Methods in Engineering, vol. 47, no. 1–3, pp. 663–684, 2000.

H.-S. Shen, “Nonlinear bending response of functionally graded plates subjected to transverse loads and in thermal environments,” International Journal of Mechanical Sciences, vol. 44, no. 3, pp. 561–584, 2002.

L. Librescu, P. Marzocca, and W. A. Silva, “Linear/nonlinear supersonic panel flutter in a high-temperature field,” Journal of Aircraft, vol. 41, no. 4, pp. 918–924, 2004.

Y. X. Hao, W. Zhang, J. Yang, and S. Y. Li, “Nonlinear dynamic response of a simply supported rectangular functionally graded material plate under the time-dependent thermalmechanical loads,” Journal of Mechanical Science and Technology, vol. 25, no. 7, pp. 1637–1646, 2011.

Y. X. Hao, W. Zhang, and J. Yang, “Periodic and chaotic motions of FGM thin plate with two simply supported opposite and two free edges,” International Journal of Bifurcation and Chaos, vol. 21, no. 6, pp. 1737–1753, 2011.