The flutter and thermal buckling behavior of laminated composite panels embedded with shape memory alloy (SMA) wires are studied in this research. The classical plate theory and nonlinear von-Karman strain-displacement relation are employed to investigate the aeroelastic behavior of the smart laminated panel. The thermodynamic behaviors of SMA wires are simulated based on one-dimensional Brinson SMA model. The aerodynamic pressure on the panel is described by the nonlinear piston theory. Nonlinear governing partial differential equations of motion are derived for the panel via the Hamilton principle. The effects of ply angle of the composite panel, SMA layer location and orientation, SMA wires temperature, volume fraction and prestrain on the buckling, flutter boundary, and amplitude of limit cycle oscillation of the panel are analyzed in detail.
National Natural Science Foundation of China912161061. Introduction
Thin panel is a common and useful form of structure component, which has been applied significantly in high-speed vehicles, such as aircraft, spacecrafts, and rockets. Panel flutter, which occurs over a critical velocity under the coupling actions of elastic, inertia, and aerodynamic force induced by the transonic, supersonic, or hypersonic airflow, is one kind of self-excited vibration. Flutter raises the amplitude of vibration as well as the aerodynamic stress dramatically, which results in the failure of the structure. The flutter phenomena were observed on the V-2 rockets for the first time during World War II; since then, lots of studies have been carried out using different structural and aerodynamic theories. Excellent surveys of early researches were presented by Dowell [1, 2]. Composite materials with the advantages of high strength, light weight, and low coefficient of thermal expansion have been broadly employed in the design of thin-panel structures. Mei et al. [3] presented a recent survey about analytical methods of nonlinear panel flutter in supersonic air flow. Birman and Librescu [4] analyzed aeroelastic instability of laminated composite panels for which the shear deformation is considered in the modeling process under supersonic flow. Kouchakzadeh et al. [5] studied the nonlinear aeroelasticity problem of a laminated composite plate under supersonic airflow, where the classic plate theory was adopted to establish the structure dynamic model, and the supersonic airflow was simulated via linear piston theory. Results obtained in [5] indicated that the ply angle had important influence on the flutter behavior. Zhao and Cao [6] considered the aerodynamic nonlinearity in the modeling process of a stiffened laminate composite panel under supersonic flow. Numerical results in [6] indicated that the position, thickness, and width of the stiffener had significant effect on the aeroelastic behavior.
High-speed aircrafts are subjected to aerodynamic pressure and aerodynamic heating, which should be taken into account when solving the aeroelastic problems. Abbas et al. [7] analyzed the aerothermoelastic behavior of the isotropic and orthotropic panel. Shiau et al. [8] analyzed the influence of temperature gradient on flutter phenomena of a composite laminated plate based on the finite element method (FEM). Xie et al. [9] employed the proper orthogonal decomposition (POD) method to analyze the flutter behavior of the panel with uniform thermal loadings. Xue and Mei [10] analyzed aeroelastic problem of the plate under nonuniform thermal loadings by the finite element frequency domain method. Li and Song [11] employed the assumed mode approach and FEM to investigate aerothermoelastic behavior of a composite panel.
In the last two decades, numerous works have been involved in suppressing panel flutter by active or passive control utilizing smart materials. Among these materials, SMA, which is most suitable for active control of composite structures, has been extensively studied. Birman [12] presented a review about various applications of SMA in industry. By heating up to the austenite finish temperature, SMA is able to recover large prestrain totally. SMA can generate large recovery stresses, when the prestrain is restrained. Through embedding SMA in the form of wires within the laminated composites, the recovery stresses can modify stiffness of the structures. This could improve structure characteristics of composite panel such as vibration, thermal bulking, impact loading, flutter, and acoustic. Park et al. [13] employed FEM to study the influences of SMA fibers on flutter and buckling behavior of the plate. Ostachowicz et al. [14] used FEM to analyze the buckling and flutter behaviors of composite plates and the recovery stress generated by SMA fibers embedded in the plate is obtained from experimental data. In accordance with the first-order shear deformation plate theory, Barzegari et al. [15] studied the aeroelastic behavior of rectangular cantilever isotropic wings embedded with SMA wires, where the aerodynamic loading is estimated by linear piston theory. Asadi et al. [16] studied the problem of vibration and thermal buckling for a laminated composite beam in which the SMA fibers are embedded symmetrically and asymmetrically. Kuo et al. [17] used FEM to research flutter behavior of buckled SMA reinforced laminates. The nonlinear aeroelastic behavior of a SMA hybrid composite plate was studied by Ibrahim et al. [18] utilizing a novel FEM.
However, studies mentioned above are focused on thermal buckling and flutter of composite panels embedded with SMA wires employing the FEM without giving consideration to the aerodynamic nonlinearity in the modeling process. Moreover, the well-developed one-dimensional Brinson model for the thermodynamic behaviors of SMA is seldom utilized in the flutter analysis. In this study, in order to analyze the smart laminated panel’s dynamic characteristics, the classical plate theory is employed to derive the nonlinear governing differential equations of motion, in which the nonlinear von-Karman strain-displacement relation is adopted. In the modeling process, the thermodynamic behaviors of SMA wires are simulated based on one-dimensional Brinson SMA model, while the aerodynamic pressure is calculated by nonlinear piston theory. The Galerkin method is adopted to derive the system discrete dynamic model, which can be solved numerically by the Runge-Kutta method. The numerical results are utilized to show the effect of SMA wires on the natural frequency, buckling, flutter, and amplitude of limit cycle oscillation (LCO) of a laminated composite panel.
2. Aeroelastic Model
Consider an 8-layer symmetric composite panel embedded with SMA wires in Cartesian coordinate system with thickness h, length a, and width b as displayed in Figure 1. SMA wires are aligned in fiber direction in arbitrary layer. The supersonic flow is along the positive x direction and SMA wires are embedded in the second and seventh layer symmetrically as shown in Figure 1.
Schematic illustration of a smart laminated panel.
2.1. Structure Modeling
The classical plate theory is employed to describe the displacements in the panel; that is,(1)ux,y,z,t=u0x,y,t-z∂w0x,y,t∂x,vx,y,z,t=v0x,y,t-z∂w0x,y,t∂y,wx,y,z,t=w0x,y,t,where u, v, and w denote the displacements in the x, y, and z directions, respectively. Subscript “0” stands for the midplane displacement. The nonlinear strain-displacement relations, according to the von-Karman assumption, are given by(2)εxεyγxy=∂u0∂x+12∂w0∂x2-z∂2w0∂x2∂v0∂y+12∂w0∂y2-z∂2w0∂y2∂u0∂y+∂v0∂x+∂w0∂x∂w0∂y-2z∂2w0∂x∂y.
Based on Hooke’s law, the constitutive equation of the nth layer of the laminated composite panel under thermal loads is (3)σxσyτxyn=T-1QmnT-Tεx-αxΔTnεy-αyΔTnγxy-αxyΔTn=Q-mnεx-αxΔTnεy-αyΔTnγxy-αxyΔTn,where the subscripts “n” and “m” indicate the layer number and the composite matrix, ΔT denotes the temperature variation, and α stands for the thermal expansion coefficient. The transform matrix T and stiffness matrix Qm are defined as [23](4)T=cos2θsin2θ2sinθcosθsin2θcos2θ-2sinθcosθ-sinθcosθsinθcosθcos2θ-sin2θ,(5)Qm=Q11mQ12m0Q12mQ22m000Q66m,respectively. Here, Q11m, Q22m, Q12m, and Q66m are the stiffness coefficients which are defined as (6)Q11m=E1m1-υ12mυ21m,Q22m=E2m1-υ12mυ21m,Q12m=υ12mE2m1-υ12mυ21m,Q66m=G12m,in which G12m, E1m, and E2m stand for the shear modulus and Young modulus and υ12m and υ21m are Poisson ratios.
And for the kth layer of the laminate composite panel embedded with SMA wires, the constitutive equation is as follows [18]:(7)σk=Q-kε+T-1σrkVsk-Q-mkαΔTkVmk,where subscripts “s” and “k” indicate SMA wires and layer number. σ is in-plane stress vectors and σr, to be described in Section 2.2, is the SMA recovery stress vectors under the temperature T. Q- is transformed reduced stiffness matrix of the smart layer. V stands for the volume fractions. The elastic properties used in Q- are expressed as(8)E1=EsVs+E1mVm,υ12=υsVs+υ12mVm,E2=EsE2mEsVm+E2mVs,G12=GsG12mGsVm+G12mVs.
2.2. Description of the Stress Model of SMA Wires
According to one-dimensional model of SMA proposed by Brinson [24] and assuming that all SMA wires are fully constrained, one can get the following expressions for the recovery stress of SMA wires:(9)σr=σ0+ΘT-T00≤T≤Asσσ1+Esξ-Esξ0ε0+Ωξξs-Ωξ0ξs0+ΘT-AsσAsσ≤T≤Afσσ2+ΘT-AfσT≥Afσ.
The two constants in (9) are(10)σ1=σ0+ΘAsσ-T0,σ2=σ1+EA-Esξ0ε0+Ωξ0ξs0+ΘAfσ-Asσ.
The phase transformation coefficient Ω(ξ) and the elastic modulus Es(ξ) have the following expressions:(11)Ωξ=-εLEsξ,Esξ=EA+ξEM-EA.
The austenite start temperature Asσ and finish temperature Afσ under stress can be expressed, respectively, as(12)Asσ=aAAs+bAΘT0-σ0aA+bAΘ,Afσ=aAAs+π-bAEA-Esξ0ε0-Ωξ0ξs0+σ0-ΘT0aA+bAΘ,where ξ denotes the martensite fraction, Es(ξ) denotes the elastic modulus of SMA, Θ represents the thermal elastic modulus, T represents temperature of SMA wires and T0 denotes the reference temperature, subscript “0” denotes initial conditions, and ξs stands for the martensite fraction induced by stress. The first expression in (9) is used for SMA in the initial martensite state, the second one is used for SMA in the phase transformation state, and the third one is used for SMA in 100% austenite state.
The relationships of SMA characteristics transformed from martensite to austenite are given as(13)ξ=ξ02cosaAT-As-σrcA+1,ξs=ξs0-ξs0ξ0ξ0-ξ,ξT=ξT0-ξT0ξ0ξ0-ξ,ξ0=ξs0+ξT0,aA=πAf-As,bA=-aAcA,aM=πMs-Mf,bM=-aMcM,where ξT denotes martensite fraction induced by temperature and cM and cA are phase transformation constants. Figure 2 shows the computed SMA recovery stress versus various temperatures with four prestrain levels. It is shown that, in the phase transformation state, SMA wires can generate large recovery stresses. Moreover, for higher prestrain, an increase in the temperature will generate larger recovery stresses.
Recovery stress generated by SMA.
2.3. Aerodynamic Pressure Modeling
The aerodynamic load, Δp, is described by the third-order piston theory as(14)Δp=2qMw˙0v∞+∂w0∂x+γ+1M4w˙0v∞+∂w0∂x2+γ+1M212w˙0v∞+∂w0∂x3,where, ρa, v∞, γ, and M denote the density, velocity, ratio of specific heats, and Mach number of airflow and q=ρav∞2/2.
2.4. Governing Equations of Motion
The partial differential equations of motion for the panel can be obtained by utilizing Hamilton principle.(15)∫t0t1δT+δW-δUdt=0,where the variations of kinetic energy δT and virtual work δW as well as variation of strains energy δU are given by(16)δT=∫0b∫0a∫-h/2h/2ρu˙δu˙+v˙δv˙+w˙δw˙dzdxdy,δW=-∫0b∫0aΔpδwx,y,0.5hdxdy,δU=∫0b∫0a∫-h/2h/2σxδεx+σyδεy+τxyδγxydzdxdy,where ρ stands for density of the smart laminated panel.
Substituting (1)–(5), (7), (8), and (16) into (15) and setting the coefficients of δu, δv, and δw to be zero, one has(17)∂Nxx∂x+∂Nxy∂y=I0u¨0,∂Nyy∂y+∂Nxy∂x=I0v¨0,∂Nxx∂x+∂Nxy∂y∂w0∂x+∂Nyy∂y+∂Nxy∂x∂w0∂y+Nxx∂2w0∂x2+Nyy∂2w0∂y2+2Nxy∂2w0∂x∂y+∂2Mxx∂x2+∂2Myy∂x2+2∂2Mxy∂x∂y-Δp=I0w¨0-I2∂2w¨0∂x2+∂2w¨0∂y2,where (I0,I2)=∫-h/2h/2ρ(1,z2)dz and the force resultants operators are obtained as(18)NxxNyyNxy=∫-h/2h/2σxσyτxydz,MxxMyyMxy=∫-h/2h/2σxσyτxyzdz.
The boundary conditions for a simply supported panel are(19)u0x=0,a=0,u0y=0,b=0,v0x=0,a=0,v0y=0,b=0,w0x=0,a=0,w0y=0,b=0,∂w02∂x2x=0,a=0,∂w02∂y2y=0,b=0.
The displacements of the panel satisfying the boundary conditions are written as follows:(20)u0x,y,t=∑i=1l∑j=1maijtΦijx,y,v0x,y,t=∑i=1l∑j=1mbijtΦijx,y,w0x,y,t=∑i=1l∑j=1mcijtΦijx,y,where the mode shapes are taken as(21)Φijx,y=siniπxasinjπyb.
The numerical results provided by Dowell [25] show that for getting reasonable results at least four modes are needed to study the panel flutter behavior. In present study, four streamwise modes and one spanwise mode are reserved in the following calculation. Substituting (20) into (17) and after that integrating over the panel area, the discrete dynamic model is given as(22)a¨i=4abI0∫0a∫0b∂Nxx∂x+∂Nxy∂yΦix,ydxdy,b¨i=4abI0∫0a∫0b∂Nyy∂y+∂Nxy∂xΦix,ydxdy,c¨i=4aba2b2I0+i2b2+a2π2I2∫0a∫0b∂Nxx∂x+∂Nxy∂y∂w0∂x+∂Nyy∂y+∂Nxy∂x∂w0∂y+∂2Mxx∂x2+∂2Myy∂y2+2∂2Mxy∂x∂y+Nxx∂2w0∂x2+Nyy∂2w0∂y2+2Nxy∂2w0∂x∂y-ΔpΦix,ydxdy.Let(23)y=a1,a˙1,…,a4,a˙4,b1,b˙1,…,b4,b˙4,c1,c˙1,…,c4,c˙4T,and then (22) can be expressed as(24)y˙=Ay+gy,where A is the Jacobian matrix at the equilibrium point y=0 and g(y) respects the nonlinear terms induced by geometric and aerodynamic nonlinearity.
3. Results and Discussion
As the dynamic pressure λ reaches the critical flutter dynamic pressure λcr, the motion of the panel changes from flat condition to flutter condition based on the nonlinear theory. If λ>λcr, because of the existence of the geometric and aerodynamic nonlinearities, the amplitude of vibration of the panel will increase with time and eventually converge at a limit cycle. On the contrary, the amplitude of vibration will decrease with time if λ<λcr. The general solution of (24) can be written as (25)yt=y0eφjt,φj=αj±iβj,where y0 and φj stand for the eigenvector and eigenvalues of the matrix A. The natural frequencies of the panel can be obtained as (26)ωj=βj2.
As the real part of arbitrary eigenvalue turns from negative to positive, the flutter will happen. By examining the maximal real part of eigenvalues μ, λcr can be obtained. μ can be described as(27)μ=maxReφj=maxαj.
In the present study, the laminated panel with width of 0.4m, length of 0.5m, and thickness of 0.0025m is taken for analysis (except Section 3.3 where the effect of length-to-width ratios on λcr is discussed). The transverse displacements are plotted at the point (x,y)=(0.75a,0.5b). The material parameters are given as follows.
SMA Wires Properties(29)EA=67GPa,EM=26.3GPa,Ms=18.4C°,Mf=9C°,As=34.5C°,Af=49C°,CM=8MPa/C°,CA=13.8MPa/C°,Θ=0.55,εl=0.067,σ0=0,T0=20C°,ρs=6450kg/m3,ξT0=0,υs=0.33.
3.1. Validations of the Present Method
The numerical calculations are performed via MATLAB software. Based on the derived formulation, the Runge-Kutta method is utilized to investigate nonlinear flutter and buckling behavior of isotropic panel and laminated composite panel, respectively. The critical thermal buckling temperature is compared with the results provided by Matsunaga [22], Zhao et al. [19], Shiau et al. [20], and Shi et al. [21]. Tables 1 and 2 list the results of the former analyses and the present works. It can be observed from Tables 1 and 2 that the results obtained here have a good agreement with those results in the literature. In addition, the amplitude of LCO of the panel with thermal effect obtained here has a good agreement with those results in [26] as shown in Figure 3.
Critical buckling temperature of the isotropic panel.
Zhao et al. [19]
Shiau et al. [20]
Shi et al. [21]
Present
ΔTcr (°C)
126.5
127.1
126.4
126.43
Comparison of critical buckling temperature for composite panel.
Ply angle
ΔTcr (°C)
Matsunaga [22]
Present
[0/90/90/0]s
6.8
6.79
[0/45/−45/90]s
7.6
7.59
Comparison of limit cycle amplitudes.
3.2. The Effects of Orientation and Position of Layer Embedded with SMA Wires
The influences of orientation of SMA wires on λcr are investigated first. The fiber orientation of the composite panel is assumed to be [90/−45/45/θSMA]s with SMA wires embedded in the fourth and fifth layer symmetrically. The SMA wires have prestrain of 0.5%, volume fraction of 1%, and temperature of 50°C. The variation of λcr versus the angle of SMA wires is depicted in Figure 4. It can be observed from Figure 4 that changing the angle of layer embedded with SMA wires from 0° to 90° decreases λcr of the panel. Also, it is shown that the panel has the highest λcr with the ply angle of [90/−45/45/0SMA]s. Consequently, orientation of the layer embedded with SMA wires is the most important parameters for designing and optimizing the smart laminated panel. In the following analysis, the orientation of layer embedded with SMA wires is designed to be zero.
Flutter boundary versus orientation angles of SMA wires.
As for the influences of position of the layer embedded with SMA wires of the 8-layer symmetric laminated panel on the λcr change, it is found in Figure 5 that five cases of laminated panels with different position of the layer embedded with SMA wires have been studied. From Figure 5 one can see that when changing the positon of the layer embedded with SMA wires from outer layer [0SMA/−60/60/−60]s to inner layer [60/−60/60/0SMA]s in sequence, λcr may decrease, which indicates that embedding SMA wires in the outer layer is more significant for enhancing λcr of the laminated panel.
Flutter boundary versus position of layer embedded with SMA wires (i=0, [60/−60/60/−60]s;i=1, [0SMA/−60/60/−60]s;i=2, [60/0SMA/60/−60]s;i=3, [60/−60/0SMA/−60]s; i=4, [60/−60/60/0SMA]s.).
3.3. The Effects of Length-to-Width Ratios
Figure 6 shows the flutter boundary versus fiber orientation [θ/−θ/θ/−θ]s under different length-to-width ratios of the panel. For the case when a/b≤2, increasing the fiber orientation from 0° to 90° results in decreasing the stiffness of the panel in the x direction, which will lower the flutter boundary of the panel. However, for the case when a/b>2, λcr increases as θ increases initially and decreases afterwards. In addition, Figure 6 reveals that an increase in the panel length-to-width ratios leads to a stiffer composite panel.
Flutter boundary versus layer angle with different length-to-width ratios.
3.4. The Effects of Thermal Loads
The curves of frequencies ω versus λ for the panel without SMA wires under different temperature variation are shown in Figure 7. It can be found in Figure 7 that as the dynamic pressures increasing the first and the second natural frequencies gradually approach each other and finally overlap, then the panel will be in a limit cycle oscillation condition. Also, the natural frequencies and λcr are reduced with the raise of the temperature variation. When the temperature variation of the panel is raised over the critical buckling temperature ΔTcr, the panel will be buckling (but dynamically stable) under small dynamic pressure and the first-order natural frequency will be zero. As the dynamic pressure increases up to the critical value λ1, the motion of the panel will turn to a flat (and stable) condition as shown in Figure 7. So it seems that the process for flutter is contrary to that for thermal buckling.
Natural frequencies versus dynamic pressure under different temperature variations.
The curves of frequencies ω versus λ for the panel embedded with/without SMA wires under different temperature variation are shown in Figure 8. The SMA wires have volume fractions of 1%, prestrain of 0.5%, and temperature of 50°C. It can be revealed from Figure 8 that the panel embedded with SMA wires has the same trend as the panel without SMA wires; however, frequencies are enhanced due to the recovery stresses caused by the SMA wires. Moreover, the SMA wires can suppress both the flutter and thermal buckling of the panel.
Natural frequencies versus dynamic pressure under different temperature variations and volume fraction of SMA wires.
Figure 9 shows stability margins for the panel embedded with/without SMA wires under combined thermal loads and dynamic pressure load. The SMA wires have volume of fractions of 1%, prestrain of 0.5%, and temperature of 50°C. The panel has four types of motion: in region (I), under small λ and ΔT, the panel is flat and stable; in region (II), for small λ and moderate ΔT, thermal buckling occurs; in region (III), for moderate λ and ΔT, LCO occurs; in region (IV), for sufficiently high λ and ΔT, chaotic motion occurs.
Stability margins.
Figures 10–14 present the time history responses of the panel at points A–E in Figure 9, respectively. Figure 10 shows that the transverse vibration of the panel corresponding to point A in district (I) decreases as time increases. Also, the vibration amplitude of the panel embedded with SMA wires converges quickly compared to the panel without SMA wires. The motions of points B and C in Figure 9 are depicted in Figures 11 and 12, respectively. In the case of low dynamic pressure, increasing the temperature variation of the panel above ΔTcr, the motion will be changed from stable flat to buckled condition. The panel without SMA wires is in a buckled state and the panel embedded with SMA wires is stable flat as shown in Figure 11. Also the buckling deflection of the smart laminated panel is smaller than the conventional laminated panel as observed from Figure 12. Therefore, the SMA wires can suppress the thermal buckling of the panel and reduce the buckling deflection significantly for a given thermal load.
Time history for λ=70 and ΔT=0.5ΔTcr.
Time history for λ=20 and ΔT=1.1ΔTcr.
Time history for λ=20 and ΔT=1.4ΔTcr.
Time history for λ=90 and ΔT=0.5ΔTcr.
Time history for λ=90 and ΔT=1.4ΔTcr.
The panel is stable flat at a higher dynamic pressure, and it will flutter as the temperature increased. Figures 13 and 14 plot the motions of points D and E in Figure 9. The panel embedded with SMA wires becomes convergence when the panel without SMA wires has LCOs as shown in Figure 13. Also, the amplitude of LCOs of the smart laminated panel is smaller than the conventional laminated panel as observed from Figure 14. As a consequence, the SMA wires can suppress the flutter of the panel and the amplitude of the LCOs can be significantly reduced for a given dynamic pressure.
3.5. The Influences of SMA Wires Temperature, Prestrain, and Volume Fraction
The influences of SMA wires temperature and volume fraction on λcr are depicted in Figure 15. The SMA wires have prestrain of 0.7% and the panel temperature variation is assumed to be zero. When the SMA wires temperature is higher than As, phase transformation from martensite to austenite occurs. During this transformation, SMA wires can generate large recovery stress until the temperature is higher than Af. Thus λcr is enhanced via raising the SMA wires temperature as shown in Figure 15. Moreover, as shown in Figure 15, when SMA wires temperature is higher than As, increasing the SMA wires volume fraction leads to an improvement on λcr. Thus, for the purpose of enhancing the load-carrying capacity of the smart laminated panel, the parameters of SMA wires must be chosen carefully.
Curves of λcr versus T with various Vs.
The heated SMA wires can generate recovery stresses and then produce additional stiffness that will change the dynamic response of the panel. SMA wires have prestrain of 0.7% and temperature of 55°C. Temperature variation of the panel is assumed to be 1.2ΔTcr. When the temperature variation is 1.2ΔTcr, the panel is buckling. With the increase of dynamic pressure, the panel will become stable from the buckling state, and then the flutter will happen as shown in Figures 16 and 17. It is shown from Figure 16 that raising the volume fraction is able to enhance the stability margin of the panel. Figure 17 demonstrates the buckling deflection and amplitude of LCO versus dynamic pressure with different volume fractions of SMA wires. As displayed, using SMA wires can reduce both the buckling deflection and the amplitude of LCO. The results clearly indicate that the recovery stress introduced by SMA wires leads to a more stiffened panel for a wide range of dynamic pressure and thus a higher critical flutter dynamic pressure. Specifically, for a higher volume fraction as Vs=0.03, the buckling of the panel will not happen for a temperature variation 1.2ΔTcr.
Influence of Vs on stability boundary.
Influence of Vs on buckling, flat, and flutter phenomena.
Figures 18 and 19 reveal the influences of SMA wires temperature on ΔTcr, λcr, and transverse vibration of the panel. The temperature variation of the panel is assumed to be 1.2ΔTcr and the SMA wires have prestrain of 1% and volume fraction of 0.01. As shown in Figure 18, ΔTcr and λcr increase with raising the SMA wires temperature. Also, as displayed in Figure 19, the amplitude of LCO can be reduced by increasing the SMA wires temperature.
Influence of T on stability boundary.
Influence of T on buckling, flat, and flutter phenomena.
The influence of prestrain on the buckling and flutter behavior of the panel is displayed in Figures 20 and 21. The temperature variation of the panel is assumed to be 1.2ΔTcr and the SMA wires have volume fraction of 0.01 and temperature of 90°C. Raising the SMA wires prestrain can augment equal stiffness of the panel and therefore λcr and ΔTcr are increased as demonstrated in Figure 20. Figure 21 demonstrates the influence of prestrain on the buckling deflection and amplitude of LCO of the panel. It can be concluded that the buckling deflection and amplitude of LCO can be alleviated by increasing the SMA wires prestrain.
Influence of ε0 on stability boundary.
Influence of ε0 on buckling, flat, and flutter phenomena.
4. Conclusions
Thermal buckling and flutter behaviors of a laminated composite panel embedded with SMA wires subjected to nonlinear aerodynamic loading and thermal load have been analyzed in this paper. The von-Karman large deflection plate theory for structures, one-dimensional Brinson model for SMA wires, and the nonlinear piston theory for aerodynamics are used to derive the nonlinear governing equation of motion for the panel embedded with SMA wires. The system discrete dynamic model is obtained by employing the Galerkin method. A composite panel with a set of typical material constants and geometrical parameters is taken as an example to illustrate the method proposed in this study. The Runge-Kutta method has been employed to solve the system. The numerical results show the following:
When a/b is lower than 2, the aeroelastic stability of laminated panel is decreasing with the increase of the ply angle [θ/−θ/θ/−θ]s; however, when a/b is more than 2, the aeroelastic stability will increase as the ply angle increases initially and decrease afterwards.
Embedding SMA wires in the layers of composite panel can improve the aeroelastic stability boundary, and the most efficient design is to embed SMA wires in the airflow direction.
It is more significant to embed SMA wires in the outer layer of the laminated panel than in the inner layer for the flutter characteristics.
The critical flutter dynamic pressure and critical thermal buckling temperature can be enhanced by heating the SMA wires, increasing SMA wires volume fraction or prestrain. Therefore, the critical flutter dynamic pressure and critical thermal buckling temperature of the smart laminated panel can be greatly increased and the amplitude of LCO can be significantly reduced for a given flutter dynamic pressure.
The theoretical results presented in this paper can be applied in the practical engineering problem involving shape memory alloy wires. It is helpful for the aeroelastic analysis and vibration control of supersonic structures.
Nomenclature(A) Composite Parametersa:
Length of the panel
b:
Width of the panel
D110:
Value of the mass moment of inertia of the panel when all fibers are aligned with the x-axis
E1m,E2m:
Young modulus of matrix in 1 and 2 directions
Es:
Young modulus of SMA wire
h:
Thickness of the composite panel
M:
Mach number
Mkl,Nkl:
Force resultants operators
q:
Dynamic pressure
Q:
Lamina stiffness matrix
u,v,w:
Displacements of the panel in the x,y,z directions
u0,v0,w0:
Displacements of the midplane
Vm:
Volume fractions of the matrix
W:
Transverse vibration amplitude
v∞:
Relative free airstream velocity
Δp:
Aerodynamic pressure
ΔT:
Panel temperature variation
δT:
Virtual kinetic energy
δU:
Virtual strains energy
δW:
Virtual work done by the aerodynamic pressure
φj:
Eigen values of a matrix A
λ:
Dimensionless dynamic (=2qa3/MD110)
μ:
The maximum real part of φ
ρm,ρs:
Density of the matrix and SMA wires
υ12m,υs:
Poisson ratio of matrix and SMA wires
ω:
Panel natural frequency.
(B) SMA Brinson Model Parametersσr:
Recovery stress of SMA
Θ:
Thermoelastic modulus
EM:
Martensite Young modulus
EA:
Austenite Young modulus
εL:
Maximum residual strain
ε0:
Prestrain
Mf:
Martensite finish temperature
Ms:
Martensite start temperature
As:
Austenite start temperature
Af:
Austenite finish temperature
cA,cM:
Stress influence coefficient
ξ:
Total martensite volume fraction
ξs:
Stress induced martensite volume fraction
ξs0:
Initial stress induced martensite volume fraction
ξT:
Temperature induced martensite volume fraction
ξT0:
Initial temperature induced martensite volume fraction.
Competing Interests
The authors declare that they have no competing interests.
Acknowledgments
The authors gratefully acknowledge the National Natural Science Foundation of China (Grant no. 91216106) for the financial support of this work.
DowellE. H.Nonlinear oscillations of a fluttering plate19664712671275DowellE. H.Panel flutter, a review of the aeroelastic stability of plates and shells19708338539910.2514/3.56802-s2.0-0014752051MeiC.Abdel-MotagalyK.ChenR.Review of nonlinear panel flutter at supersonic and hypersonic speeds1999521032133210.1115/1.3098919BirmanV.LibrescuL.Supersonic flutter of shear deformable laminated composite flat panels1990139226527510.1016/0022-460X(90)90887-62-s2.0-0025703081KouchakzadehM. A.RasekhM.HaddadpourH.Panel flutter analysis of general laminated composite plates201092122906291510.1016/j.compstruct.2010.05.0012-s2.0-77954144722ZhaoH.CaoD.A study on the aero-elastic flutter of stiffened laminated composite panel in the supersonic flow2013332194668467910.1016/j.jsv.2013.04.0062-s2.0-84879014188AbbasJ. F.IbrahimR. A.GibsonR. F.Nonlinear flutter of orthotropic composite panel under aerodynamic heating19933181478148810.2514/3.117982-s2.0-0027642910ShiauL.-C.KuoS.-Y.LiuY.-P.Aerothermoelastic analysis of composite laminated plates20129461982199010.1016/j.compstruct.2012.01.0222-s2.0-84862787143XieD.XuM.DaiH.DowellE. H.Proper orthogonal decomposition method for analysis of nonlinear panel flutter with thermal effects in supersonic flow201533726328310.1016/j.jsv.2014.10.0382-s2.0-84916623534XueD. Y.MeiC.Finite element nonlinear panel flutter with arbitrary temperatures in supersonic flow199331115416210.2514/3.11332ZBL0779.730752-s2.0-0027334457LiF.-M.SongZ.-G.Flutter and thermal buckling control for composite laminated panels in supersonic flow2013332225678569510.1016/j.jsv.2013.05.0322-s2.0-84880923027BirmanV.Review of mechanics of shape memory alloy structures1997501162964510.1115/1.31016742-s2.0-0031272492ParkJ.-S.KimJ.-H.MoonS.-H.Thermal post-buckling and flutter characteristics of composite plates embedded with shape memory alloy fibers200536862763610.1016/j.compositesb.2004.11.0072-s2.0-23944465524OstachowiczW.KrawczukM.ŻakA.Dynamics and buckling of a multilayer composite plate with embedded SMA wires2000481–316316710.1016/s0263-8223(99)00090-22-s2.0-0033392145BarzegariM. M.DardelM.FathiA.GhadimiM.Aeroelastic characteristics of cantilever wing with embedded shape memory alloys20127918920210.1016/j.actaastro.2012.04.0232-s2.0-84862218645AsadiH.BodaghiM.ShakeriM.AghdamM. M.An analytical approach for nonlinear vibration and thermal stability of shape memory alloy hybrid laminated composite beams20134245446810.1016/j.euromechsol.2013.07.0112-s2.0-84883758508KuoS.-Y.ShiauL.-C.LaiC.-H.Flutter of buckled shape memory alloy reinforced laminates201221303502010.1088/0964-1726/21/3/035020IbrahimH. H.TawfikM.NegmH. M.Thermal buckling and nonlinear flutter behavior of shape memory alloy hybrid composite plates201117332133310.1177/1077546309353368MR28087872-s2.0-79951937235ZhaoX.LeeY. Y.LiewK. M.Mechanical and thermal buckling analysis of functionally graded plates200990216117110.1016/j.compstruct.2009.03.0052-s2.0-67349258477ShiauL.-C.KuoS.-Y.ChenC.-Y.Thermal buckling behavior of composite laminated plates201092250851410.1016/j.compstruct.2009.08.0352-s2.0-70349829007ShiY.LeeR. Y. Y.MeiC.Thermal postbuckling of composite plates using the finite element modal coordinate method199922659561410.1080/0149573992807792-s2.0-0032624043MatsunagaH.Thermal buckling of cross-ply laminated composite and sandwich plates according to a global higher-order deformation theory200568443945410.1016/j.compstruct.2004.04.0102-s2.0-11144251833BirmanV.2011New York, NY, USASpringerBrinsonL. C.One-dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable19934222924210.1177/1045389x93004002132-s2.0-0027575167DowellE. H.Nonlinear oscillations of a fluttering plate. II19675101856186210.2514/3.4316XueD. Y.MeiC.Finite element nonlinear flutter and fatigue life of two-dimensional panels with temperature effects1993306993100010.2514/3.464442-s2.0-0027698762