Published: 30 June 2016

An analytical model for flutter behavior of composite panels with shape memory alloy fibers

Yongsheng Ren1
Chenggang Du2
Xiaoqing Ma3
1, 2, 3College of Mechanical and Electronic Engineering, Shandong University of Science and Technology, Qindao 266590, China
Corresponding Author:
Yongsheng Ren
Views 119
Reads 58
Downloads 1310

Abstract

An analytical model for predicting the effect of shape memory alloys (SMA) on the flutter behavior of composite panels is developed in the frequency and time domains. The laminated plate theory and piston theory are employed to model the aeroelastic response of the composite panels to aerodynamic loads. A thermo-mechanical constitutive equation of SMA proposed by Brinson et al. is used to calculate the recovery stress of the constrained SMA fibers. The approximate solution is obtained for supersonic flutter analysis of the composite panels based on the Galerkin approach. The parametric study is carried out to display the effect of the actuation temperature, volume fraction, the initial strain of SMA fibers and the length to width ratio of the panels on flutter boundaries.

1. Introduction

The surface panel of the vehicles made of lightweight composite may result in flutter or self-exited oscillation during high-speed flight. Although panel flutter does not indicate structural failure as happened in the aircraft wing, this may cause fatigue of panel and affect its aerodynamic characteristics, leading to poor flight performance. In the recent years, the studies of the flutter suppression using shape memory alloy materials have received much attention. It is only natural that SMA is embedded into the composite panel to increase the aeroelastic stability of vehicles. This is because composite panel is thin thanks to their high stiffness/weight ratio. On the other hand, SMA has a unique ability known as shape memory effect which can result in the large tensile recovery stress if the SMA is restrained. Large in-plane forces can be induced in the composite panel at temperatures higher than the austenite start temperature which increase the composite panel stiffness.

Weisshaar et al. [1-2] showed that the flutter boundary of the panel with partial SMA layer is significantly increased. Park et al. [3] and Duan et al. [4] investigated the control of the thermal post-buckling deflection, critical temperature and critical dynamic pressure by using the SMA. Ostachowicz et al. [5] studied the effect of SMA fiber on the flutter of the composite plate with delaminations. Recently, Kuo et al. [6] presented a study of the influence of SMA on the flutter of buckled composite laminates by varying the SMA fiber spacing.

However, most of the previous studies were performed based on finite element modeling of structures [3-6]. It can be seen that few researches presented analytical solutions [1-2] for flutter analysis of the composite panels with SMA. Furthermore, it is found that the thermo-mechanical behavior was described by the expression of approximate fitting or use of the approximate data obtained from the curves of SMA [3-4, 6]. And some important factors affecting the flutter behaviour of the composite panels with SMA have not studied enough and the time domain response was not provided.

In this work, an analytical solution is presented to investigate the effects of SMA fibers on the flutter behavior of composite panels. The thermo-mechanical constitutive equation of SMA proposed by Brinson [7] and the mixture rule for evaluating the elastic properties of a SMA hybrid composite layer [8] are employed to derive the constitutive equation of the SMA reinforced composite panels. The composite panels model uses laminated plate theory. The Galerkin approach is carried out to reduce the motion equations of composite panels with SMA fibers to a system of ordinary differential equations. Based on these approximate equations the flutter behavior of the composite panels are then calculated in the frequency and time domains. Finally, the influence of parameters including SMA activation on flutter behavior of the composite panels is examined.

2. Formulations

As shown in Fig. 1, a composite panel with length a, width b and thickness h is considered. The SMA fibers are embedded parallel to the fibers of composite medium in the arbitrary layers. Supersonic airflow with air density ρa, flow velocity Ua, Mach number M is imposed on the composite panel. We assume that aerodynamic pressure P is passing over the top surface of the composite panel along the x-axis.

Fig. 1Geometry of a composite panel embedded with SMA fibers.

Geometry of a composite panel embedded with SMA fibers.

The aerodynamic pressure obtained from quasi-steady supersonic aerodynamic theory has the form:

1
P=-ρaUa2M2-1wx+1UaM2-2M2-1wt,

where w is the transverse displacement component in the mid-plane of the composite panel.

The effect of the w/t term, which causes aerodynamic damping, on the flutter boundary is insignificant [9] and is not included in this work.

By using the laminated plate theory and the one-dimensional model of SMA proposed by Brinson [7], the equation of motion for the composite panels with SMA fibers can be written as:

2
m02wt2+D114wx4+2D12+2D664wx2y2+D224wy4+4D164wx3y
+4D264wxy3 +Nx2wx2+Ny2wy2+2Nxy2wxy=p,

where m0 is the panel mass per unit area, the panel bending stiffnesses Dij are defined as:

3
Dij=-h/2h/2Q-ijz2dz, i,j=1,2,6,

where Q-ij is the transformed stiffness matrix of the layer, and:

4
Q¯=TQ(T)T,

in which:

5
T=n2m2-2mnm2n22mnmn-mnn2-m2,

where m and n are the cosine and sine functions of the ply angle in each lamina, respectively, and the stiffness matrix Q has following form:

6
Q=Q11Q120Q21Q22000Q66,

where:

7
Q11=E11-ν12ν21,Q12=Q21=ν12E21-ν12ν21,Q22=E21-ν12ν21,Q66=G12.

The modulus and Poisson’s ratio of the composite panel can be determined by considering the mixture rule of composition [8]:

8
E1=E1m1-Vs+EsVs,E2=E2mEsE2mVs+Es(1-Vs),G12=G12mGsG12mVs+Gs(1-Vs),ν12=ν12m1-Vs+νsVs,ν21=E2ν12E1,

where, the subscript m and s denote the composite matrix and SMA fibers, respectively, V is volume fraction.

Nx, Ny and Nxy in Eq. (2) are the resultant in-plane forces induced by SMA fibers and temperature, which can be given as:

9
(Nx,Ny,Nxy)=-h/2h/2σxr-σxT,σyr-σyT,σxyr-σxyTdz,

where the superscripts ‘r’ and ‘T’denote the in-plane stresses induced by recovery stresses of SMA fibers and temperature, respectively.

The recovery stresses of SMA fibers and the stresses induced by temperature of the k layer after coordinate transformation have the following form, respectively:

10
σxrσyrσxyr=m2n2-2mnn2m22mnmn-mnm2-n2σrVs00,
11
σxTσyTσxyT=Q-kcαxαyαxyT,

where, T denotes temperature increment, Q-kc denotes the elastic constant of composite medium, αx, αy and αxy are thermal expansion coefficients of the SMA hybrid panel which can be obtained as:

12
αxαyαxy=m2n2-mnn2m2mn2mn-2mnm2-n2α1α20.

Based on the one-dimensional model of SMA proposed by Brinson [7] and assuming all SMA fibers are fully constrained, the expressions for the recovery stress of SMA fibers during heating and cooling can be derived as follows, respectively [10].

During heating:

13
σr=θT-T0+σ0,0TAsσ,Eξε0-Eξ0ε0+Ωξξs-Ωξ0ξs0+θT-Asσ+σ1,AsσTAfσ,θT-Afσ+σ2,T Afσ.

During cooling:

14
σr=θT-Afσ+σ2,MsσTAfσ,Eξε0-Eξ0ε0+Ωξξs+θT-Msσ+σ3,MfσTMsσ,

where, ξ denotes the martensite fraction, E(ξ) denotes the elastic modulus of SMA, θ denotes the thermal elastic modulus, T denotes temperature and T0 is the reference temperature, Ω(ξ) is the phase transformation coefficient. The subscript 0 denotes initial state. Asσ and Afσ denote the start and finish temperatures of austenite in stress, and Msσ and Mfσ denote the start and finish temperatures of martensite in stress. ξs is the martensite fraction induced by stress. The first expression in Eqs. (13) is used for SMA in the initial martensite state, while the third expression among them is used for SMA in 100 % austenite state, and the second expression among them is used for SMA in the phase transformation state from martensite to austenite.

The transformation process from martensite to austenite (M→A) and from austenite to martensite (A→M) are given as follows, respectively:

15
ξMA=ξs02cosπAf-AsT-As-σrCA+1,
16
ξAM=1-ξs02cosπσscr-σfcrσr-σfcr-CM(T-Ms)+1+ξs02,

where, Msand Mf denote the start and finish temperature points during the phase transformation from austenite to martensite, and Asand Af denote the start and finish temperature points during the phase transformation from martensite to austenite, respectively. CA and CM are material constant of SMA fibers which determine the influence of the stress on the transformation temperature. σscr and σfcr denote the start and finish stress points of phase transformation, respectively, and the subscript “s0” denotes martensite fraction induced by stress corresponding to initial state.

In Eqs. (13) and (14), σ1 and σ2 denote the start and finish stress for the transformation (M→A), σ3 denotes the start stress for the transformation (A→M).

The phase transformation coefficient Ωξ and the elastic modulus Es(ξ) can be expressed respectively as:

17
Eξ=EA+ξEM-EA,
18
Ωξ=-εLEξ,

where, EM and EA denote Young’s modulus in the martensite and austenite phase, respectively, εL is recovery strain limit.

3. Solution of the governing equation

The composite panel is simply supported and the mode shapes that satisfy the displacement boundary conditions can be assumed as:

19
wx,y,t=i=1NΦix,yqit=i=1Nsinmiπxasinniπybqi(t).

Substituting Eq. (19) into the governing Eq. (1) and applying Galerkin procedure, the following governing equations in matrix form can be found:

20
Mq¨(t)+K+λBqt=0,

where q(t)=q1t,,qNt is generalized coordinate vector, and dot denotes the derivative with respect to time t. M, K and B are the mass and stiffness and aerodynamic influence matrices, respectively. They can be expressed as:

21
Mij=m0a44D110,i=j,Mij=0,ij,
22
Kii=π44D11D110mi4+D22D110γ4ni4+2D12D110+4D66D110γ2mi2ni2+π24N¯xπ2mi2+N¯yπ2γ2ni2,i=j,Kij=0,ij,
23
Bij=0,i=j,Bij=mi41mi-mj+1mi+mj-cosπmi-mjmi-mj-cosπ(mi+mj)mi+mj, ij,

where γ denotes the length to width ratio, D110 denotes the bending stiffness of the composite panel when all the fibers of the composite layers are aligned in the airflow x direction and Vs equal to zero.

Nondimensional aerodynamic pressure is given as:

24
λ=ρaUa2(M2-1)1/2/D110a4.

Nondimensional in-plane forces are defined as:

25
N¯x=β3NxD-110, N¯y=β3NyD-110,

where β is the length to thickness ratio.

4. Numerical results

4.1. The critical flutter dynamic pressure

Flutter behavior of the composite panels with embedded SMA fibers is studied. The composite matrix material used in the numerical examples is Graphite/Epoxy AS-3002. Its properties are listed in Table 1. The [0/90/–90/90]s and [0/45/–45/90]s laminates are used to examine the effect of the lay-up on flutter behavior. The material properties of SMA fibers are listed in Table 2.

Table 1The material parameters of Graphite/Epoxy AS-3002

E1m (GPa)
26.5 E2m
G12m (GPa)
1.184 E2m
v12m
0.21
α1m (1/°C)
–0.07×10-6
α2m (1/°C)
30.1×10-6

Table 2Material properties of SMA wires

SMA (Nitinol)
EA / MPa
67×103
EM / MPa
26.3×103
σscr / MPa
100
σfcr / MPa
170
εL
0.067
Θ / (MPa/°C)
0.55
Mf / °C
9
Ms / °C
18.4
As / °C
34.5
Af / °C
49
CM / (MPa/°C)
8
CA / (MPa/°C)
13.8
ρ / (kg/m3)
6450

The thickness of the composite panel is 0.002 (m), β= 200. The width can be obtained by b=βh/γ so long as γ is given.

Fig. 2Frequency coalescence curves of a SMA composite panel with increase of number of vibration modes N(ε0= 0.01, Vs= 0.4, T= 55°C, γ= 2)

Frequency coalescence curves of a SMA composite panel with increase of number of vibration modes N(ε0= 0.01, Vs= 0.4, T= 55°C, γ= 2)

Fig. 3Flutter boundaries of a SMA composite panel with increase of number of vibration modes N(ε0= 0.01, Vs= 0.4, γ= 2)

Flutter boundaries of a SMA composite  panel with increase of number of vibration  modes N(ε0= 0.01, Vs= 0.4, γ= 2)

It is necessary to examine if sufficient mode shape functions are used in the solution of the equation to obtain converged results in the case of using Galerkin method. Fig. 2 shows the frequency coalescence curves of the SMA composite panel with [0/90/–90/90]s laminates for increase of number of vibration modes N. When λ= 0, the eigenvalues of all modes of the panel are real. λ may be increased from zero until the eigenvalues of two modes coalesce or become identical. This value of λ that creates coalescence is the critical value λcr, where flutter occurs.

It is seen from the figure that, fourteen modes of the composite panel are enough to obtain converged solution.

Fig. 3 shows the critical flutter dynamic pressure vs. temperature curves for the SMA composite panel with [0/90/–90/90]s laminates for increase of number of vibration modes N. The result shows that the approximate solution methodology developed in this paper can obtain good convergence solutions.

Based on mode convergence examination it is found that N= 14 gives suitably converged results. So for all results given in this paper N= 14 unless otherwise noted.

4.1.1. Symmetrical orthotropic panels

In this section, the symmetrical orthotropic panels with SMA fibers are considered. The stacking sequence is [0/90/–90/90]s.

Fig. 4 presents the flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of Vs(ε0= 0.01, γ= 2). As shown in Fig. 4, The critical flutter dynamic pressure of the composite panel without SMA fibers (Vs= 0) is seen decrease linearly with increase of temperature. This is because there is not the recovery stress of the SMA fiber in this case, as a result of the stiffness softened, therefore, λcr is reduced. However, it is seen from Fig. 4 that SMA fiber embedding result in significant increase of λcr compared with the result of the composite panels without SMA. The phase transformation M→A is induced during heating, and the phase transformation A→M is induced during cooling. The curves of temperature-critical flutter dynamic pressure (λcr-T) are typical hysteresis loops in a thermal cycle. λcr increases with the increase in temperature during heating, but decline with the decrease in temperature during cooling. Furthermore, it can be seen that any increase of the volume fraction of SMA fibers Vs, leads to an obvious increase of λcr.

Fig. 5 presents the effect of the length to width ratio on the critical flutter dynamic pressure of the composite panels. And the increase of γ is accompanied by the shift of these hysteresis loops towards higher critical flutter dynamic pressure. It seems that increasing γ results in a more dynamically stable panels.

Fig. 6 presents the influence of the initial strain of the SMA fibers on the critical flutter dynamic pressure of the composite panels. As can be seen, the initial strain of SMA fiber appears to have marginal effect on λcr.

Fig. 4Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of Vs for [0/90/–90/90]s laminates (ε0= 0.01, γ= 2)

Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of Vs for  [0/90/–90/90]s laminates (ε0= 0.01, γ= 2)

Fig. 5Flutter boundaries of SMA/Graphite/Epoxy AS-3002 with different values of γ for [0/90/–90/90]s laminates (ε0= 0.01, Vs= 0.4)

Flutter boundaries of SMA/Graphite/Epoxy AS-3002 with different values of γ for  [0/90/–90/90]s laminates (ε0= 0.01, Vs= 0.4)

Fig. 6Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of ε0 for [0/90/–90/90]s laminates (Vs= 0.4, γ= 2)

Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of ε0 for [0/90/–90/90]s laminates (Vs= 0.4, γ= 2)

Fig. 7Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of Vs for [0/45/–45/90]s laminates (ε0= 0.01, γ= 2)

Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of Vs for  [0/45/–45/90]s laminates (ε0= 0.01, γ= 2)

4.1.2. Symmetrical anisotropic panels

In this section, the symmetrical anisotropic panels with SMA fibers are considered. The stacking sequence is [0/45/–45/90]s. Fig. 7 presents the flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels. It is plotted with different values of Vs. Like the symmetrical orthotropic panels, the critical flutter dynamic pressure increases as the volume fraction of SMA fibers increases.

Figs. 8 and 9 present the flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels. They are plotted with different values of γ and ε0, respectively. The results are similar to the symmetrical orthotropic case. However, it can be evidently seen by comparing results in Figs. 4-6 and Figs. 7-9 that the stacking sequence of the composite panels with embedded SMA fibers significantly affects the critical flutter dynamic pressure.

Fig. 8Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of γ for [0/45/–45/90]s laminates (ε0= 0.01, Vs= 0.4)

Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of γ for  [0/45/–45/90]s laminates (ε0= 0.01, Vs= 0.4)

Fig. 9Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of ε0 for [0/45/–45/90]s laminates (Vs= 0.4, γ= 2)

Flutter boundaries of SMA/Graphite/Epoxy AS-3002 panels with different values of ε0 for [0/45/–45/90]s laminates (Vs= 0.4, γ= 2)

4.2. Dynamical response

Fig. 10 presents the time history responses and phase plane plots of the geometric center of the panel for various dynamic pressures (λ<λcr, λ>λcr, λ=λcr).This figure is plotted for the SMA fibers in the initial martensite state during heating (T= 30 °C). It clearly shows that when dynamic pressure is larger than the critical flutter dynamic pressure, an unstable response occurs.

Fig. 11 presents the time history response and phase plane plots of the geometric center of the panel for various dynamic pressures (λ<λcr, λ>λcr, λ=λcr). This figure is plotted for the SMA fibers in the phase transformation state from austenite to martensite during cooling (T= 60 °C). The results are similar to the case of symmetrical orthotropic in Fig. 9. It also shows that when dynamic pressure is smaller than the critical flutter dynamic pressure, the response of the composite panel is stable.

Fig. 10Displacement response curves with time and phase plane plots during heating ([0/90/–90/90]s, T= 30 °C)

Displacement response curves with time  and phase plane plots during heating ([0/90/–90/90]s, T= 30 °C)

a)λ<λcr

Displacement response curves with time  and phase plane plots during heating ([0/90/–90/90]s, T= 30 °C)

b)λ=λcr

Displacement response curves with time  and phase plane plots during heating ([0/90/–90/90]s, T= 30 °C)

c)λ>λcr

5. Conclusions

An analytical model was presented for the study of the flutter behavior of the composite panels embedded with SMA fibers. The presented model was used to predict the critical flutter dynamic pressure and the aeroelastic stability. Emphasis is placed on the study of the effects of SMA activation. From the present analysis and the numerical results, the following main conclusions can be obtained.

Fig. 11Displacement response curves with time and phase plane plots during cooling ([0/45/-45/90]s, T= 60 °C)

Displacement response curves with time  and phase plane plots during cooling ([0/45/-45/90]s, T= 60 °C)

a)λ<λcr

Displacement response curves with time  and phase plane plots during cooling ([0/45/-45/90]s, T= 60 °C)

b)λ=λcr

Displacement response curves with time  and phase plane plots during cooling ([0/45/-45/90]s, T= 60 °C)

c)λ>λcr

1) The effects of temperature on flutter behavior during the phase transformation are important. The curves of the critical flutter dynamic pressure vs. temperature are typical hysteresis loops in a thermal cycle. The reason underlying this phenomenon is the well-known recovery stress-temperature hysteresis behavior of SMA fibers itself.

2) The recovery stress resulted due to SMA fiber actuation significantly increases the critical flutter dynamic pressure. The effects of the recovery stress on the critical flutter dynamic pressure can be enhanced by increasing of both the initial strain and the volume fraction of SMA fibers.

3) The lay-up affects flutter behavior of the composite panels significantly, whereas the initial strain of SMA fiber appears to have marginal effect on flutter behavior of the composite panels.

References

  • Weisshaar Sadlowski T. A. M. J. The Effect of SMA Micro-Actuation on Panel Flutter. AAE Report 92-100. School of Aeronautics and Astronautics, Purdue University, West Lafayette, 1992.
  • Scott R. C., Weisshaar T. A. Panel flutter suppression using adaptive material actuators. Journal of Aircraft, Vol. 31, Issue 1, 1994, p. 213-222.
  • Park J. S., Kim J. H., Moon S. H. Thermal post-buckling and flutter characteristics of composite plates embedded with shape memory alloy fibers. Composites: Part B, Vol. 36, 2005, p. 627-636.
  • Duan B., Abdel-Motagaly K., Guo X. Y., et al. Flutter and thermal deflection suppression of composite plates using shape memory alloy. AIAA Journal, Vol. 43, Issue 9, 2005, p. 2015-2023.
  • Ostachowicz W. M., Kaczmarczyk S. Vibrations of composite plates with SMA fibers in a gas stream with defects of the type delamination. Composite Structures, Vol. 54, 2001, p. 305-311.
  • Kuo S. Y., Shiau L. C., Lai C. H. Flutter of buckled shape memory alloy reinforced laminates. Smart Materials and Structures, Vol. 21, 2012, p. 1-11.
  • Brinson L. C. One-dimensional constitutive behavior of shape memory alloys: thermmechanical derivation with non-constant material functions and redefined variable. Journal of Intelligent Material Systems and Structures, Vol. 4, 1993, p. 229-242.
  • Zhong Z. W., Chen R. R., Mei C., Pates III C. S. Buckling and postbuckling of shape memory alloy fiber-reinforced composite plates. Proceedings of the ASME Aerospace Division, Mech Eng Congress and Exp, Chicago, 1994, p. 115-132.
  • Ashley H., Zartarian G. Piston theory-a new aero-dynamic tool for the aeroelastician. Journal of the Aeronautical Sciences, Vol. 23, Issue 12, 1956, p. 1109-1118.
  • Ren Y. S., Dai Q. Y., An R. J., Zhu Y. F. Modeling and dynamical behavior of rotating composite shafts with SMA wires. Shock and Vibration, 2014, p. 765875.

Cited by

Investigation of dynamic behavior of Hollow Rotor Shaft System through FEM
Ashish Gupta | Viswajeet Singh Yadav | Sandeep Chauhan
(2021)

About this article

Received
23 December 2015
Accepted
24 May 2016
Published
30 June 2016
SUBJECTS
Flow induced structural vibrations
Keywords
shape memory alloys
composite panel
flutter
Acknowledgements

The research is funded by the National Natural Science Foundation of China (Grant No. 11272190) and Graduate Innovation Project of Shandong University of Science and Technology of China (Grant No. YC150325).