Nonlinear Vibrations of FGM Cylindrical Panel with Simply Supported Edges in Air Flow
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.
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 . 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.  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 . 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  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.  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  used the numerical and analytical methodologies to study postcritical and vibration behaviors of the FGM panels in a supersonic air flow. Ibrahim et al.  analyzed nonlinear flutter and thermal buckling of an FGM panel under the combined effect of elevated temperature conditions and aerodynamic loading. Lee and Kim  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.  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 . Numerous studies on the cylindrical shells focused on their flutter. Marzocca et al.  reviewed the corresponding researches that analyse the dynamic behavior of curved and flat panels exposed to supersonic flow fields. Hosseini et al.  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.  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.  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 . 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.
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  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. , 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. . 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  and later used as an aeroelastic tool by Ashley and Zartarian , 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  and Cao and Zhao  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  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 : 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  and Kadoli and Ganesan  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 , 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 ) 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.
In this study, Al2O3 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 . 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. . 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.
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 .
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.
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.
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.
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
M. Amabili, Nonlinear Vibrations and Stability of Shells and Plates, Cambridge University Press, New York, NY, USA, 2008.View at: Publisher Site | MathSciNet
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
M. J. Lighthill, “Oscillating airfoils at high Mach number,” Journal of the Aeronautical Sciences, vol. 20, no. 6, pp. 402–406, 1953.View at: Google Scholar | MathSciNet
H. Ashley and G. Zartarian, “Piston theory—a new aerodynamic tool for the aeroelastician,” vol. 23, pp. 1109–1118, 1956.View at: Google Scholar | MathSciNet
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
J. N. Reddy, “Analysis of functionally graded plates,” International Journal for Numerical Methods in Engineering, vol. 47, no. 1–3, pp. 663–684, 2000.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Publisher Site | Google Scholar