Flutter Characteristic Study of Composite Sandwich Panel with Functionally Graded Foam Core

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 panel-wing clamped at the root, with three simple types of grading strategies for FGMcore: (1) linear grading strategy in the chord-wise direction, (2) linear grading strategy in the span-wise direction, and (3) 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).


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 high-speed aircraft, many researchers [3][4][5][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 low-speed 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 plate-like wing using both material grading and geometric grading and concluded that an FGM can be used to improve the performance of a plate-like wing compared with a uniform material design.Dunning et al. fill a gap in the literature pertaining to FGM-based aeroelastic tailoring by considering a metallic wing in low-speed (subsonic) conditions.
Based on author's previous research [9], this paper considers the low-speed flutter behavior of a clamped plate-wing, 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]: (1) linear grading strategy in the chord-wise direction, (2) linear grading strategy in the 2 International Journal of Aerospace Engineering span-wise direction, and (3) bilinear grading of properties of foam core across the panel.For a typical sandwich panel with a thick core ( face / core ≪ 1), 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.

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 V core = 0.25.
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 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,  4 = 0. Considering the dependency, the feasible domain for the inplane lamination parameters is known to be defined by [11]

Gradient Material Properties of Foam Core
For the density gradient foam material, to calculate Young's modulus of the closed-cell 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 closed-cell 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  s = 80 kg/m 3 .The flutter speed is always affected by the span-wise bending stiffness and chord-wise 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 chord-wise direction is as follows: where   and   are the mass densities of closed-cell 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 span-wise direction is as follows: where   and   are the mass densities of closed-cell 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 closed-cell foam of upper and lower bounds, respectively.Submitting for   into (6), Young's modulus of DG closed-cell foam is obtained.

Flutter Speed of Composite Sandwich Panel
4.1.Aeroelastic Solution Methodology.PK-method method is adopted for flutter speed calculation [14].The principal advantage of the PK-method 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.

Flutter Speed Contours for Various Laminate Configurations.
Keeping the total mass of the foam core constant and varying the densities of  min core and  max core , 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  min core and  max core in Table 2.

Flutter Velocity Contour for Orthotropic Panel with
Gradient Foam.Orthotropic sandwich panel ( 3 = 0 for facesheets) is investigated firstly, shown in Figure 3.The bending-torsional 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 quasi-isotropic panel ( 1 =  2 =  3 = 0) 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 chord-wise 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.

Flutter Velocity Contour for Nonorthotropic Panel with
Gradient Foam.When  3 is positive, the bending-torsional flutter mode still remains the same as the orthotropic sandwich panel [9].When  3 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: 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 chord-wise 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  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.carried out.The design variables include thickness and membrane lamination parameters of facesheet ( face ,  1 ,  2 ,  3 ), thickness, and density distribution of foam core ( core ,  min core ,  max core ) because the density grading strategy in the chord-wise 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  1 ,  2 , and  3 are described in (11) [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.

Numerical Optimization Example
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; V denotes the velocity and  denotes the position;  1 and  2 are uniformly distributed random numbers in the interval [0, 1];  1 and  2 are the acceleration constants;  is the inertia weight;   is the best position attained by the particle  in the swarm so far and   −1 is the global best position attained by the swarm at iteration  − 1.The position of each particle at iteration  is calculated using the formula:

Figure 1 :
Figure7shows the histories of the damping and frequency of the first four modes for the two cases.Figure8details the flutter speed of the two sets with various gradient factors.For the two sets, density grading strategy in the chord-wise 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

Table 3 :
Average increases of the flutter speeds for various GS and GF compared with the sandwich panel with no gradient change.
, and the bound constraints of the  face ,  core ,  min core , and  max core are assumed.The optimization problem can be stated as follows: min: mass face + mass core s.t.V critical = min [V  , V  ] ≥ 130 m/s;  core ,  face ,  min core ,  max core ,  1 ,  2 ,  3 .