Research Article  Open Access
Flutter Characteristic Study of Composite Sandwich Panel with Functionally Graded Foam Core
Abstract
This paper attempts to investigate the flutter characteristic of sandwich panel composed of laminated facesheets and a functionally graded foam core. The macroscopic properties of the foam core change continuously along this direction parallel to the facesheet lamina. The model used in the study is a simple sandwich panelwing clamped at the root, with three simple types of grading strategies for FGM core: linear grading strategy in the chordwise direction, linear grading strategy in the spanwise direction, and bilinear grading of properties of foam core across the panel. The results show that use of FGM core has the potential to increase the flutter speed of the sandwich panel. Finally, a minimum weight design of composite sandwich panel with lamination parameters of facesheet and density distribution of foam core as design variables is conducted using particle swarm optimization (PSO).
1. Introduction
Aeroelastic phenomena are the result of the mutual interaction of elastic and aerodynamic forces, the occurrence of which during the flight can be destructive and lead to the loss of the aircraft [1]. The stiffness and mass distribution of the wing structure both have effects on the aeroelastic properties. Functionally graded materials have continuously varying properties by spatially varying the distribution of two (or more) materials. FGMs enable changes in structural stiffness and mass distribution without necessarily requiring a geometric change in the structural geometry, so FGM can be a novel aeroelastic tailoring concept [2]. The idea of using functionally graded materials for aeroelastic tailoring is not totally new. For highspeed aircraft, many researchers [3–6] have studied on the thermal flutter analysis of functionally graded materials (FGMs) plates with mechanical properties varying through the thickness. This is because the distribution of temperature is important in affecting the thermal flutter behavior, so FGMs have been shown in these papers to improve the aerothermoelastic panel flutter boundaries [2]. For lowspeed structures, Kuo [7] investigated the flutter analysis of composite plates with variable fiber spacing. Librescu and Maalawi [8] researched the torsional divergence stability of a cantilevered wing using material fiber volume fraction grading of composite materials to optimize the material distribution. Dunning et al. [2] studied the efficiency of FGM for the aeroelastic tailoring of a metallic cantilever platelike wing using both material grading and geometric grading and concluded that an FGM can be used to improve the performance of a platelike wing compared with a uniform material design. Dunning et al. fill a gap in the literature pertaining to FGMbased aeroelastic tailoring by considering a metallic wing in lowspeed (subsonic) conditions.
Based on author’s previous research [9], this paper considers the lowspeed flutter behavior of a clamped platewing, whose structure is composed of laminated facesheets and a functionally graded foam core. The present study aims to consider the flutter characteristic of sandwich panel with the macroscopic properties of the foam core changing continuously along this direction parallel to the facesheet lamina. The model used in the study is a simple sandwich panel clamped at the root, with three simple types of grading strategies for FGM core according to [2]: linear grading strategy in the chordwise direction, linear grading strategy in the spanwise direction, and bilinear grading of properties of foam core across the panel. For a typical sandwich panel with a thick core (), the bending and transverse shear stiffness of the facesheet itself can be neglected. In the present study, composite sandwich stiffness matrix is calculated using lamination parameters of composite facesheets for a thick sandwich panel, and the flutter speed of composite sandwich panel is studied using lamination parameters, with influence of foam density distribution and shear modulus distribution on flutter speed considered.
2. Equivalent Stiffness Matrix of Sandwich
In a cantilevered composite sandwich panel (1.2 m × 0.4 m) with one end fixed, thicknesses of facesheet and foam core are 0.0005 m and 0.01 m. The material properties of T300/QY8911 used for composite facesheet are listed in Table 1. Poisson’s ratio of isotropic material property of foam core is .

The sandwich stiffness matrices can be written in terms of the facesheet and core stiffness as follows [10]:
For composite facesheet, the membrane stiffness matrix can be expressed in terms of lamination parameters and material stiffness invariants:
The material stiffness invariants are calculated as follows:where , , , , and .
The membrane lamination parameters are given by the following integrals:
The lamination parameters cannot be chosen arbitrarily. For a laminate of 0, 90, 45, and −45 degree fiber angles, . Considering the dependency, the feasible domain for the inplane lamination parameters is known to be defined by [11]
3. Gradient Material Properties of Foam Core
For the density gradient foam material, to calculate Young’s modulus of the closedcell foam material, Xiao et al. [12] gave a modified equation based on the Gibson equation [13]:where and are the Young’s modulus and mass density of closedcell foam, respectively; and are the Young’s modulus and mass density of the solid material of the cell wall, respectively. The material properties of the solid material are assumed: = 100 MPa and = 80 kg/m^{3}.
The flutter speed is always affected by the spanwise bending stiffness and chordwise torsional stiffness for a wing panel. In order to research the flutter speed with grading bending and torsional stiffness, three types of grading strategies [2] are adopted to consider the grading transverse shear of foam core for composite sandwich flutter characteristics.
Grading Strategy 1 (GS1, shown as Figure 2(a)). Density grading in the chordwise direction is as follows:where and are the mass densities of closedcell foam of the leading edge and trailing edge, respectively; denotes the width of the panel.
Grading Strategy 2 (GS2, shown as Figure 2(b)). Density grading in the spanwise direction is as follows:where and are the mass densities of closedcell foam of the root edge and tip edge, respectively; denotes the length of the panel.
Grading Strategy 3 (GS3, shown as Figure 2(c)). Density diagonal grading across the plane is as follows:where and are the mass densities of closedcell foam of upper and lower bounds, respectively.
Submitting for into (6), Young’s modulus of DG closedcell foam is obtained.
4. Flutter Speed of Composite Sandwich Panel
4.1. Aeroelastic Solution Methodology
PKmethod method is adopted for flutter speed calculation [14]. The principal advantage of the PKmethod is that it produces results directly for given values of velocity where , , and denote modal mass matrix, modal stiffness matrix, and modal damping matrix, respectively. , , , and denote the eigenvalue, density, reference chord length, and velocity. and are modal aerodynamic force matrix and modal aerodynamic damping matrix, which are both functions of Mach number (Ma) and reduced frequency (). The divergence speed can be evaluated when tends to zero.
4.2. Flutter Speed Contours for Various Laminate Configurations
Keeping the total mass of the foam core constant and varying the densities of and , the influence of the density gradient changes of the foam cores on the flutter speed is investigated. To describe the changing range of the mass density, the gradient factor (GF) is defined as the difference between and in Table 2.
 
density of 50 kg/m^{3} (nongradient). 
4.2.1. Flutter Velocity Contour for Orthotropic Panel with Gradient Foam
Orthotropic sandwich panel ( for facesheets) is investigated firstly, shown in Figure 3. The bendingtorsional flutter of the orthotropic panel is due to coupling between the second bending mode and the first torsional mode, with the second mode as the critical vibration mode. The histories of the damping and frequency of the first four mode for the quasiisotropic panel () are illustrated in Figure 4. The critical flutter mode is the second mode.
Figure 5 gives the flutter velocity contours with various gradient change for orthotropic sandwich panel. For the three types of grading strategies, it is concluded that the flutter speed of the sandwich panel increases when the gradient changes of the foam cores (GF) increase, though the total masses of the cores keep constant. Also, the velocity contours are continuous within the feasible lamination parameter region, all with the second mode as the critical vibration mode. Table 3 gives the average increases of the flutter speeds for various grading strategies and grading factors compared with the sandwich panel with no gradient change. Moreover, density grading strategy in the chordwise direction (GS1) has the most obvious influence on the flutter speed, and density diagonal grading across the panel (GS3) has less influence compared with the other two grading strategies. And the diagonal density grading strategy can lead to minor decrease of the flutter speed.

4.2.2. Flutter Velocity Contour for Nonorthotropic Panel with Gradient Foam
When is positive, the bendingtorsional flutter mode still remains the same as the orthotropic sandwich panel [9]. When is negative, there may be another flutter mode within part of the lamination parameter region (shown as Figure 6), with the third mode as the critical vibration mode. Two sets of lamination parameters in Figure 6 are chosen to examine the effect of laminate configurations on the flutter modes: Point B: , , , Point C: , , .
Figure 7 shows the histories of the damping and frequency of the first four modes for the two cases.
Figure 8 details the flutter speed of the two sets with various gradient factors. For the two sets, density grading strategy in the chordwise direction (GS1) has the most obvious influence on the flutter speed. For Point C, the flutter speed of the sandwich panel increases when the gradient changes of the foam cores (GF) increase with the total masses of the cores keeping constant for GS1 and GS2, while, for point B, there occurs a decrease of the flutter speed compared with the nongradient panel for GS1 and GS2. Gradient change has the potential to increase the flutter speed. And the diagonal density grading (GS3) strategy can lead to minor decrease of the flutter speed for both sets. This is because the flutter characteristic is dependent on several other factors: wing sweep, aspect ratio, and the stiffness distribution. However, it shows that the proper gradient of the foam has the potential the flutter speed of the sandwich pane.
5. Numerical Optimization Example
5.1. Numerical Example Description
Minimum weight design of cantilevered sandwich panel as shown in Figure 1 with the flutter and divergence velocity constraints is carried out. The design variables include thickness and membrane lamination parameters of facesheet (), thickness, and density distribution of foam core () because the density grading strategy in the chordwise direction (GS1) has the most obvious influence on the flutter speed, so the GS1 is chosen as density grading strategy for foam core. The critical speed is assumed not less than 130 m/s. The constraints of lamination parameters , , and are described in (11), and the bound constraints of the , , , and are assumed.
(a) Density grading in the chordwise direction
(b) Density grading in the spanwise direction
(c) Density diagonal grading across the plane
(a) GF = 1 (GS1)
(b) GF = 2 (GS1)
(c) GF = 3 (GS1)
(d) GF = 1 (GS2)
(e) GF = 2 (GS2)
(f) GF = 3 (GS2)
(g) GF = 1 (GS3)
(h) GF = 2 (GS3)
(i) GF = 3 (GS3)
(a) Point B
(b) Point C
(a) Point B
(b) Point C
The optimization problem can be stated as follows:
5.2. Optimization Algorithm
Because of the ability to deal with the continuous global optimization problem with a nonlinear objective function, particle swarm optimization (PSO) is an evolutionary global algorithm and has become more and more popular. PSO was first proposed by Kennedy and Eberhart [15, 16]. It is observed that a swarm of birds or insects search for food in a very typical manner. If one member of the swarm finds a desirable path to go, the rest of the particles will follow quickly. Each particle searches for the best in its locality with learning from its own experience. Additionally, each member learns from the others, typically from the best performer among them. PSO has been successfully applied to some engineering and structural problems.
The basic steps in the PSO algorithm are as follows.
Step 1. Initialize the swarm with random position values and random initial velocities.
Step 2. Determine the velocity vector for each particle in the swarm using the knowledge of the best position obtained by each particle and the swarm as a whole and also the previous position of each particle in the swarm.
Step 3. Modify the current position of each particle using the velocity vector and the previous position of each particle.
Step 4. Repeat from Step 2 until the stop criterion is achieved.
The velocity vector of each particle is calculated as follows:where the superscript denotes the particle and the subscript denotes the iteration number; denotes the velocity and denotes the position; and are uniformly distributed random numbers in the interval ; and are the acceleration constants; is the inertia weight; is the best position attained by the particle in the swarm so far and is the global best position attained by the swarm at iteration .
The position of each particle at iteration is calculated using the formula:
5.3. Results
According to (911) of [9], with the optimum facesheet thickness (0.125 mm) and optimum lamination parameters ), the corresponding layer thickness of facesheet can be obtained (). And the density distribution of the foam core is kg/m^{3}. So the optimum sandwich configuration is obtained as , shown in Table 4. The speeds of flutter and divergence of the optimum sandwich configuration are 130.00 m/s and 213.95 m/s, respectively, shown as Figure 9.
 
denotes a 45 degree ply with 0.05775 mm, denotes a −45 degree ply with 0.06725 mm, and denotes a core with thickness of 5 mm. means that the type of sandwich is symmetric. 
6. Conclusions
In this paper, the flutter speed of a cantilevered sandwich panel is studied using lamination parameters of composite facesheet, with influence of density gradient on sandwich panel considered. The results show that the flutter speed of the sandwich panel has the potential to increase when the proper gradient changes of the foam cores (GF) increase for chordwise strategy and spanwise strategy with the total mass of the core constant. A minimum weight design of composite sandwich panel with lamination parameters of facesheet and density distribution of foam core as design variables was conducted using PSO. The results show that the PSO algorithm can be effective for the lamination parametersbased optimizaiton problem.
The paper conducts a series of flutter analysis for a range of functionally graded core of composite sandwich. The previous studies in this area have already shown that the flutter characteristics are sensitive to the wing geometry (e.g., aspect ratio and sweep) as well as material properties and that the conclusion drawn in this paper is a little simplistic. Further research will be focused on careful analysis considering other key parameters based on the existing literature and offer a more comprehensive discussion.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The support for this work by the National Natural Science Foundation of China (11402204) is gratefully acknowledged.
References
 A. Attaran, D. L. Majid, S. Basri, A. S. Mohd Rafie, and E. J. Abdullah, “Structural optimization of an aeroelastically tailored composite flat plate made of woven fiberglass/epoxy,” Aerospace Science and Technology, vol. 15, no. 5, pp. 393–401, 2011. View at: Publisher Site  Google Scholar
 P. D. Dunning, B. K. Stanford, H. A. Kim, and C. V. Jutte, “Aeroelastic tailoring of a plate wing with functionally graded materials,” Journal of Fluids and Structures, vol. 51, pp. 292–312, 2014. View at: Publisher Site  Google Scholar
 T. Prakash and M. Ganapathi, “Asymmetric flexural vibration and thermoelastic stability of FGM circular plates using finite element method,” Composites Part B: Engineering, vol. 37, no. 78, pp. 642–649, 2006. View at: Publisher Site  Google Scholar
 H. H. Ibrahim, M. Tawfik, and M. AlAjmi, “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
 H. H. Ibrahim, M. Tawfik, and M. AlAjmi, “Nonlinear panel flutter for temperaturedependent functionally graded material panels,” Computational Mechanics, vol. 41, no. 2, pp. 325–334, 2008. View at: Publisher Site  Google Scholar
 H. M. Navazi and H. Haddadpour, “Aerothermoelastic stability of functionally graded plates,” Composite Structures, vol. 80, no. 4, pp. 580–587, 2007. View at: Publisher Site  Google Scholar
 S.Y. Kuo, “Flutter of rectangular composite plates with variable fiber pacing,” Composite Structures, vol. 93, no. 10, pp. 2533–2540, 2011. View at: Publisher Site  Google Scholar
 L. Librescu and K. Y. Maalawi, “Material grading for improved aeroelastic stability in composite wings,” Journal of Mechanics of Materials and Structures, vol. 2, no. 7, pp. 1381–1394, 2007. View at: Publisher Site  Google Scholar
 P. Jin, B. Song, X. Zhong, T. Yu, and F. Xu, “Aeroelastic tailoring of composite sandwich panel with lamination parameters,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 230, no. 1, pp. 105–117, 2016. View at: Publisher Site  Google Scholar
 V. Balabanov, O. Weckner, M. Epton, G. Mabson, E. Cregger, and A. W. Blom, “Optimal design of a composite sandwich structure using lamination parameters,” in Proceedings of the 53rd AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics and Materials Conference, Honolulu, Hawaii, USA, April 2012. View at: Google Scholar
 H. Fukunaga and G. N. Vanderplaats, “Stiffness optimization of orthotropic laminated composites using lamination parameters,” AIAA journal, vol. 29, no. 4, pp. 641–646, 1991. View at: Publisher Site  Google Scholar
 D. B. Xiao, L. Mu, and G. P. Zhao, “The influence of correlating material parameters of gradient foam core on free vibration of sandwich panel,” Composites Part B: Engineering, vol. 77, pp. 153–161, 2015. View at: Publisher Site  Google Scholar
 L. J. Gibson and M. F. Ashby, Cellular Solids: Structure and Properties, Cambridge University Press, Cambridge, UK, 2nd edition, 1997.
 W. P. Rodden and E. H. Johnson, MSC/Nastran Aeroelastie Analysis User's Guide V68, MSC, Software Corporation, Los Angeles Calif, USA, 1994.
 J. Kennedy and R. Eberhart, “Particle swarm optimization,” in Proceedings of the IEEE International Conference on Neural Networks, pp. 1942–1948, Perth, Australia, December 1995. View at: Google Scholar
 R. Eberhart and J. Kennedy, “New optimizer using particle swarm theory,” in Proceedings of the 6th International Symposium on Micro Machine and Human Science (MHS '95), pp. 39–43, Nagoya, Japan, October 1995. View at: Google Scholar
Copyright
Copyright © 2016 Peng Jin and Xiaoping Zhong. 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.