Chonghui Shao 1 and Dengqing Cao 1 and Yuqian Xu 1 and Hai Zhao 2
Academic Editor:Linda L. Vahala
1, School of Astronautics, Harbin Institute of Technology, P.O. Box 137, Harbin 150001, China
2, Resource Department, Southwest China Research Institute of Electronic Equipment, Chengdu 610000, China
Received 25 July 2016; Accepted 25 September 2016
This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. 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.
Figure 1: Schematic illustration of a smart laminated panel.
[figure omitted; refer to PDF]
2.1. Structure Modeling
The classical plate theory is employed to describe the displacements in the panel; that is, [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF]
Based on Hooke's law, the constitutive equation of the nth layer of the laminated composite panel under thermal loads is [figure omitted; refer to PDF] 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] [figure omitted; refer to PDF] respectively. Here, Q11m , Q22m , Q12m , and Q66m are the stiffness coefficients which are defined as [figure omitted; refer to PDF] in which G12m , E1m , and E2m stand for the shear modulus and Young modulus and [upsilon]12m and [upsilon]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]: [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF]
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: [figure omitted; refer to PDF]
The two constants in (9) are [figure omitted; refer to PDF]
The phase transformation coefficient Ω(ξ) and the elastic modulus Es (ξ) have the following expressions: [figure omitted; refer to PDF]
The austenite start temperature Asσ and finish temperature Afσ under stress can be expressed, respectively, as [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] 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.
Figure 2: Recovery stress generated by SMA.
[figure omitted; refer to PDF]
2.3. Aerodynamic Pressure Modeling
The aerodynamic load, Δp, is described by the third-order piston theory as [figure omitted; refer to PDF] 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. [figure omitted; refer to PDF] where the variations of kinetic energy δT and virtual work δW as well as variation of strains energy δU are given by [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] where (I0 ,I2 )=∫-h/2h/2 ρ(1,z2 )dz and the force resultants operators are obtained as [figure omitted; refer to PDF]
The boundary conditions for a simply supported panel are [figure omitted; refer to PDF]
The displacements of the panel satisfying the boundary conditions are written as follows: [figure omitted; refer to PDF] where the mode shapes are taken as [figure omitted; refer to PDF]
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 [figure omitted; refer to PDF] Let [figure omitted; refer to PDF] and then (22) can be expressed as [figure omitted; refer to PDF] 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 [figure omitted; refer to PDF] where y0 and [straight phi]j stand for the eigenvector and eigenvalues of the matrix A. The natural frequencies of the panel can be obtained as [figure omitted; refer to PDF]
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 [figure omitted; refer to PDF]
In the present study, the laminated panel with width of 0.4 m, length of 0.5 m, and thickness of 0.0025 m 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.
Composite Lamina Properties [figure omitted; refer to PDF]
SMA Wires Properties [figure omitted; refer to PDF]
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.
Table 1: Critical buckling temperature of the isotropic panel.
| Zhao et al. [19] | Shiau et al. [20] | Shi et al. [21] | Present |
Δ T c r (°C) | 126.5 | 127.1 | 126.4 | 126.43 |
Table 2: Comparison of critical buckling temperature for composite panel.
Ply angle | Δ T c r (°C) | |
Matsunaga [22] | Present | |
[ 0 / 90 / 90 / 0 ] s | 6.8 | 6.79 |
[ 0/45/-45/90]s | 7.6 | 7.59 |
Figure 3: Comparison of limit cycle amplitudes.
[figure omitted; refer to PDF]
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.
Figure 4: Flutter boundary versus orientation angles of SMA wires.
[figure omitted; refer to PDF]
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.
Figure 5: 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 .).
[figure omitted; refer to PDF]
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.
Figure 6: Flutter boundary versus layer angle with different length-to-width ratios.
[figure omitted; refer to PDF]
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.
Figure 7: Natural frequencies versus dynamic pressure under different temperature variations.
[figure omitted; refer to PDF]
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.
Figure 8: Natural frequencies versus dynamic pressure under different temperature variations and volume fraction of SMA wires.
[figure omitted; refer to PDF]
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.
Figure 9: Stability margins.
[figure omitted; refer to PDF]
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.
Figure 10: Time history for λ=70 and ΔT=0.5ΔTcr .
[figure omitted; refer to PDF]
Figure 11: Time history for λ=20 and ΔT=1.1ΔTcr .
[figure omitted; refer to PDF]
Figure 12: Time history for λ=20 and ΔT=1.4ΔTcr .
[figure omitted; refer to PDF]
Figure 13: Time history for λ=90 and ΔT=0.5ΔTcr .
[figure omitted; refer to PDF]
Figure 14: Time history for λ=90 and ΔT=1.4ΔTcr .
[figure omitted; refer to PDF]
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.
Figure 15: Curves of λcr versus T with various Vs .
[figure omitted; refer to PDF]
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 .
Figure 16: Influence of Vs on stability boundary.
[figure omitted; refer to PDF]
Figure 17: Influence of Vs on buckling, flat, and flutter phenomena.
[figure omitted; refer to PDF]
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.
Figure 18: Influence of T on stability boundary.
[figure omitted; refer to PDF]
Figure 19: Influence of T on buckling, flat, and flutter phenomena.
[figure omitted; refer to PDF]
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.
Figure 20: Influence of [straight epsilon]0 on stability boundary.
[figure omitted; refer to PDF]
Figure 21: Influence of [straight epsilon]0 on buckling, flat, and flutter phenomena.
[figure omitted; refer to PDF]
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:
(1) 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.
(2) 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.
(3) 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.
(4) 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 Parameters
[figure omitted; refer to PDF] :
Length of the panel
[figure omitted; refer to PDF] :
Width of the panel
[figure omitted; refer to PDF] :
Value of the mass moment of inertia of the panel when all fibers are aligned with the [figure omitted; refer to PDF] -axis
[figure omitted; refer to PDF] :
Young modulus of matrix in 1 and 2 directions
[figure omitted; refer to PDF] :
Young modulus of SMA wire
[figure omitted; refer to PDF] :
Thickness of the composite panel
[figure omitted; refer to PDF] :
Mach number
[figure omitted; refer to PDF] :
Force resultants operators
[figure omitted; refer to PDF] :
Dynamic pressure
[figure omitted; refer to PDF] :
Lamina stiffness matrix
[figure omitted; refer to PDF] :
Displacements of the panel in the [figure omitted; refer to PDF] directions
[figure omitted; refer to PDF] :
Displacements of the midplane
[figure omitted; refer to PDF] :
Volume fractions of the matrix
[figure omitted; refer to PDF] :
Transverse vibration amplitude
[figure omitted; refer to PDF] :
Relative free airstream velocity
[figure omitted; refer to PDF] :
Aerodynamic pressure
[figure omitted; refer to PDF] :
Panel temperature variation
[figure omitted; refer to PDF] :
Virtual kinetic energy
[figure omitted; refer to PDF] :
Virtual strains energy
[figure omitted; refer to PDF] :
Virtual work done by the aerodynamic pressure
[figure omitted; refer to PDF] :
Eigen values of a matrix [figure omitted; refer to PDF]
[figure omitted; refer to PDF] :
Dimensionless dynamic (= [figure omitted; refer to PDF] )
[figure omitted; refer to PDF] :
The maximum real part of [straight phi]
[figure omitted; refer to PDF] :
Density of the matrix and SMA wires
[figure omitted; refer to PDF] :
Poisson ratio of matrix and SMA wires
[figure omitted; refer to PDF] :
Panel natural frequency.
(B) SMA Brinson Model Parameters
[figure omitted; refer to PDF] :
Recovery stress of SMA
[figure omitted; refer to PDF] :
Thermoelastic modulus
[figure omitted; refer to PDF] :
Martensite Young modulus
[figure omitted; refer to PDF] :
Austenite Young modulus
[figure omitted; refer to PDF] :
Maximum residual strain
[figure omitted; refer to PDF] :
Prestrain
[figure omitted; refer to PDF] :
Martensite finish temperature
[figure omitted; refer to PDF] :
Martensite start temperature
[figure omitted; refer to PDF] :
Austenite start temperature
[figure omitted; refer to PDF] :
Austenite finish temperature
[figure omitted; refer to PDF] :
Stress influence coefficient
[figure omitted; refer to PDF] :
Total martensite volume fraction
[figure omitted; refer to PDF] :
Stress induced martensite volume fraction
[figure omitted; refer to PDF] :
Initial stress induced martensite volume fraction
[figure omitted; refer to PDF] :
Temperature induced martensite volume fraction
[figure omitted; refer to PDF] :
Initial temperature induced martensite volume fraction.
Acknowledgments
The authors gratefully acknowledge the National Natural Science Foundation of China (Grant no. 91216106) for the financial support of this work.
[1] E. H. Dowell, "Nonlinear oscillations of a fluttering plate," The American Institute of Aeronautics and Astronautics , vol. 4, no. 7, pp. 1267-1275, 1966.
[2] E. H. Dowell, "Panel flutter, a review of the aeroelastic stability of plates and shells," AIAA Journal , vol. 8, no. 3, pp. 385-399, 1970.
[3] C. Mei, K. Abdel-Motagaly, R. Chen, "Review of nonlinear panel flutter at supersonic and hypersonic speeds," Applied Mechanics Reviews , vol. 52, no. 10, pp. 321-332, 1999.
[4] V. Birman, L. Librescu, "Supersonic flutter of shear deformable laminated composite flat panels," Journal of Sound and Vibration , vol. 139, no. 2, pp. 265-275, 1990.
[5] M. A. Kouchakzadeh, M. Rasekh, H. Haddadpour, "Panel flutter analysis of general laminated composite plates," Composite Structures , vol. 92, no. 12, pp. 2906-2915, 2010.
[6] H. Zhao, D. Cao, "A study on the aero-elastic flutter of stiffened laminated composite panel in the supersonic flow," Journal of Sound and Vibration , vol. 332, no. 19, pp. 4668-4679, 2013.
[7] J. F. Abbas, R. A. Ibrahim, R. F. Gibson, "Nonlinear flutter of orthotropic composite panel under aerodynamic heating," AIAA Journal , vol. 31, no. 8, pp. 1478-1488, 1993.
[8] L.-C. Shiau, S.-Y. Kuo, Y.-P. Liu, "Aerothermoelastic analysis of composite laminated plates," Composite Structures , vol. 94, no. 6, pp. 1982-1990, 2012.
[9] D. Xie, M. Xu, H. Dai, E. H. Dowell, "Proper orthogonal decomposition method for analysis of nonlinear panel flutter with thermal effects in supersonic flow," Journal of Sound and Vibration , vol. 337, pp. 263-283, 2015.
[10] D. Y. Xue, C. Mei, "Finite element nonlinear panel flutter with arbitrary temperatures in supersonic flow," AIAA journal , vol. 31, no. 1, pp. 154-162, 1993.
[11] F.-M. Li, Z.-G. Song, "Flutter and thermal buckling control for composite laminated panels in supersonic flow," Journal of Sound and Vibration , vol. 332, no. 22, pp. 5678-5695, 2013.
[12] V. Birman, "Review of mechanics of shape memory alloy structures," Applied Mechanics Reviews , vol. 50, no. 11, pp. 629-645, 1997.
[13] J.-S. Park, J.-H. Kim, S.-H. Moon, "Thermal post-buckling and flutter characteristics of composite plates embedded with shape memory alloy fibers," Composites Part B: Engineering , vol. 36, no. 8, pp. 627-636, 2005.
[14] W. Ostachowicz, M. Krawczuk, A. Zak, "Dynamics and buckling of a multilayer composite plate with embedded SMA wires," Composite Structures , vol. 48, no. 1-3, pp. 163-167, 2000.
[15] M. M. Barzegari, M. Dardel, A. Fathi, M. Ghadimi, "Aeroelastic characteristics of cantilever wing with embedded shape memory alloys," Acta Astronautica , vol. 79, pp. 189-202, 2012.
[16] H. Asadi, M. Bodaghi, M. Shakeri, M. M. Aghdam, "An analytical approach for nonlinear vibration and thermal stability of shape memory alloy hybrid laminated composite beams," European Journal of Mechanics A--Solids , vol. 42, pp. 454-468, 2013.
[17] S.-Y. Kuo, L.-C. Shiau, C.-H. Lai, "Flutter of buckled shape memory alloy reinforced laminates," Smart Materials and Structures , vol. 21, no. 3, 2012.
[18] H. H. Ibrahim, M. Tawfik, H. M. Negm, "Thermal buckling and nonlinear flutter behavior of shape memory alloy hybrid composite plates," Journal of Vibration and Control , vol. 17, no. 3, pp. 321-333, 2011.
[19] X. Zhao, Y. Y. Lee, K. M. Liew, "Mechanical and thermal buckling analysis of functionally graded plates," Composite Structures , vol. 90, no. 2, pp. 161-171, 2009.
[20] L.-C. Shiau, S.-Y. Kuo, C.-Y. Chen, "Thermal buckling behavior of composite laminated plates," Composite Structures , vol. 92, no. 2, pp. 508-514, 2010.
[21] Y. Shi, R. Y. Y. Lee, C. Mei, "Thermal postbuckling of composite plates using the finite element modal coordinate method," Journal of Thermal Stresses , vol. 22, no. 6, pp. 595-614, 1999.
[22] H. Matsunaga, "Thermal buckling of cross-ply laminated composite and sandwich plates according to a global higher-order deformation theory," Composite Structures , vol. 68, no. 4, pp. 439-454, 2005.
[23] V. Birman Plate Structures (Solid Mechanics and Its Applications) , Springer, New York, NY, USA, 2011.
[24] L. C. Brinson, "One-dimensional constitutive behavior of shape memory alloys: thermomechanical derivation with non-constant material functions and redefined martensite internal variable," Journal of Intelligent Material Systems and Structures , vol. 4, no. 2, pp. 229-242, 1993.
[25] E. H. Dowell, "Nonlinear oscillations of a fluttering plate. II," AIAA Journal , vol. 5, no. 10, pp. 1856-1862, 1967.
[26] D. Y. Xue, C. Mei, "Finite element nonlinear flutter and fatigue life of two-dimensional panels with temperature effects," Journal of Aircraft , vol. 30, no. 6, pp. 993-1000, 1993.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright © 2016 Chonghui Shao et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
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.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer