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

Authors: Peng Jin; Xiaoping Zhong
Journal: International Journal of Aerospace Engineering
 (2016)

---

## 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 panel-wing clamped at the root, with three simple types of grading strategies for FGM core:(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).


---

## Body

## 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 high-speed 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 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 span-wise direction, and (3) bilinear grading of properties of foam core across the panel. For a typical sandwich panel with a thick core (tface/tcore≪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.


## 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 Table1. Poisson’s ratio of isotropic material property of foam core is vcore=0.25.Table 1
Material property of T300/QY891.
E 11 (GPa) E 22 (GPa) μ 12 G 12 (GPa) ρ f a c e (kg/m3) t p l y (mm) 127.56 13.03 0.3 6.41 1577.76 0.125The sandwich stiffness matrices can be written in terms of the facesheet and core stiffness as follows [10]:(1)ASandwich=2AFace,Dsandwich=2AFacetcore+tface22,GSandwich=Gcore.For composite facesheet, the membrane stiffness matrix can be expressed in terms of lamination parameters and material stiffness invariants:(2)
A
F
a
c
e
=
A
11
A
22
A
12
A
66
A
16
A
26
=
t
f
a
c
e
A
11
∗
A
22
∗
A
12
∗
A
66
∗
A
16
∗
A
26
∗
=
t
f
a
c
e
U
1
V
1
V
2
U
1
-
V
1
V
2
U
4
0
-
V
2
U
5
0
-
V
2
0
V
3
2
V
4
0
V
3
2
V
4
1
U
2
U
3
.The material stiffness invariants are calculated as follows:(3)U1U2U3U4U5=18323440-401-21-4161-41-214Q11Q12Q22Q66,where Q11=E11/1-ν12ν21, Q22=E22/1-ν12ν21, Q21=Q12, Q66=G12, and ν21=ν12E22/E11.The membrane lamination parameters are given by the following integrals:(4)V1234=1h∫-h/2h/2cos⁡2θcos⁡4θsin⁡2θsin⁡4θdz.The lamination parameters cannot be chosen arbitrarily. For a laminate of 0, 90, 45, and −45 degree fiber angles,V4=0. Considering the dependency, the feasible domain for the in-plane lamination parameters is known to be defined by [11](5)2V1-V2-1≤0,2V3+V2-1≤0.


## 3. 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]:(6)Efc=1-ρfcρs2ρfcρs2+ρfcρs2Es,where Efc and ρfc are the Young’s modulus and mass density of closed-cell foam, respectively; Es and ρs 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: Es = 100 MPa and ρs = 80 kg/m3.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:(7)ρfc=ρl-xρl-ρtb,where ρl and ρt are the mass densities of closed-cell foam of the leading edge and trailing edge, respectively; b 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:(8)ρfc=ρr-yρr-ρta,where ρr and ρt are the mass densities of closed-cell foam of the root edge and tip edge, respectively; a denotes the length of the panel.Grading Strategy  3 (GS3, shown as Figure 2(c)). Density diagonal grading across the xy plane is as follows:(9)ρfc=b-xba-yaρu+yaρl+xba-yaρl+yaρu,where ρu and ρl are the mass densities of closed-cell foam of upper and lower bounds, respectively.Submitting forρfc into (6), Young’s modulus of DG closed-cell foam is obtained.


## 4. 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 (10)Mhhp2+Bhh-14ρcVQhhIkp+khh-12ρV2QhhRuh=0,where Mhh, khh, and Bhh denote modal mass matrix, modal stiffness matrix, and modal damping matrix, respectively. p, ρ, c, and V denote the eigenvalue, density, reference chord length, and velocity. QhhR and QhhI are modal aerodynamic force matrix and modal aerodynamic damping matrix, which are both functions of Mach number (Ma) and reduced frequency (k). The divergence speed VD can be evaluated when k 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ρcoremin and ρcoremax, 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 ρcoremin and ρcoremax in Table 2.Table 2
Gradient factor.
Gradient factor (GF) 0 1 2 3 Density variation range [ρcoremin,ρcoremax] 50,50 a [ 40,60 ] [ 30,70 ] [ 20,80 ] ⁢ C o n s t a n t   a density of 50 kg/m3 (nongradient).


#### 4.2.1. Flutter Velocity Contour for Orthotropic Panel with Gradient Foam

Orthotropic sandwich panel (V3=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 (V1=V2=V3=0) are illustrated in Figure 4. The critical flutter mode is the second mode.Figure5 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.Table 3
Average increases of the flutter speeds for various GS and GF compared with the sandwich panel with no gradient change.
GF1 GF2 GF3 GS1 2.33% 4.07% 26.58% GS2 1.94% 3.69% 5.32% GS3 −0.10% −0.58% −1.45%


#### 4.2.2. Flutter Velocity Contour for Nonorthotropic Panel with Gradient Foam

WhenV3 is positive, the bending-torsional flutter mode still remains the same as the orthotropic sandwich panel [9]. When V3 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:  V1=-0.4, V2=-0.2, V3=-0.6,Point C:  V1=0.4, V2=-0.2, V3=-0.6.Figure7 shows the histories of the damping and frequency of the first four modes for the two cases.Figure8 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 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.


## 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 (10)Mhhp2+Bhh-14ρcVQhhIkp+khh-12ρV2QhhRuh=0,where Mhh, khh, and Bhh denote modal mass matrix, modal stiffness matrix, and modal damping matrix, respectively. p, ρ, c, and V denote the eigenvalue, density, reference chord length, and velocity. QhhR and QhhI are modal aerodynamic force matrix and modal aerodynamic damping matrix, which are both functions of Mach number (Ma) and reduced frequency (k). The divergence speed VD can be evaluated when k 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ρcoremin and ρcoremax, 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 ρcoremin and ρcoremax in Table 2.Table 2
Gradient factor.
Gradient factor (GF) 0 1 2 3 Density variation range [ρcoremin,ρcoremax] 50,50 a [ 40,60 ] [ 30,70 ] [ 20,80 ] ⁢ C o n s t a n t   a density of 50 kg/m3 (nongradient).


### 4.2.1. Flutter Velocity Contour for Orthotropic Panel with Gradient Foam

Orthotropic sandwich panel (V3=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 (V1=V2=V3=0) are illustrated in Figure 4. The critical flutter mode is the second mode.Figure5 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.Table 3
Average increases of the flutter speeds for various GS and GF compared with the sandwich panel with no gradient change.
GF1 GF2 GF3 GS1 2.33% 4.07% 26.58% GS2 1.94% 3.69% 5.32% GS3 −0.10% −0.58% −1.45%


### 4.2.2. Flutter Velocity Contour for Nonorthotropic Panel with Gradient Foam

WhenV3 is positive, the bending-torsional flutter mode still remains the same as the orthotropic sandwich panel [9]. When V3 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:  V1=-0.4, V2=-0.2, V3=-0.6,Point C:  V1=0.4, V2=-0.2, V3=-0.6.Figure7 shows the histories of the damping and frequency of the first four modes for the two cases.Figure8 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 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.


## 4.2.1. Flutter Velocity Contour for Orthotropic Panel with Gradient Foam

Orthotropic sandwich panel (V3=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 (V1=V2=V3=0) are illustrated in Figure 4. The critical flutter mode is the second mode.Figure5 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.Table 3
Average increases of the flutter speeds for various GS and GF compared with the sandwich panel with no gradient change.
GF1 GF2 GF3 GS1 2.33% 4.07% 26.58% GS2 1.94% 3.69% 5.32% GS3 −0.10% −0.58% −1.45%


## 4.2.2. Flutter Velocity Contour for Nonorthotropic Panel with Gradient Foam

WhenV3 is positive, the bending-torsional flutter mode still remains the same as the orthotropic sandwich panel [9]. When V3 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:  V1=-0.4, V2=-0.2, V3=-0.6,Point C:  V1=0.4, V2=-0.2, V3=-0.6.Figure7 shows the histories of the damping and frequency of the first four modes for the two cases.Figure8 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 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 Figure1 with the flutter and divergence velocity constraints is carried out. The design variables include thickness and membrane lamination parameters of facesheet (tface,V1,V2,V3), thickness, and density distribution of foam core (tcore,ρcoremin,ρcoremax) 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 V1, V2, and V3 are described in (11), and the bound constraints of the tface, tcore, ρcoremin, and ρcoremax are assumed.Figure 1
Cantilevered sandwich panel.Figure 2
Material grading across the sandwich panel [2].
(a)
Density grading in the chord-wise direction (b)
Density grading in the span-wise direction (c)
Density diagonal grading across thexy planeFigure 3
Flutter velocity contours for orthotropic sandwich panel with no gradient change (V3=0, GF = 0).Figure 4
V - g and v-f curves of quasi-isotropic panel (V1=V2=V3=0).Figure 5
Flutter velocity contours for orthotropic sandwich panel (V3=0) with various gradient changes.
(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)Figure 6
Flutter velocity contours for nonorthotropic sandwich panel with no gradient change (V3=-0.6, GF = 0).Figure 7
V - g and v-f curves of the two points.
(a)
Point B (b)
Point CFigure 8
Flutter speed of the two sets with various GS and GF.
(a)
Point B (b)
Point CThe optimization problem can be stated as follows:(11)min:massface+masscores.t.vcritical=min⁡vf,vD≥130m/s;2V1-V2-1≤0;2V3+V2-1≤0;tface∈0.000125m,0.000625mtcore∈0.005m,0.01mρcoremax∈20kg/m3,80kg/m3;ρcoremin∈20kg/m3,80kg/m3;ρcoremin≤ρcoremaxdesign  variables:tcore,tface,ρcoremin,ρcoremax,V1,V2,V3.


### 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 Step2 until the stop criterion is achieved.The velocity vector of each particle is calculated as follows:(12)vki=wvk-1i+c1r1pi-xk-1i+c2r2pk-1g-xk-1i,where the superscript i denotes the particle and the subscript k denotes the iteration number; v denotes the velocity and x denotes the position; r1 and r2 are uniformly distributed random numbers in the interval [0,1]; c1 and c2 are the acceleration constants; w is the inertia weight; pi is the best position attained by the particle i in the swarm so far and pk-1g is the global best position attained by the swarm at iteration k-1.The position of each particle at iterationk is calculated using the formula:(13)xki=xk-1i+vki.


### 5.3. Results

According to (9-11) of [9], with the optimum facesheet thickness (0.125 mm) and optimum lamination parameters ([V1,V2,V3]=[0,-1,-0.076]), the corresponding layer thickness of facesheet can be obtained ([n0,n90,n45,n-45]=[0,0,0.05775mm,0.06725mm]). And the density distribution of the foam core is [ρcoremin,ρcoremax]=[20,37.64] kg/m3. So the optimum sandwich configuration is obtained as 450.05775mm/-450.06725mm/foam5mms, 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.Table 4
Optimum design variables and laminate configuration.
t c o r e (mm) t f a c e (mm) ρ c o r e m i n (kg/m3) ρ c o r e m a x (kg/m3) V 1 V 2 V 3 Sandwichconfigurationa 10 0.125 20 37.64 0 −1 −0.076 45 0.05775  m m / - 45 0.06725  m m / f o a m 5  m m s ⁢ a 45 0.05775  m m denotes a 45 degree ply with 0.05775 mm, -450.06725mm denotes a −45 degree ply with 0.06725 mm, and foam5mm denotes a core with thickness of 5 mm. s means that the type of sandwich is symmetric.Figure 9
V - g and v-f curves of the optimum result.


## 5.1. Numerical Example Description

Minimum weight design of cantilevered sandwich panel as shown in Figure1 with the flutter and divergence velocity constraints is carried out. The design variables include thickness and membrane lamination parameters of facesheet (tface,V1,V2,V3), thickness, and density distribution of foam core (tcore,ρcoremin,ρcoremax) 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 V1, V2, and V3 are described in (11), and the bound constraints of the tface, tcore, ρcoremin, and ρcoremax are assumed.Figure 1
Cantilevered sandwich panel.Figure 2
Material grading across the sandwich panel [2].
(a)
Density grading in the chord-wise direction (b)
Density grading in the span-wise direction (c)
Density diagonal grading across thexy planeFigure 3
Flutter velocity contours for orthotropic sandwich panel with no gradient change (V3=0, GF = 0).Figure 4
V - g and v-f curves of quasi-isotropic panel (V1=V2=V3=0).Figure 5
Flutter velocity contours for orthotropic sandwich panel (V3=0) with various gradient changes.
(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)Figure 6
Flutter velocity contours for nonorthotropic sandwich panel with no gradient change (V3=-0.6, GF = 0).Figure 7
V - g and v-f curves of the two points.
(a)
Point B (b)
Point CFigure 8
Flutter speed of the two sets with various GS and GF.
(a)
Point B (b)
Point CThe optimization problem can be stated as follows:(11)min:massface+masscores.t.vcritical=min⁡vf,vD≥130m/s;2V1-V2-1≤0;2V3+V2-1≤0;tface∈0.000125m,0.000625mtcore∈0.005m,0.01mρcoremax∈20kg/m3,80kg/m3;ρcoremin∈20kg/m3,80kg/m3;ρcoremin≤ρcoremaxdesign  variables:tcore,tface,ρcoremin,ρcoremax,V1,V2,V3.


## 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 Step2 until the stop criterion is achieved.The velocity vector of each particle is calculated as follows:(12)vki=wvk-1i+c1r1pi-xk-1i+c2r2pk-1g-xk-1i,where the superscript i denotes the particle and the subscript k denotes the iteration number; v denotes the velocity and x denotes the position; r1 and r2 are uniformly distributed random numbers in the interval [0,1]; c1 and c2 are the acceleration constants; w is the inertia weight; pi is the best position attained by the particle i in the swarm so far and pk-1g is the global best position attained by the swarm at iteration k-1.The position of each particle at iterationk is calculated using the formula:(13)xki=xk-1i+vki.


## 5.3. Results

According to (9-11) of [9], with the optimum facesheet thickness (0.125 mm) and optimum lamination parameters ([V1,V2,V3]=[0,-1,-0.076]), the corresponding layer thickness of facesheet can be obtained ([n0,n90,n45,n-45]=[0,0,0.05775mm,0.06725mm]). And the density distribution of the foam core is [ρcoremin,ρcoremax]=[20,37.64] kg/m3. So the optimum sandwich configuration is obtained as 450.05775mm/-450.06725mm/foam5mms, 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.Table 4
Optimum design variables and laminate configuration.
t c o r e (mm) t f a c e (mm) ρ c o r e m i n (kg/m3) ρ c o r e m a x (kg/m3) V 1 V 2 V 3 Sandwichconfigurationa 10 0.125 20 37.64 0 −1 −0.076 45 0.05775  m m / - 45 0.06725  m m / f o a m 5  m m s ⁢ a 45 0.05775  m m denotes a 45 degree ply with 0.05775 mm, -450.06725mm denotes a −45 degree ply with 0.06725 mm, and foam5mm denotes a core with thickness of 5 mm. s means that the type of sandwich is symmetric.Figure 9
V - g and v-f curves of the optimum result.


## 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 chord-wise strategy and span-wise 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 parameters-based 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.



---

*Source: 7971435-2016-11-10.xml*