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

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.


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][4][5][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.

Formulations
As shown in Fig. 1, a composite panel with length , width and thickness ℎ is considered.The SMA fibers are embedded parallel to the fibers of composite medium in the arbitrary layers.Supersonic airflow with air density , flow velocity , Mach number is imposed on the composite panel.We assume that aerodynamic pressure ∆ is passing over the top surface of the composite panel along the -axis.The aerodynamic pressure obtained from quasi-steady supersonic aerodynamic theory has the form: where is the transverse displacement component in the mid-plane of the composite panel.The effect of the ⁄ 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: where is the panel mass per unit area, the panel bending stiffnesses are defined as: where is the transformed stiffness matrix of the layer, and: in which: where and are the cosine and sine functions of the ply angle in each lamina, respectively, and the stiffness matrix has following form: where: The modulus and Poisson's ratio of the composite panel can be determined by considering the mixture rule of composition [8]: where, the subscript m and s denote the composite matrix and SMA fibers, respectively, is volume fraction., and in Eq. ( 2) are the resultant in-plane forces induced by SMA fibers and temperature, which can be given as: where the superscripts ' ' and '∆ '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 layer after coordinate transformation have the following form, respectively: where, ∆ denotes temperature increment, denotes the elastic constant of composite medium, , and are thermal expansion coefficients of the SMA hybrid panel which can be obtained as: 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: During cooling: where, denotes the martensite fraction, ( ) denotes the elastic modulus of SMA, denotes the thermal elastic modulus, denotes temperature and is the reference temperature, Ω( ) is the phase transformation coefficient.The subscript 0 denotes initial state.and denote the start and finish temperatures of austenite in stress, and and denote the start and finish temperatures of martensite in stress.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: where, and denote the start and finish temperature points during the phase transformation from austenite to martensite, and and denote the start and finish temperature points during the phase transformation from martensite to austenite, respectively.and are material constant of SMA fibers which determine the influence of the stress on the transformation YONGSHENG REN, CHENGGANG DU, XIAOQING MA temperature.and denote the start and finish stress points of phase transformation, respectively, and the subscript " " denotes martensite fraction induced by stress corresponding to initial state.
In Eqs. ( 13) and ( 14), and denote the start and finish stress for the transformation (M→A), denotes the start stress for the transformation (A→M).
The phase transformation coefficient Ω( ) and the elastic modulus ( ) can be expressed respectively as: where, and denote Young's modulus in the martensite and austenite phase, respectively, is recovery strain limit.

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: Substituting Eq. ( 19) into the governing Eq. ( 1) and applying Galerkin procedure, the following governing equations in matrix form can be found: where ( ) = ( ), … , ( ) is generalized coordinate vector, and dot denotes the derivative with respect to time .
, and are the mass and stiffness and aerodynamic influence matrices, respectively.They can be expressed as: where denotes the length to width ratio, denotes the bending stiffness of the composite panel when all the fibers of the composite layers are aligned in the airflow direction and equal to zero.
Nondimensional aerodynamic pressure is given as: Nondimensional in-plane forces are defined as: where is the length to thickness ratio.

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.The thickness of the composite panel is 0.002 (m), = 200.The width can be obtained by = ℎ/ so long as is given.YONGSHENG REN, CHENGGANG DU, XIAOQING MA 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 , 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 .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 = 14 gives suitably converged results.So for all results given in this paper = 14 unless otherwise noted.

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 ( = 0.01, = 2).As shown in Fig. 4, The critical flutter dynamic pressure of the composite panel without SMA fibers ( = 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, is reduced.However, it is seen from Fig. 4 that SMA fiber embedding result in significant increase of 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 ( − ) are typical hysteresis loops in a thermal cycle.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 , leads to an obvious increase of .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 .).This figure is plotted for the SMA fibers in the phase transformation state from austenite to martensite during cooling ( = 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.

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.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

Fig. 1 .
Fig. 1.Geometry of a composite panel embedded with SMA fibers.
2060.AN ANALYTICAL MODEL FOR FLUTTER BEHAVIOR OF COMPOSITE PANELS WITH SHAPE MEMORY ALLOY FIBERS.

Fig. 2 .Fig. 3 .
Fig. 2. Frequency coalescence curves of a SMA composite panel with increase of number of vibration modes ( = 0.01, = 0.4, = 55°C, = 2) Fig. 3. Flutter boundaries of a SMA composite panel with increase of number of vibration modes ( = 0.01, = 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.2shows the frequency coalescence curves of the SMA composite panel with [0/90/-90/90] s laminates for increase of number of vibration modes .When = 0, the eigenvalues of all modes of the panel

Fig. 11 .
Fig. 11.Displacement response curves with time and phase plane plots during cooling ([0/45/-45/90]s, = 60 °C)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 2060.AN ANALYTICAL MODEL FOR FLUTTER BEHAVIOR OF COMPOSITE PANELS WITH SHAPE MEMORY ALLOY FIBERS.YONGSHENG REN, CHENGGANG DU, XIAOQING MA

Table 1 .
The